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EXECUTIVE  SUMMARY 


This  is  the  final  report  submitted  to  the  Naval  Research  Laboratory  which  contains 
the  results  and  conclusions  of  research  performed  under  grant  #N00014-91-J-2023 
between  September  of  1991  to  June  of  1995.  This  effort  benefited  from  close 
collaboration  between  the  primary  institution,  Worcester  Polytechnic  Institute,  and  the 
Factory  Mutual  Research  Corporation  as  a  subcontractor.  A  Ph.D.  student  was  employed 
full-time  during  91-94  period  and  another  graduate  student  and  a  number  of 
undergraduate  students  participated  in  the  effort. 

This  research  resulted  in  the  development  of  a  computer  model  for  prediction  of 
Horizontal  Flame  Spread  on  Compartment  Surfaces  (HFSCS)  complemented  by  a 
Combustion  Model.  The  HFSCS  was  developed  to  predict  flame  spread  on  real  material 
(whose  properties  are  not  always  known  or  well-defined).  The  intent  of  the  work  was  to 
provide  the  Navy  with  a  validated  and  practical  numerical  model  that  can  be  used  either  as 
an  stand  alone  for  material  flame  spread  characterization  or  ideally  be  incorporated  into  a 
compartment  fire  model,  such  as  C-FAST  developed  by  NIST. 

The  research  program  can  be  described  by  the  following  key  components: 

1.  Development  of  the  core  of  the  HFSCS  for  charring  and  non-charring  material  by 
relying  on  a  pragmatic  integral  model.  This  model  requires  key  material  flammability 
properties  and  char  thermal  properties  applied  to  a  transient  integral  heat  transfer  model 
with  surface  reradiation.  The  need  for  known  properties  was  minimized  through  careful 
theoretical  analysis. 

2.  Development  and  validation  of  a  methodology  to  obtain  the  key  properties  referred 
to  in  (1)  using  a  numerical  model  and  measurements  from  a  standard  Cone  Calorimeter 
test  method  (ASTM  1354)  with  (or  without)  additional  measurements. 

3.  In  order  for  HFSCS  to  be  applicable  to  real  material  and  be  practical  (i.e. 
executed  on  a  Unix  Workstation  as  well  as  a  PC  and  in  a  short  time)  one  could  not 
attempt  to  solve  a  complex  2  or  3  dimensional  heat  transfer  problem  (including  radiation) 
with  combustion,  surface  kinetics  and  pyrolysis.  We  opted  for  using  the  surface 
temperature  to  track  the  flame  front  position  (i.e.  when  the  surface  temperature  reached 
the  pyrolysis  temperature  we  assumed  the  flame  fi'ont  to  have  reached  that  point)  and  thus 
being  able  to  rely  on  a  thermal  model  which  uses  key  experimental  results.  The  analysis  of 
the  problem  revealed  that  the  flame  heat  fluxes  need  to  be  characterized,  especially  the 

gas-phase  convective  flux  at  the  leading  edge.  We  identified  the  maximum  flux,  q  g ,  and 
the  characteristic  thermal  length,  6g,  to  be  the  key  parameters.  These  two  parameters 
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account  for  the  complex  heat  transfer,  finite  chemistry  and,  flow  dynamics  at  the  flame 
leading  edge  as  well  as  the  ambient  conditions.  Thus,  we  developed  a  methodology  to 
obtain  these  two  parameters  using  flame  spread  data  obtained  fi'om  a  constant  speed 
flame  spread  experiment. 

4.  We  developed  a  Constant  Speed  Horizontal  Flame  Spread  (CSHFS) 

experiment  to  both  evaluate  the  HFSCS  predictions  and  construct  a  methodology  to 

•  " 

determine  the  two  key  parameters,  qg  and  6g,  described  in  (3). 

5.  An  integral  combustion  model  was  developed  and  validated.  This  model 
was  first  examined  for  a  jet  flame  (to  obtain  some  confidence)  and  then  applied  to  a  pool 
fire  configuration.  The  combustion  model  determines  the  flame  height,  major  species 
yield,  flame  temperature  and  velocity  distribution  and  most  importantly  (for  the  HFSCS) 
the  radiation  distribution  per  unit  height  of  the  flame.  This  radiation  distribution  can  be 
used  to  calculate  the  flame  radiative  flux  ahead  of  the  flame  fi'om  as  well  as  to  the  fuel 
burning  surface.  The  latter  with  the  developed  pyrolysis  model  can  be  used  to  calculate 
the  mass  pyrolysis  rate  and  thus  the  heat  release  rate.  However,  calculating  the  flame 
radiation  surface  to  the  burning  surface  is  not  trivial  as  the  soot  absorption  needs  to  be 
accounted  for.  It  is  also  possible  to  use  the  flame  temperature  distribution  with  species 
and  soot  concentration  predictions  to  calculate  radiation  fluxes  using  ray-tracing.  The 
combustion  model  provides  the  important  coupling  needed  between  the  flame  spread 
model  and  the  increasing  burning  area  as  the  flame  spreads.  We  found  the  model  to  be 
numerically  intensive,  unless  the  centerline  flame  velocity  calculation  could  be  avoided  by 
using  an  empirical  relationship.  The  combustion  model  uses  state  relationship  for  species 
concentrations  that  need  to  be  obtained  for  these  real  burning  material.  This  is  an  area  of 
major  future  research.  However,  one  can  use  empirical  results  and  data  to  execute 
HFSCS  without  coupling  to  the  combustion  model.  At  this  time,  this  coupling  has  not 
been  implemented  because  of  the  numerical  intensity  of  the  combustion  model,  lack  of 
data  for  state  relationships  for  real  surfaces  and  difficulties  in  modelling  the  radiative  heat 
flux  feed-back  to  the  burning  surface. 

The  EDFSCS  has  been  validated  using  a  constant  speed  flame  spread  experiments 
that  simulated  a  radiating  flame  by  an  external  radiant  panel.  The  flame  spread  model 
developed  here  has  the  following  major  segments: 

1)  An  integral  heat-up  and  pyrolysis  model,  which  has  been  previously  developed  and 
validated  both  analytically  and  experimentally  for  charring  and  non-charring 
material,  as  the  core  component  of  the  horizontal  flame  spread  model.  The  flame 
spread  model  has  been  tested  for  sudden  increase  in  heat  flux  which  is 
representative  of  gas-phase  convective  heat  flux  and  a  new  formulation  has  been 
used  to  accomplish  this. 

The  integral  heat-up  model  uses  as  its  input,  the  material  flammability  properties 
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(e.g.  k,  p,  c,  L)  and  heat  fluxes  at  the  fuel  surface,  to  determine  surface  temperature  rise 
ahead  of  the  flame  front  from  which  the  flame  front  location  can  be  determined. 

2)  A  combustion  and  radiation  model  has  been  developed  and  validated  through 
limited  experimental  data. 

The  combustion  model  uses  a  direct  correlation  of  the  fluctuation  of  a  conserved 
scalar  and  state  relationship  for  species  concentration.  The  soot  formation  is  simulated  for 
an  optically  thin  flame  using  the  smoke-point  heat  release. 

The  constant  speed  flame  spread  experiment  was  designed  to  provide  measurement 
of  forward  gas-phase  convective  heat  flux  which  is  an  essential  mechanism  for  creeping 
flame  spread.  In  this  experiment,  the  flame  and  pyrolysis  fronts  were  stationary  as  the  fuel 
bed  was  moved  at  a  constant  speed  by  a  slide  towards  the  flame.  Surface  and  gas 
temperatures  as  well  as  heat  flux  histories  were  measured  at  three  stations  over  samples  of 

PMMA  and  particle  board  as  they  approached  the  stationary  pyrolysis  front  in  turn.  An 

•  • 

exponential  decay  type  gas-phase  convective  heat  flux  profile  (maximum:  ^gand  length 

scale:  5g)  has  been  derived  from  the  experiment.  The  results  are:  for  particle  board,  qg 

w30  kW/m^  and  5g«l  mm;  for  PMMA,  »70  kW/m^  and  5g«1.4  mm.  The  uncertainty  of 
these  results  from  a  representative  case  is  about  14%. 

The  horizontal  flame  spread  code  has  been  assembled  including  burning  and 
radiation  from  the  spreading  fire.  The  numerical  model  was  validated  by  comparison 
predicted  flame  spread  rates  and  surface  temperatures  with  measurements  from  several 
experiments  which  included  effects  of  external  and  flame  radiative  heat  fluxes.  Good 
agreements  between  the  numerical  predictions  and  the  experimental  data  were  obtained. 

The  theoretical  details  of  the  models,  subroutines  and  analysis  are  contained  in 
Appendixes  A-F.  The  computer  code  for  the  HFSCS  and  the  combustion  model  are 
provided  on  magnetic  disk.  A  listing  of  the  HFSCS  and  the  combustion  model  along  with 
their  input  and  description  of  the  input  are  included  in  Appendix  (G-J).  Extensive 
presentation  of  the  data  and  results  are  included  in  the  Appendixes  L  and  M. 

This  work  has  resulted  in  a  number  of  refereed  publications  and  other 
presentations  during  1991-1995  which  are  listed  here  for  completeness.  Additional 
paper(s)  have  and  will  be  submitted  for  publication. 
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Chapter  1  Introduction  and  Background 

Significant  progress  has  been  made  in  enclosure  fire  modeling  over  the  last  decade. 
Several  zone  models  have  been  developed  (e.g.  CFAST  (Jones,  1990),  BRI2  (Tanaka,  1983)). 
Though  these  fire  models  generally  predict  the  environment  generated  fi'om  fires  (i.e. 
temperature  and  velocity  of  fire  induced  flow)  and  transport  of  fire  products,  they  rely  on 
empirical  correlations  or  experimental  data  for  heat  release  rates  which  are  inputs  to  such 
models  because  predictive  capability  of  fire  spread  and  fire  growth  is  limited.  The  use  of  heat 
release  rates  fi'om  specific  test  (e.g.  furniture  calorimeter)  in  the  fire  models  to  simulate  fire 
growth  histories  in  real  fire  scenarios  having  various  compartment  geometries  and  ambient 
conditions  makes  it  (using  prescribed  heat  release  rates)  questionable. 

One  would  argue  that  a  robust  model  should  be  able  to  account  for  ambient  conditions 
(e.g.  heat  flux  from  the  environment),  predict  fire  spread  and  growth  and  heat  release  rate  for 
an  arbitrary  material  (charring  and  non-charring),  and  for  an  arbitrary  burning  scenario 
(upward,  downward,  horizontal  flame  spread,  etc.).  Such  a  flame  spread  model  should  also 
have  the  capability  of  estimating  products  from  the  burning  of  pyrolysis  gases  (CO,  COj,  etc). 
The  ability  to  predict  fire  growth  and  consequently  smoke  transport  together  with  appropriate 
engineering  methods  can  potentially  lead  to  improved  ship  designs  and  fire  control/fighting 
strategies. 

Upward  (vertical)  flame  spread  have  been  studied  extensively  and  several  numerical 
simulation  models  have  been  developed  (Delichatsios,  M.M.  etal,  1990;  1991;Mitler,  1990). 
These  upward  flame  spread  models  can  match  the  needs  of  current  compartment  fire  models. 
The  focus  in  this  research  is  opposed  flow  horizontal  flame  spread  although  a  lot  of  work  has 
been  done  in  this  area.  Numerical  models  of  horizontal  flame  spread  are  either  too  complex 
to  use  (di  Blasi,  1987,  Chen,  1990)  because  their  models  require  large  computational  time  and 
assume  kinetics  for  gas  and  solid  phases  or  are  not  general  enough  (Atreya,  1984,  Quintiere, 
1981).  The  latter  approaches  have  merits  for  practical  usage  since  they  do  not  require  large 
computational  time  as  the  former.  In  fact,  most  widely  used  method  to  quantify  horizontal 
(lateral)  flame  spread  is  described  in  ASTM  1354  (Quintiere,  1981;  Quintiere  et  al.,  1984). 
This  method  is  not  general  because  1)  it  is  suitable  for  thermally  thick  material,  but  not  valid 
for  thermally  thin  and  intermediate  materials,  2)  it  can  be  used  for  ambient  conditions,  but  not 
apply  to  other  environmental  conditions.  It  is  the  goal  of  this  project  to  develop  a  generalized 
model  that  uses  material  thermal  properties  which  can  be  obtained  from  small  scale 
experiment  for  predicting  horizontal  flame  spread. 

Horizontal  flame  spread  is  a  complex  problem  and  requires  modeling  of  a  number  of 
phenomena  (Atreya,  1984),  namely  1),  heat  transfer  inside  the  fuel  including  pyrolysis  of  the 
fuel,  2),  combustion  of  the  pyrolyzing  gases  and  heat  transfer  from  the  flame  to  the  fuel 
surface  (radiative  and  convective).  In  light  of  this,  several  numerical  models  are  needed  to 
simulate  horizontal  flame  spread: 
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1.  A  model  for  prediction  of  material  pyrolysis  during  flaming  combustion  for  both 
charring  and  non-charring  materials. 

2.  A  combustion  model  that  couples  the  pyrolysis  with  burning  in  the  gas  phase. 

3.  A  model  for  the  preheating  of  the  unbumed  material  prior  to  flame  spread.  This 
preheating  must  be  modeled  using  a  complete  heat  transfer  treatment  that  includes  all 
three  modes  of  heat  transfer:  conduction,  convection  and  radiation,  including  surface 
reradiation. 

4.  A  model  that  quantifies  gas  phase  conductive  heat  flux  profile  using  experimental  data 
obtained  fi’om  a  constant  speed  horizontal  flame  spread  (CSHFS)  experiment.  This 
profile  would  take  into  consideration  of  ambient  flow  condition. 


Convective 


Figure  1-1  The  mechanism  of  horizontal  flame  spread 

Flame  spread  occurs  when  the  fuel  surface  is  heated  to  ignition  or  pyrolysis 
temperature.  Therefore,  surface  temperature  histories  have  to  be  evaluated.  In  order  to  do 
that,  heat  fluxes  at  the  fuel  surfece  have  to  be  determined.  There  are  a  number  of  heat  fluxes 
involved  in  heating  the  fuel  surface  ahead  of  flame  fi-ont  (see  Figure  1-1):  1)  external  heat 
flux(es),  that  comes  fi'om  other  sources  (such  as  hot  ceiling  layers)  and  other  flames,  2) 
radiative  heat  flux  fi’om  the  flame  to  the  fuel  surfece,  qy^j,  3)  natural  convection  loss  from  the 
fuel  surface,  q^,  4)  surface  reradiation  loss,  and  5)  close  to  the  flame  fi'ont,  there  is  an 
additional  heat  flux  fi’om  the  flame  to  the  fUel  surface,  qy^n^,  which  is  called  gas-phase 
convective  (conductive)  heat  flux. 
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A  discussion  of  each  heat  flux  and  its  status  concerning  current  work  is  followed  in  the  next 
paragraphs. 

1 .  A  integral  heat-up  and  pyrofysis  model,  that  uses  material  flammability  properties  and 
heat  fluxes  at  the  fuel  surfece  described  above  has  been  developed  and  validated  both 
analytically  and  experimentally  (Chen  et  al.,  1991,  1993a).  This  model  calculates  fuel 
surface  temperature  and  pyrolysis  rate  histories  by  including  convective  heat  loss  and 
surface  reradiation,  q"  and  q"„.  Natural  convection  loss  is  treated  using  an  empirical 
model  for  natural  convection  over  a  horizontal  plate  which  can  be  found  in  the 
literature  readily.  The  surface  reradiation  loss  can  be  simply  determined  as; 

(1-1) 


where  e  is  surface  emissivity  and  a  is  Stefan-Boltzmann  constant,  T,  and  To 
are  the  surface  and  ambient  temperatures  respectively.  The  integral  heat-up 
model  was  developed  earlier  (Chen  et  al.,  1991,  1993a)  and  was  not  in  the 
scope  of  the  current  work  and  it  will  not  be  discussed  further  here.  This 
model,  has  been  incorporated  into  a  comprehensive  horizontal  flame  spread 
model  in  current  work.  This  incorporation  requires  the  use  of  the  heat-up  and 
pyrofysis  model  in  discrete,  time  dependent  manner.  To  accomplish  this,  the 
heat-up  model  had  to  be  to  able  to  account  for  the  sudden  increase  in  the  heat 
flux  at  the  flame  front  and  transition  to  the  pyrolysis  model. 

2.  The  external  heat  flux(es)  towards  the  fuel  surface,  q'ex,,  can  be  evaluated  fi'om  a 
compartment  fire  model  or  measurements  and  it  is  not  the  concern  here.  The  model 
developed  in  this  work  simply  uses  it  as  an  input. 

3.  The  radiative  heat  flux  fi'om  the  flame  to  the  fUel  surface,  can  be  determined  by 
a  combustion,  radiation  and  heat  flux  model.  A  combustion  and  radiation  model  has 
been  developed  in  this  work  and  part  of  it  has  been  published  (Chen  et  al.,  1993, 
Delichatsios  et  al.,  1994).  This  model  uses  the  pyrolysis  rate  determined  by  the 
P5n'olysis  model  as  its  input  to  estimate  combustion  of  the  pyrolysis  gases  and 
radiation  distribution  along  the  flame  length. 

4.  Another  important  component  of  horizontal  flame  spread  is  forward  gas-phase 
conduction  (convection),  q"o„v  fi-om  the  flame  to  the  foel  surface.  Because  of  the 
complex  flow  involved  near  the  flame  fi-ont,  this  heat  flux  is  difficult  to  simulate  or 
measure.  Therefore,  a  quasi-steady  horizontal  flame  spread  experiment  has  been 
designed  and  conducted  to  determine  this  gas-phase  conductive  heat  flux. 

A  comprehensive  model,  that  incorporates  all  these  components,  has  been  developed 
in  this  research.  Numerical  results  have  been  compared  with  experimental  data  fi'om  several 
tests  and  good  agreement  is  observed.  The  model  has  been  documented  and  has  output  that 
can  be  readily  incorporated  into  compartment  fire  models  such  as  C-FAST  (Jones,  1990). 
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This  model  has  been  structured  such  that  its’  future  incorporation  into  C-FAST  would  be 
possible  with  relative  ease. 

In  this  report,  a  literature  review  of  opposed  flow  horizontal  flame  spread  research 
and  efforts  to  model  diffusion  flames  (including  fires)  and  their  radiation  are  presented  in 
Ch^ter  2.  The  physics,  numerical  model  as  well  as  comparison  between  the  numerical  model 
and  an  analytical  solution  is  discussed  in  Chapter  3.  Chapter  4  describes  formulation,  coding 
and  validation  of  an  integral  combustion  model.  A  description  of  the  Constant  Speed 
Horizontal  Flame  Spread  (CSHFS)  experiment  is  included  in  Chapter  5.  Chapter  6  contains; 
1)  analysis  of  temperature  and  heat  flux  measurement,  2)  derivation  of  gas  phase  conductive 
heat  flux,  and  3)  comparison  of  numerical  results  with  the  CSHFS  experiments.  Finally,  the 
concluding  remarks  and  fiiture  work  is  presented  in  Chapter  7.  A  number  of  appendixes  to 
this  report  contain  the  details  of  the  combustion  mode,  analytical  solution  for  a  constant 
horizontal  flame  spread  velocity,  details  of  the  interpretation  of  the  experimental  results, 
numerical  schemes  to  account  for  sudden  flux  changes  (incident  on  the  fuel  surface),  detailed 
comparison  of  numerical  and  experimental  data,  supporting  data  and  the  computer  code. 
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Chapter  2  Previous  Work  and  Literature  Review 


2.1  AN  OVERVIEW  OF  HORIZONTAL  FLAME  SPREAD  RESEARCH 


Modes  of  flame  spread  can  be  characterized  by  flame  spreading  over  the  surface  of 
a  combustible  material  in  a  gaseous  oxidizer  that  flows  (either  naturally  induced  or  forced) 
in  a  direction  opposite  or  concurrent  to  that  of  the  flame  propagation  direction  (Femandez- 
Pello,  1984).  Opposed  flame  spread  includes  horizontal,  downward  and  lateral  flame  spread 
while  concurrent  mode  of  flame  spread  includes  wind-assisted  and  upward  flame  spreads. 
This  review  work  is  limited  to  opposed  flow  flame  spread. 

2.1.1  Theoretical  Work 

Modem  work  in  horizontal  flame  spread  started  with  the  work  of  de  Ris  (1969)  where 
previous  references  on  flame  spread  are  cited.  He  pointed  out  that  Emmons  (1965)  discussed 
the  general  subject  of  fire  spread  and  Tarifa  and  Torralbo  (1967)  solved  solid  phase 
conduction  equations  with  heat  flux  from  gas  phase  to  solid  surface  as  boundary  conditions, 
de  Ris  (1969)  investigated  flame  spread  over  solid  fuels  in  a  uniform  flow  of  velocity,  U,  in 
the  direction  opposite  the  flame  spread.  He  assumed  infinite  Damkoler  number,  Da  (the  ratio 
of  a  characteristic  chemical  reaction  rate  to  a  characteristic  convective  or  diffusive  transport 
rate),  constant  thermodynamic  and  transport  properties,  and  equal  diffusivities  of  thermal 
energy,  fuel,  and  oxidant  (i.e.  Lewis  number  (Le)  of  unity).  He  derived  analytical  solutions 
for  flame  spread  both  for  thermally  thin  (approximate  analytical  solution)  and  thermally  thick 
(exact  analytical  solution)  conditions.  Delichatsios  (1986)  also  derived  an  exact  analytical 
solution  of  creeping  flame  spread  over  a  thermally  thin  solid.  Other  anal3^cal  solutions 
become  also  available  (Wichman  and  Williams,  1983 ). 

Wichman  and  Williams  (1983)  assume  that  the  flame-sheet  is  located  at  the  fuel 
surface,  consider  non-unity  Lewis  number  and  solve  coupled  gas  and  solid  phase  equations. 
A  very  interesting  result  is  obtained:  for  Lewis  number  unity  the  analysis  produces  the  same 
flame  spread  formula  as  that  obtained  by  de  Ris  ( 1 969) .  They  accurately  pointed  out  that  de 
Ris  (1969),  by  linearizing  the  boundary  conditions  for  the  fuel  concentration  at  the  surface, 
actually  forces  the  flame-sheet  to  be  located  at  the  fuel  surface. 

Frey  and  Tien  (1979)  developed  a  comprehensive  model  of  opposed  flow  flame 
spread.  The  analysis  uses  the  Oseen  approximation  (constant  flow  velocity)  and  includes  a) 
an  Arrhenius  type  chemical  kinetics  in  the  gas  phase,  b)  transient  heating  and  c)  Arrhenius 
pyrolysis  kinetics.  Their  numerical  solutions  predict  gas  temperature  and  species 
concentration  profiles  as  well  as  the  existence  of  a  quenching  layer  near  the  fuel  surface  in  the 
region  of  the  inception  of  the  flame.  Another  interesting  aspect  of  this  model  is  that  it  is 
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capable  of  predicting  the  conditions  (flow  velocity,  pressure,  oxygen  concentration)  for  which 
flame  spread  stops. 

Other  full  numerical  solutions  are  also  available.  Chen  (1990)  solved  two-dimensional 
Navier-Stokes  momentum,  energy  and  species  equations  with  one-step  overall  chemical 
reaction  and  second-order,  finite-rate  Arrhenius  kinetics  in  gas  phase.  The  solid  phase  is 
thermally  thin  and  the  analysis  consists  of  an  energy  balance  coupled  with  the  heat  flux  from 
the  gas  phase  and  a  mass  balance  including  Arrhenius  pyrolysis  kinetics.  The  results  show 
that  flame  spread  speed  decreases  as  the  opposed  flow  velocity  is  increased.  His  model  can 
predict  flame  blowoff  limit  but  experimental  validation  is  not  complete.  Other  interesting 
phenomena  observed  from  his  analysis  are:  ahead  of  the  flame  front,  a  flow  recirculation  zone 
is  found  for  every  case  of  computation  and  the  flame  has  both  premixed  and  diffusion 
characteristics.  Di  Blasi  et  al.  (1987, 1989)  numerically  simulated  opposed  flow  flame  spread 
over  thin  paper  (charring  material).  Their  numerical  results  were  compared  with  experimental 
data  and  good  agreement  was  observed.  Bhattacharjee  and  Altenkirch  (1990)  numerically 
solved  opposed  flow  flame  spread  over  thin  solid  in  a  microgravity  environment.  Specifically, 
they  solved  gas-phase  radiation  using  three  simple  parameters  (a  Planck  mean  absorption 
coefficient,  the  fraction  of  total  emission  directed  towards  the  surface  and  a  shape  function 
that  represents  the  distribution  of  gas  to  surface  radiation).  At  low  opposing  flow  velocities, 
the  radiation  effects  become  significant.  Flame  spread  decreases  with  decreasing  opposed 
flow  velocity.  This  phenomenon  was  observed  in  drop-tower  experiments.  Recently,  West 
et  al.  (1994)  numerically  investigated  the  surface  radiation  effects  on  flame  spread  for 
thermally  thick  solid  in  an  opposed  flow.  They  claimed  that  the  ratio  of  the  rate  of  heat 
transfer  by  reradiation  from  the  surface  to  that  by  conduction  from  the  gas  to  the  solid  is 
proportional  to  the  length  over  which  heat  can  be  conducted  forward  of  the  flame  to  sustain 
spreading.  For  thick  solid,  this  length  decreases  with  increasing  flow  velocity  at  lower 
velocities  and  is  determined  by  gas-phase  process.  But  at  higher  velocities,  the  conduction 
length  is  determined  by  solid  phase  processes  and  is  independent  of  the  gas-phase  flow. 

Greenberg  and  Ronney  (1993)  studied  the  effect  of  Lewis  number  (Le)on  flame 
spread.  They  modified  de  Ris’  analysis  (Le  is  unity)  using  their  experimental  data  to  include 
the  influences  of  nonunity  Lewis  numbers  of  oxidant  and  fuel  vapors.  The  deduced  flame 
spread  relationship  shows  that  the  Lewis  number  of  the  oxidant  affects  the  spread  rate  only 
through  its  effect  on  the  flame  temperature  and  the  fuel  Lewis  number  plays  no  role.  They 
have  conpared  their  theoretical  results  with  the  flame  spread  data  over  thin  paper  samples  in 
a  variety  of  02-diluent  atmosphere  and  good  agreements  were  obtained.  They  concluded  that 
the  assumption  of  unity  Le  can  lead  to  rather  inaccurate  spread  rate  and  flame  temperature 
predictions. 

2.1.2  Experimental  Work 

De  Ris  (1969)  reviewed  McAlevy  and  Magee'  (1967)  work,  where  they  did  extensive 
experiments  for  thermally  thick  materials.  Atreya  ( 1984)  conducted  horizontal  flame  spread 
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experiments  for  wood  and  also  developed  a  numerical  simulation  model  to  predict  horizontal 
flame  spread.  This  model  included  aU  necessary  components  of  flame  spread  phenomena, 
namely,  combustion  of  pyrolyzing  gases,  conduction  and  pyrolysis  in  the  solid,  radiation  and 
convection  from  the  flame  to  the  fuel  surface.  He  verified  that:  1)  forward  radiation  heat 
transfer  from  the  flame  is  the  primary  accelerating  mechanism,  2)  forward  gas  phase  heat 
transfer  is  a  local  phenomenon  independent  of  fire  size  (only  material  dependent).  He  solved 
transient  partial  differential  equations  in  solid  phase  which  requires  large  computational  time. 

Quintiere  (1981, 1983  and  1988)  used  experimental  data  of  flame  spread  and  ignition 
from  the  LEFT  jqjparatus  to  generate  material  properties  data  to  calculate  lateral  flame  spread 
rate.  He  assumed  a  constant  heat  flux  from  the  flame  to  the  fuel  surface  ahead  of  pyrolysis 
front  extended  over  a  small  area.  The  gas  phase  conduction  from  the  flame  to  the  fuel  surface 
is  the  major  part  of  heat  flux  for  the  test  condition  of  the  LIFT  apparatus  because  flame  size 
is  usually  small.  An  analytical  solution  and  a  material  heat  flux  characteristic  property  for 
flame  spread  are  derived  for  calculating  flame  spread  speed. 

Femandez-Pello  and  Hirano  (1983)  reviewed  experimental  studies  on  opposed  and 
concurrent  flow  flame  spread.  They  realized  that  flame  spread  in  opposed  gas  flows 
frequently  occurs  at  near  extinction  or  non-propagating  condition  because  fires  usually  spread 
in  atmosphere  or  vitiated  ambient  oxygen  concentration.  The  oxygen  concentrations  are  close 
to  the  extinction  or  non-propagation  limit.  They  pointed  out  that  theoretical  description,  of 
the  flame  spread  process  based  solely  on  heat  transfer  analysis  are  inadequate  to  predict  rates 
of  spread  near  extinctions,  unless  steps  are  taken  to  empirically  introduce  the  effect  of  gas 
phase  chemical  kinetics  through  one  or  several  of  the  parameters.  Therefore  gas  phase 
chemical  kinetics  plays  a  critical  role  while  in  concurrent  flow  conditions,  heat  transfer  from 
the  flame  to  the  imbumt  fuel  is  the  primary  controlling  mechanism  of  flame  spread.  Ray  and 
Glassman  (1983)  further  investigated  the  variation  of  flame  propagation  rate  for  a  thermally 
thick  solid  as  a  function  of  opposed  flow  velocity.  They  concluded  from  their  analysis  and 
earlier  experimental  work  (Femandez-Pello  et  al.,  1980)  that  the  opposed  flow  flame  spread 
comprises  of  three  regimes:  at  low  opposed  flow  velocities  which  is  dominated  by  the 
naturally  induced  flow,  the  flame  spread  rate  shows  veiy  little  variation;  the  second  regime 
shows  a  linear  increase  of  propagation  rate  with  opposed  flow  and  is  dominated  by  thermal 
process  alone.  In  this  regime,  the  laminar  premixed  flame  spread  speed  at  the  lean  limit  is 
larger  than  opposed  flow  velocity.  The  flame  spread  speed  is  determined  solely  by  production 
of  fuel  vapor;  when  the  laminar  flame  spread  speed  is  less  than  opposed  flow  velocity,  as  in 
the  case  of  the  third  regime,  the  chemical  effect  dominates,  the  flame  spread  speed  declines. 
Opposed  gas  velocity  profile  has  significant  influence  on  the  flame  spread  (Di  Blasi  et  al., 
1989). 


Femandez-Pello  (1984)  reviewed  modeling  efforts  of  the  flame  spread  (both  opposed 
and  wind  assisted  flame  spreads).  He  pointed  out  that  most  of  the  early  work  of  flame  spread 
models  are  based  on  the  opposed  flow  mode  of  flame  spread  because  most  available 
experimental  information  was  for  this  configuration. 
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Ito  and  Kashiwagi  (1987)  conducted  PMMA  flame  spread  experiments  at  various 
angles  (-90  to  90°),  with  0°  being  horizontal  orientation.  Holographic  interferometry  was 
used  to  measure  the  temperature  distributions  along  the  sample  surface  (flame  spread 
direction)  and  towards  the  solid.  The  net  heat  flux  from  the  gas  phase  to  sample  surface  at 
the  pyrolysis  front  of  the  PMMA  sample  was  therefore  obtained.  They  concluded  that  the 
total  net  heat  transfer  rate  into  the  sample  from  the  gas  phase  is  about  56%  of  the  total  heat 
transfer  input  to  the  sample  at  0=-9O°,  78%  at  0=0°,  87%  at  0=+lO°,  and  99%  at  0=90°, 
where  -90°  is  downward  flame  spread  and  +90°  means  upward  flame  spread. 

Bhatanger  et  al.  (1990)  conducted  horizontal  flame  spread  experiments  of  insulation 
fibre  boards.  The  san5)les  were  moved  discretely  to  insure  that  the  flame  front  remains  at  the 
same  location  with  respect  to  the  external  radiating  surface.  A  constant  flame  spread  speed 
was  obtained,  effectively  eliminating  the  transient  effects  of  flame  spread.  They  also 
investigated  edge  effects  on  flame  spread  and  found  that  flame  spread  speeds  for  samples 
with  uninhibited  edges  (edges  are  not  wrapped  with  non-flammable  material)  are  higher  than 
that  of  samples  with  inhibited  edges.  Fangrat  and  Wolanski  (1991)  conducted  horizontal 
flame  spread  for  textile  floor  coverings  made  of  polypropylene  or  polyethylene).  The  flame 
spread  speed  as  a  function  of  external  heat  flux  was  obtained.  The  surface  and  bottom 
ten5)eratures  of  the  sample  were  measured.  Both  direct  and  schlieren  photographs  of  flame 
front  were  taken.  They  clearly  observed  from  their  schlieren  photographs  a  hot  zone  just 
ahead  of  the  flame  front.  This  zone  just  extends  a  comparatively  small  distance  ahead  of  the 
flame  front.  They  used  exponential  profile  for  radiative  and  convective  heat  fluxes  from  the 
flame  to  the  fuel  surface.  The  coefficient  of  these  heat  flux  distributions  were  first  fitted  with 
their  experimental  data.  A  general  numerical  solution  and  analytical  solution  for  a  thermally 
thin  material  (such  as  floor  coverings)  were  developed  to  predict  horizonal  flame  spread.  The 
result  gives  right  trends  though  the  predictions  do  not  agree  with  the  experiments  well 
because  they  did  not  accurately  measure  the  flame  heat  flux,  especially  the  gas-phase 
convective  heat  flux,  in  additional  to  over  simplicity  of  their  numerical  model. 

Zhang  et  al.  (1991)  experimentally  investigated  Lewis  number  effect  on  flame  spread. 
They  measured  flame  spread  rate  of  thin  paper  samples  and  flame  temperature  as  a  function 
of  diluent  type,  O2  mole  fraction,  pressure  and  fuel  bed  thickness.  Since  Lewis  number  has 
significant  influence  on  temperature  of  diffusion  flame,  the  flame  temperature  in  turn  affect 
the  rate  of  finite  chemistry,  consequently,  the  flame  spread  rate  would  be  altered  from  Lewis 
number  of  unity  as  described  in  de  Ris'  (1969)  work.  Experimental  data  of  Zhang  et  al. 
(1991)  show  that  the  flame  spread  rate  increases  with  increasing  oxygen  concentrations  for 
thin  paper  samples. 

Bhattachaijee  and  Altenkirch  (1992)  conducted  flame  spread  experiments  in  a 
microgravity  environment  using  a  sealed  chamber  filled  with  50/50  (volumetric) 
oxygen/nitrogen  mixture.  The  test  samples  were  thin  filter  papers.  Solid  and  gas-phase 
temperatures  were  measured.  The  experimental  results  show  that  the  fuel  preheat  zone 
upstream  of  the  gas-phase  flame  is  followed  by  a  constant  temperature  pyrolysis  zone,  which 
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is  followed  by  a  large  zone  of  surface  reactions.  They  have  also  compared  their  numerical 
results  (Bhattachaijee  and  Altenkirch,  1990)  with  the  experimental  data  and  suggested  that 
radiation  has  to  be  included  in  their  model  to  accurately  predict  flame  spread. 

2.2  COMBUSTION  AND  RADIATION  IN  FIRES 

An  important  component  for  modeling  and  estimating  fire  spread  and  fire  growth  is 
the  heat  flux  inposed  on  unpyrolyzed  and  pyrolyzing  materials.  Aside  from  heat  fluxes  from 
a  hot  upper  layer  or  other  sources,  heat  fluxes  originate  primarily  from  the  flame  under 
consideration  or  other  fires  adjacent  to  a  material.  Such  heat  fluxes  are  mainly  due  to 
luminous  soot  radiation  from  the  flames,  because  fires  in  general  produce  more  soot  than 
controlled  combustion  systems.  Thus,  it  is  essential  for  fire  spread  and  growth  predictions 
to  be  able  to  estimate  radiative  heat  fluxes  from  fires  produced  from  the  burning  of  the 
pyrolyzing  gases  of  a  given  material.  For  this  purpose,  predictive  methods  for  combustion 
in  fires  are  needed.  Such  predictive  methods  are  difficult  to  develop  because  fires  are 
turbulent  diffusion  flame  systems  dominated  by  buoyancy.  Despite  the  remarkable  progress 
in  turbulent  combustion  over  the  last  decade,  current  engineering  turbulent  models,  such  as 
k-e-g  models,  can  not  provide  reliable  tools  for  predicting  fluctuation  in  turbulent  buoyant 
combustion,  luminous  flame  radiation  and  products  of  combustion  issuing  from  the  burning 
of  an  arbitrary  material  (Bilger,  1976,  1977)  (note  that  k-e-g  models  have  been  shown  to  be 
adequate  to  account  for  the  dispersion  of  heat  and  combustion  products  after  combustion  is 
complete)  (Koylu  et  al.,  1990). 

Moreover,  application  of  detailed  k-e-g  models  in  compartment  fire  growth  and 
smoke  movement  codes  would  be  rather  prohibitive,  especially  for  multi-room  or  multi-story 
fires. 


Because  of  such  limitations,  i.e.  scientific  and  practical,  we  have  put  an  effort  in  this 
project  into  the  development  of  new  integral  type  models  for  turbulent  fires  (such  as  models 
by  Steward  (1970);  Tamanini  (1981);  Delichatsios  and  Mathews  (1989);  Cook  (1991)).  The 
major  innovation  in  this  new  integral  model  concerns  the  modeling  of  a)  turbulent  scalar 
fluctuations  and  b)  luminous  flame  radiation,  as  described  in  Chapter  4  and  Appendix  A. 
Briefly,  for  the  fluctuations  in  turbulent  buoyant  fires,  we  have  developed  a  general 
correlation  by  combining  analysis  and  experimental  data  (See  Appendix  B)  (note:  the 
magnitude  of  these  fluctuations  is  underpredicted  by  k-e-g  models  (Koylu  et  al.  (1990)).  For 
the  luminous  flame  radiation  in  turbulent  flames,  we  have  applied  a  soot  formation  model  ( 
Delichatsios,  1993)  that  has  been  derived  from  measurements  of  smoke-point  data  of  a  given 
fuel  material  (gaseous,  liquid  or  solid)  (See  Section  4.2).  It  should  be  pointed  out  that  the 
chemistry  of  combustion  is  reproduced,  at  it  is  commonly  done,  by  using  so-called  state 
relationships  between  chemical  species  and  mixture  fraction  (see  Section  4.2  and  Appendix 
A).  It  still  remains  a  challenge  to  obtain  these  state  relationships  from  measurements  of 
product  yields  from  practical  material  (such  measurements  can  be  done  in  a  hood  collector 
(Zukoski  et  al.,  1991,  Beyler,  1986)  or  in  a  turbulent  flow  (Tewarson  et  al.,  1993).  The 
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intended  results  of  the  integral  model  are  a)  flame  radiation  distribution  in  fires  and  b) 
products  of  combustion  including  soot  and  even  effects  of  vitiation  for  an  arbitrary  material. 
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Chapter  3 


Physics  of  Opposed-Flow  Horizontal  Flame 
Spread,  Description  of  the  Horizontal 
Flame  Spread  Model  and  Comparison  with 
an  Analytical  Solution 

3.1  PHYSICS  OF  HORIZONTAL  FLAME  SPREAD 

Flame  spread  occurs  when  fuel  surface  reaches  a  characteristic  temperature  at  which 
material  is  ignited  (Atreya,  1984;  Quintiere,  1981).  This  temperature  is  called  ignition  or 
pyrolysis  temperature  in  the  present  work  (it  is  assumed  that  pyrol5^is  temperature  is  the  same 
as  ignition  temperature  (Atreya,  1984;  Quintiere,  1981)).  A  study  of  heat  transfer  involved 
in  this  process  is  necessary  to  predict  fuel  surface  temperature  histories  and  flame  spread. 

There  are  two  types  of  flame  spreads,  namely  concurrent  and  opposed-flow  flame 
spreads.  Opposed-flow  flame  spreads  include  downward,  lateral  and  horizontal  flame  spreads 
(see  Chapter  2).  In  this  work,  we  focus  on  opposed-flow  flame  spread,  in  particular 
horizontal  flame  spread.  Heat  transport  mechanisms  in  opposed-flow  flame  spread  include 
(see  Figure  3-1);  convection/radiation  fi'om  the  flame  to  the  fuel  surface,  heat  loss  firom  the 
flame  to  the  ambient  environment  (by  radiation  and  conduction/convection), 
radiative/convective  heat  loss  from  the  fuel  surface  to  the  ambient,  conduction  normal  and 
parallel  to  the  fuel  surface.  The  relative  importance  of  the  heat  transfer  mechanisms  depends, 
among  other  things,  on  fuel  thickness,  ambient  oxygen  concentration  and  opposed  flow 
velocity.  Heat  transfer  to  the  imbumt  fuel  occurs  by  conduction  through  the  gas  fi'om  the 
flame  to  the  unbumt  region,  heat  conduction  from  behind  the  flame  through  the  solid  to  the 
solid  ahead  of  the  flame,  and  by  radiation  fi'om  the  flame.  Radiation  fi'om  the  flame  is 
relatively  unimportant  in  downward  flame  spread  because  the  view  factor  for  radiation  is 
small.  For  thermally  thin  solid,  heat  conduction  through  the  gas  is  generally  considered  as  the 
dominant  mode  of  heat  transfer  because  the  path  for  heat  transfer  through  the  solid  is  very 
limited  (VWlliams,  1976).  For  thick  solids,  both  heat  conduction  through  the  solid  phase  and 
gas  phase  are  important  (Femandez-Pello  and  Williams,  1975;  Femandez-Pello,  1984). 

Solid  phase  heat  transfer  is  important  for  flame  spread  over  a  solid  fuel  (di  Blasi  et  al., 
1989;  Ray  and  Classman,  1983).  Heat  is  conducted  at  directions  both  normal  and  parallel  to 
fuel  surface  and  by  radiation.  The  emissivity  of  the  fuel  surface,  the  material  properties  (k, 
p,  c,  L  (heat  of  pyrolysis)),  including  charring  and  boundary  conditions  at  the  back  surface, 
would  also  affect  the  solid  heat  transfer  and  flame  spread  rate  (Chen  et  al.,  1993a,  1995). 

Gas  phase  chemical  kinetics  and  gas  phase  heat  transfer  to  the  solid  surface  are  also 
important  factors  affecting  flame  spread.  Combustion  in  the  gas  phase  is  controlled  both  by 
mass  transfer  of  reactant  to  flame  as  well  as  by  chemical  kinetics.  Combustion  also  affects 
radiation  fi'om  the  flame  to  the  fiiel  surface,  flame  temperature  and  gas  phase  conduction 
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ahead  of  the  flame  front.  Chemical  kinetics  effects  depend  on  opposed  flow  velocity,  oxygen 
concentration,  pyrolysis  gas  composition  and  vapor  temperature.  It  is  difficult  to  determine 
these  effects  on  gas  phase  conduction  because  of  the  complex  flow  involved  near  the  flame 
front. 


Vf 


(1 )  Gas’phase  conduction/convection  (2)  Radiation  from  flame 

(3)  Reradiation  loss  at  the  fuel  surface  (4)  Conduction  into  the  solid 

(5)  Streamwise  Conduction  (6)  Vaporization  of  the  solid 

(7)  Opposed  flow 


Figure  3-1  Components  of  opposed  flow  horizontal  flame  spread 

Flame  spread  occurs  when  fuel  surface  reaches  pyrolysis  temperature,  Tp  (Atreya, 
1984;  Quintiere,  1981).  Therefore,  in  order  to  calculate  flame  spread  as  a  function  of  time, 
the  fuel  surfece  temperature  ahead  of  flame  front  has  to  be  determined.  The  solid  temperature 
increases  because  the  solid  fuel  receives  various  heat  fluxes  from  the  flame  and  other  sources 
as  shown  in  Figure  3-1 .  There  are  two  stages  involved  in  preheating  a  point  (the  preheating 
process  is  evaluated  by  a  transient  1-D  conduction  model)  on  the  fuel  surface  in  horizontal 
flame  spread:  1)  this  point  is  heated  by  external  heat  flux  from  other  sources  or  other  flames 
and  by  radiation  from  its  own  flame.  A  combustion  model  using  heat  release  rate  or  pyrolysis 
rate  (therefore  a  solid  pyrolysis  model)  is  thus  needed  to  calculate  the  radiation  from  its  flame 
2)  when  the  flame  front  approaches  the  point,  a  convective  heat  flux  is  applied.  This  heat 
flux  applies  to  a  very  short  distance  (O  (1  mm))  and  it  remains  to  be  determined.  We  can 
conclude  that  a  study  of  horizontal  flame  spread  should  include  the  following  components  as 
illustrated  in  Figure  3-1 : 

1 .  A  preheat  model  for  calculating  the  solid  temperature  rise, 

2.  A  pyrolysis  model  to  provide  pyrolysis  rate  from  different  materials. 
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3 .  A  combustion  model  coupled  with  the  solid  phase  to  compute  flame  radiation  to  the 
fuel  surface  and  coupling  to  the  flame  spread  phenomena  based  on  flame  growth 
which  is  affected  by  ambient  condition. 

4.  A  model  to  evaluate  gas-phase  conduction  near  the  flame  front  at  gas-solid  interface. 

Figure  3-2  illustrates  the  submodels  and  interconnections  of  the  horizontal  flame 
spread  model.  An  description  of  models  and  overview  of  work  involved  in  current  project 
are  presented  as  follows. 

A.  Heat-up  and  Pyrolysis  Model 

The  process  of  solid  temperature  rise  ahead  of  the  flame  is  predicted  by  an  integral 
heat-up  model  developed  earlier  (Chen,  1991,  Chen  et  al.,  1993a).  An  integral  pyrolysis 
model  is  also  used  to  calculate  pyrolysis  rate  when  the  fuel  surface  reaches  pyrolysis 
temperature,  Tp.  This  1-D  transient  conduction  model  was  also  developed  and  validated  in 
earlier  works  (Chen,  1991,  Chen  et  al.,  1993a)  and  has  been  developed  to  predict  pyrolysis 
of  charring  and  non-charring  materials.  The  integral  heat-up  and  pyrolysis  models  use  as 
inputs,  material  flammability  properties  such  as  k,  p,  c,  measured  in  standard  flammability 
apparatus  (Chen  et  al.,  1995,  Delichatsios  and  Saito,  1991).  Heat  fluxes  acting  on  the  fuel 
surface  are  needed  in  the  model  to  determine  temperature  rises  or  pyrolysis  rate  of  the  solid. 
The  pyrolysis  rate  is  required  input  for  the  combustion. 


Pyrolysis  Front 


Solid  Combustible  Surf. 


Figure  3-2  A  schematic  diagram  of  horizontal  flame  spread 
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There  are  several  heat  fluxes  involved  in  evaluating  surface  temperature  described  as 
follows;  1)  is  radiative  heat  flux  from  the  flame  originating  from  burning  of  the  pyrolysis 
gases.  This  radiative  heat  flux  is  quantified  by  a  radiation-combustion  model  which  uses 
pyrolysis  rate  given  by  the  pyrolysis  model  as  an  input.  2)  q"o^  is  called  gas-phase  conductive 
heat  flux.  This  heat  flux  drops  from  a  large  value  (~50  kW/m^)  to  0  over  a  very  short 
distance  (O  (Imm)).  Therefore,  it  is  very  difficult  to  measure  or  calculate.  A  constant  speed 
horizontal  flame  spread  experiment  has  been  designed  and  conducted  to  derive  this  heat  flux 
profile.  These  two  heat  fluxes  are  major  focuses  of  the  current  work.  The  following  three 
heat  fluxes  are  basically  inputs  to  the  horizontal  flame  spread  model  and  can  be  readily 
determined:  1)  q",^  is  external  heat  flux  originating  from  other  sources,  2)  4"  is  convective 
heat  loss  from  the  hot  fuel  surface  to  the  environment  which  can  be  calculated  using  an 
empirical  correlation  of  natural  convection,  3)  4'er  is  reradiation  loss  from  the  fuel  surface  and 
can  be  estimated  as; 


(3-1) 


where  T^  is  surface  temperature  and  T,,  is  ambient  temperature; 

B.  Combustion  Model 

The  pyrolysis  rate,  rh",  obtained  from  the  integral  pyrolysis  model  is  used  as  an  input 
to  an  integral  combustion  model.  The  combustion  model  uses  state  relationship  (laminar 
species  concentration  as  a  function  of  mixture  fraction),  a  correlation  to  account  for  turbulent 
fluctuations  of  the  scalar  quantities  together  with  equations  of  conservation  of  mass, 
momentum  and  energy  (in  terms  of  total  enthalpy)  to  solve  for  the  enthalpy,  species 
concentration,  centerline  velocity  and  temperature  at  discrete  intervals.  This  model  was 
validated  using  jet  flame  and  pool  fire  experimental  data  in  the  current  work.  It  simulates 
flame  radiation  and  its  relation  to  flame  turbulence  and  the  sooting  tendency  of  the  fuel. 
Chemical  reaction  and  combustion  is  modeled  using  conservation  equations  and  an  ensemble 
of  flamelets.  The  model  can  estimate,  as  a  function  of  the  height  from  the  flame  base, 
radiation  loss  distribution,  qj[(z),  radiation  fraction,  Xr  and  heat  release  rate.  Flame  length, 
tf,  and  major  species  production  rate  can  also  be  computed.  The  radiation  loss  distribution, 
heat  release  rate  and  flame  length  are  used  to  calculate  heat  flux  to  the  fuel  surface  ahead  of 
the  flame  front.  The  combustion  model  needs  to  be  coupled  with  the  flame  spread  model. 
Thus,  as  the  flame  grows,  the  energy  input  to  the  fiiel  surface  can  change  leading  to  a 
variation  ofm".  This  coupling  is  a  very  important  development  and  its  completion  requires 
availability  of  state  relations  for  real  fires  and  an  accurate  modelling  of  the  radiation  to  the 
fuel  surface.  At  the  present,  a  simple  coupling  using  known  correlations  for  mass  burning  rate 
and  heat  release  rates  (given  a  burning  area)  is  incorporated.  The  actual  predicton  of  mass 
loss  rate  depends  on  the  flame  radiation  feed  back  to  the  fuel  surface.  This  would  require 
additional  data  and  information  on  the  physics  of  the  soot  yield  and  flow  dynamics  near  the 
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base  of  the  flame.  Such  details  are  diflScult  to  model  and  the  research  in  this  field  is  still  in  a 
developing  stage.  The  integral  combustion  model ,  however,  is  discussed  in  detail  in  Chapter 
4  and  Appendix  A. 

C.  Heat  Flux  Model 


Using  the  radiation  loss  distribution,  the  heat  release  rate  at  height  intervals  and  the 
flame  length,  the  radiative  heat  flux,  q”^,  jfrom  the  flame  to  the  fuel  surface,  both  inside  and 
outside  the  flame,  can  be  detamined  by  a  heat  flux  model.  One  such  model  is  currently  being 
developed  by  Vezis  (1995)  in  our  research  group  and  the  final  model  can  be  incorporated  into 
the  horizontal  flame  spread  model.  This  model  will  not  be  discussed  here.  To  simulate  the 
radiation  ahead  of  the  flame,  we  have  used  an  exponentially  decaying  function  (discussed 
next),  which  is  a  very  good  representation  of  the  actual  phenomenon. 

D.  Constant  Speed  Horizontal  Flame  Spread  (CSHFS)  Experiment 


As  stated  before,  the  quantification  of  gas-phase  conduction  is  very  difficult.  An 
exponential  type  profile  for  the  gas-phase  conductive  (convective)  heat  flux  (this  can  be 
supported  by  the  numerical  analysis  (di  Blasi  et  al.,  1987,  1989)  and  analytical  solution 
(Atreya,  1984))  has  been  assumed: 


.  // 

^conv 


(3-2) 


where  x  is  distance  along  the  horizontal  surface  in  the  direction  of  flame  spread,  is  flame 
fi'ont  location,  6  g  is  length  scale  of  the  gas-phase  conductive  heat  flux  and  qg  is  the  magnitude 
of  this  heat  flux.  A  constant  speed  horizontal  flame  spread  experiment  was  developed  in  this 
work  in  order  to  obtain  6g  and  q^. 

In  the  CSHFS  e3q)eriment,  a  fuel  sample  was  moved  at  a  constant  speed  towards  the 
stationary  flame  or  pyrolysis  front.  The  flame  is  quenched  by  a  cooling  plate  to  maintain  a 
constant  (and  small)  burning  area  as  the  sample  was  moving.  In  this  way,  the  transient  effects 
of  flame  spread,  pyrolysis  rate  as  well  as  flame  radiation  can  be  eliminated.  Surface 
temperature,  total  heat  flux  histories  at  the  fuel  surface  and  flame  fi'ont  locations  were 
measured.  The  surface  temperature  histories  were  used  to  determine  the  gas-phase 
conductive  heat  flux.  The  description  of  the  experiment  is  presented  in  Chapter  5.  Chapter 
6  discusses  experimental  results  and  derivation  of  the  gas-phase  conduction.  The 
comparisons  for  both  flame  fi'ont  and  surface  temperature  between  experimental  data  and 
numerical  simulation  using  the  derived  gas-phase  conduction  are  also  shown  in  Chapter  6. 
More  results  can  also  be  found  in  Appendix  H. 
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By  using  the  heat-up  model  (Chen,  1993a)  (which  was  not  the  scope  of  the  current 
work,  but  has  been  incorporated  in  the  horizontal  flame  spread  model)  that  uses  material 
flammability  properties  and  heat  fluxes  described  above,  surface  temperatures  ahead  of  the 
flame  (i^rolysis  fi'ont)  can  be  obtained.  The  flame  fi'ont  location,  Xp  can  therefore  be  derived 
as  described  in  Section  3.2. 


3.2  DESCRIPnON  OF  THE  HORIZONTAL  FLAME  SPREAD  MODEL 

The  heat-up,  pyrolysis  and  combustion  models  which  was  briefly  described  here  have 
been  combined  to  form  a  comprehensive  horizontal  flame  spread  model.  This  model 
involves  several  steps  to  find  flame  fi'ont  location  and  pyrolysis  rate: 

1 .  Obtain  external  heat  fluxes  form  input  or  calculated  by  a  compartment  fire  model, 
calculate  flame  radiation  and  gas-phase  conductive  heat  fluxes  derived  based  on 
experiments  from  the  flame  to  each  node, 

2.  Using  above  heat  fluxes,  calculate  surface  temperature  for  each  node  using  either 
preheat  or  pyrolysis  model  depending  on  whether  the  node  is  at  heating  or  burning, 
sum  pyrolysis  rate  of  each  node, 

3.  For  the  pyrolysis  front,  both  heat-up  and  pyrolysis  would  be  calculated, 

4.  Find  new  pyrolysis  front  using  surface  temperature  data  obtained  in  steps  2  and  3, 

5.  The  sum  of  pyrolysis  rate  is  used  as  an  input  to  the  combustion-radiation  model  to 
evaluate  flame  radiation  to  the  fuel  surface,  then  go  back  to  step  1  for  the  next  time 
step. 

The  flow  chart  of  the  Horizontal  Flame  Spread  on  Compartment  Surface  (HFSCS) 
is  shown  in  Figure  3-3  where  the  interconnections  of  different  subroutines  are  illustrated.  The 
purpose  of  each  subroutine  is  described  as  follows: 


I. hfscs.f 
2.inputdata 

3. maincal 

4. flxte,  fixtf 

5.  preheat 

6. hpcg 

7. heatin 
S.heat 

9. qflux 

10.  vapor 

II. pyrobb 

12. pyrobv 

13.  front 
14.locxp 


The  name  of  the  HFSCS  code 

This  routine  is  used  to  read  the  physical  properties  (e.g.  k,  p,  c,  Tp,  L,  etc.)  of 
the  fuel  that  is/will  be  burning 

The  main  driver  calling  all  the  calculation  and  output  routines 
These  routines  obtain  external,  flame  heat  fluxes 
Driver  for  heat-up  calculation 

Hamming's  Predictor-Corrector  method  numerical  solver 
Heat-up  calculation  for  the  first  time  step 
Normal  heat-up  calculations 

Combine  all  the  heat  fluxes  absorbed  by  the  fuel  surface 

The  driver  for  transient  pyrolysis  process 

The  initial  pyrolysis  calculation 

Normal  pyrolysis  process  calculation 

Driver  to  call  routines  for  combustion 

Classifying  the  nodes  to  see  if  the  node  is  pyrolysis,  heat-up  or  the  section 
containing  the  node  just  starts  to  pyrolysis 
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1 5 .  getxp  Interpolation  to  obtain  pyrolysis  front 

Heat  release  rate  and  flame  height 
1 7.1octxb  Locate  burnout  front 

18. output  Output  results  (flame  front,  flame  height,  heat  release  rate,  etc.) 

This  program  has  been  successfully  run  on  Unix  and  PC  systems. 

3.3  LOCATING  THE  FLAME  (PYROLYSIS)  FRONT 

The  foel  is  divided  into  large  number  of  segments  with  n  nodes  along  the  fuel  length. 
The  surface  temperature  at  each  node  at  each  time  step  is  evaluated  using  the  preheat  or 
pyrolysis  model  (Chen  et  al.,  1993a). 

By  comparing  surface  temperature  with  the  pyrolysis  temperature,  the  flame  front 
location  can  be  bracketed  into  a  certain  segment.  This  is  called  the  pyrolysis  section.  The  first 
section  includes  the  fuel  segment  between  nodes  1  and  2,  with  node  1  as  the  beginning  of  the 
section  while  node  2  is  excluded  from  this  section.  Section  2  would  begin  with  node  2. 
Section  I  would  begin  with  node  I  and  end  with  I+l,  but  exclude  node  I+l).  The  pyrolysis 
front  is  located  between  nodes  Xj.,  and  X;  if  surface  temperatures  of  these  nodes,  T,^}.!  and  Tj  j 
are  higher  and  lower  than  pyrolysis  temperature,  respectively.  The  new  pyrolysis  front 
location  can  be  determined  either  from  the  two  nodal  surface  temperatures  or  from  the 
fictitious  temperature  of  the  pyrolysis  front,  T^p  and  Tj^j  depending  on  the  flame  front  location 
previously. 

To  find  new  pyrolysis  front  location,  Xp  one  has  to  determine  where  the  surface 
temperature  is  at  the  pyrolysis  temperature,  Tp ,  i.e. : 

9(<.'”)=0p 

where  0=Ts-To  and  0p=Tp-To,  Tq  is  initial  temperature.  is  pyrolysis  front  location  at  time 
t". 

If  X;  is  the  nodal  location  next  to  x^  where  pyrolysis  has  yet  to  start,  i.e.: 


at  time  t”^^’  one  calculates  by  using  the  heat-up  model,  the  fictitious  temperature  at  xj  (as  if 
pyrolysis  did  not  start)  and  surface  temperature  at  Xj,  if 

Q{x;r*^)>%cmdQ(x,r*^)<Q^ 
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(3-3) 


The  new  pyrolysis  front  can  be  obtained  by  using  a  linear  approximation: 
X  =X  +- - £ - ^ — (x-x) 

0(x;,r-‘)-0(x„r-i) 


In  reality,  the  process  of  finding  new  pyrolysis  front  location  is;  1)  the  model  figures 
out  whether  the  pyrolysis  front  entered  the  section  at  the  previous  time  step  and  2)  If  the 
pyrolysis  front  was  in  this  section  at  the  previous  time  step,  eqn.  (3-3)  is  used  to  locate  new 
pyrolysis  front,  otherwise,  the  fictitious  temperature  of  node  Xj.!  and  preheat  temperature  of 
node  jq  are  used  to  locate  the  new  pyrolysis  front,  eqn.  (3-3)  is  still  used,  with  jq.i  replacing 
xj  and  6(3^,  t"^*)  being  the  fictitious  temperature  at  node  I-l. 
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Flow  Chart  HFSCS 
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Figure  3-3  A  flowchart  and  interconnections  of  the  horizontal  flame  spread  on 
compartment  surfaces 
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3.4  SURFACE  TEMPERATURE  CALCULATIONS  USING  TWO  SOLID 
TEMPERATURE  PROFILES 


When  a  sudden  jump  of  heat  flux  (e.g.  gas  phase  conduction  from  flame)  starts  to 
apply  to  a  fuel  surface  which  has  been  heated  for  a  period  of  time  by  other  heat  fluxes  (e.g. 
external,  flame  radiative  heat  fluxes),  surface  tenperatures  of  the  fuel  and  flame  front  location 
can  not  be  accurately  calculated  and  estimated  by  the  flame  spread  code  developed  initially 
(see  also  Delichatsios,  M.M.  et  al.,  1990).  This  problem  arises  from  the  fact  that  there  exists 
actually  two  different  heat-up  layers:  1)  a  new  heat-up  layer  caused  by  the  sudden  jump  heat 
flux,  and  2)  the  original  heat-up  layer  caused  by  other  heat  fluxes.  The  original  model  uses 
a  single  tenperature  exponential  decay  type  profile  having  a  .single  length  scale.  This  would 
introduce  error  in  the  calculation  of  surface  temperature  as  well  as  thermal  penetration  length 
(Chen  et  al.,  1993a).  Therefore,  modification  of  the  original  code  was  made  by  observing 
that  the  heat  conduction  equation  used  in  the  integral  heat-up  model  (Chen  et  al.,  1993a): 


dd_  a^d 
at  d^z 


(3-4) 


can  be  split  into  two  equations  with  appropriate  initial  and  boundary  conditions  for  ts:ti,  when 
a  sudden  jump  heat  flux,  q^'-  q"  is  applied  (see  eqns.  (F-5a),  (F-5b),  (F-6a)  and  (F-6b)  of  the 
Appendix  F): 


S0J 

— Ua - 

dt  a\ 


(3-5a) 


(3-5b) 


802  8^02 

- =a - 

s<  a\ 


(3-5c) 


(3-5d) 
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(3-5e) 


@t=ti 


02(j;,ii)”O 


(3-5f) 


Therefore  two  tenperature  rises,  0i  and  02  can  be  evaluated  separately  by  using  two 
different  heat  conduction  equations  and  boundary  conditions.  In  other  words,  two 
temperature  rises  in  the  two  different  thermal  layers  caused  by  two  heat  fluxes  imposed  on 
fuel  surface  are  calculated  independently.  The  actual  temperature  of  the  solid  would  be: 

0-0, .02  (3.6) 


In  this  way,  the  correct  temperatures  can  be  determined.  The  original  code  has  been 
modified  to  account  for  this  situation  and  for  detail  see  Appendix  F.  Hereafter  aU  results 
from  the  model  are  from  the  modified  code  using  eqn.  (3-6)  and  the  described 
formation. 

Before  conducting  comparison  of  the  numerical  results  of  the  horizontal  flame  spread 
model  and  an  analytical  solution,  the  sensitivities  of  the  model  will  be  checked  at  first,  i.e.  the 
numerical  results  using  different  time  steps,  surface  nodal  spacings  and  heat  flux  profiles  will 
be  conqiared  and  this  will  be  presented  in  next  section.  The  accuracy  of  this  model  will  first 
be  checked  against  an  analytical  solution  which  will  be  presented  later  in  this  chapter,  and 
finally,  this  numerical  model  will  be  validated  by  comparing  the  numerical  results  with 
experimental  data  of  constant  speed  horizontal  flame  spread  (CSHFS)  which  we  will  be 
described  in  Chapter  6. 

3.5  SENSITIVITY  ANALYSIS  OF  THE  HORIZONTAL  FLAME  SPREAD 
MODEL 

Following  the  completion  of  the  numerical  code  for  horizontal  flame  spread,  a 
systematic  sensitivity  analysis  is  performed  by  varying:  1)  nodal  spacings  (uniform  mesh  sizes 
from  small  to  large  and  non-uniform  mesh  size  using  power  series),  2)  time  steps  (from  small 
to  large),  3)  heat  fluxes  (constant  external  heat  fluxes  of  different  levels;  heat  flux  distribution 
at  different  levels,  similar  to  the  experimental  conditions  which  will  be  described  in  Chapter 
5;  different  flame  heat  flux  profiles  (the  magnitudes  and  length  scale  of  convective  (gas-phase 
conductive)  heat  fluxes).  It  is  anticipated  that  the  limitation  of  the  numerical  code  will  be 
revealed  from  this  process. 
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Figures  3-4  through  3-6  fllustrates  the  effects  of  nodal  spacings,  time  steps,  and  heat 
fluxes.  The  accuracies  of  the  numerical  code,  needlessly  to  say,  depend  on  both  nodal  spacing 
and  time  step.  The  differences  between  different  nodal  spacings  are  clearly  manifested  in 
Figure  3-4.  Differences  between  results  of  Ax=0.5, 1,  2  and  4  mm  are  really  small.  Fairly 
large  difference  fi'om  Ax=10  mm  is  observed,  so  does  power  series  nodal  spacings.  It  can  be 
easily  concluded  if  nodal  spacing  is  shorter  than  the  shortest  length  scale  of  heat  fluxes  (in  this 
case,  6.67  mm  which  is  1/150  (m)),  the  results  would  not  differ  much.  Therefore  Ax=4  mm 
can  be  used  in  this  case  to  calculate  flame  spread. 

In  Figure  3-5,  the  results  from  different  time  steps  ranging  from  At=0.0625  to  At=l 
seconds  are  shown.  Figure  3-5  shows  substantial  difference  between  time  steps  of  1  and  0.5 
seconds.  On  the  other  hand,  the  difference  between  the  results  of  two  smallest  time  steps 
(At=0.0625  and  0.125  seconds)  are  fairly  small.  One  can  see  for  this  case,  that  the  numerical 
results  of  At=0.125  seconds  are  accurate  enough.  It  can  be  postulated  that  if  At=eAx/V, 
where  e  is  a  numbers  1  and  V  is  flame  spread  speed,  the  prediction  would  be  accurate.  The 
flame  spread  speed,  V,  in  this  case  is  about  5.3  mm/s,  therefore,  e=0. 125x5.3/4=0. 166.  This 
can  also  be  supported  by  Figure  3-7  which  will  be  discussed  later. 
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Figure  3-4  Sensitivity  analysis  of  the  horizontal  flame  spread  model-different  nodal 
spacings  and  At=  1  seconds 
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Finally,  the  results  from  different  heat  flux  profiles  are  demonstrated  in  Figure  3-6. 
Obviously,  increasing  heat  flux  and  heat  flux  profile  length  scale  would  increase  flame  spread 
speed.  It  has  to  be  pointed  out  that  we  use  the  external  heat  flux  distribution  experienced  in 
the  experiments  (see  Chapters  5  and  6)  to  obtain  two  sets  of  results  in  Figure  3-6.  One  may 
also  notice  that  in  Figure  3-1,  At=l  second  has  been  used  while  previous  results  indicate  that 
in  this  case,  At=0.125  seconds  should  be  used  in  order  to  obtain  accurate  results.  The 
rational  here  is  that  the  purpose  of  Figure  3-6  is  to  check  whether  the  code  would  give  right 
trend  as  one  changes  heat  fluxes,  rather  to  see  if  the  code  can  give  accurate  result  for  each 
case.  In  light  of  this,  At=l  second  is  sufficient  and  code  does  give  satisfactory  results. 

A  brief  conclusion  of  above  section  can  help  an  user  to  chose  correct  time  step  and 
nodal  spacing  in  order  to  obtain  accurate  result  from  the  horizontal  flame  spread  code.  The 
size  of  nodal  spacing.  Ax,  can  be  selected  as  no  larger  than  the  smallest  length  scale  of  heat 
fluxes  (usually  gas-phase  convective  heat  flux).  The  flame  spread  speed  can  be  guessed  next. 
A  fairly  conservative  estimate  of  flame  spread  speed  is  6  nun/s.  The  time  step.  At,  can  be 
chosen  as  At=0.166  AxfV.  These  time  step  and  nodal  spacing  are  then  used  in  the  flame 
spread  model.  The  flame  spread  speed  can  be  determined  and  compared  with  the  value  that 
is  guessed.  One  may  need  to  adjust  the  time  step  if  the  calculated  speed  is  larger  than  the 
guessed  flame  spread  speed.  Sensitivity  of  flame  spread  code  has  been  examined  for  a  similar 
problem  (Delichatsios,  M.M.,  1991).  It  can  be  concluded  from  their  analysis  that  the  stability 
requirement  is  less  stringent  than  the  accuracy  requirement  for  both  time  step  and  nodal 


Figure  3-5  Sensitivity  analysis  of  the  horizontal  flame  spread  model-different  time  steps, 
Ax=4mm 
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spacing  presented  above.  Therefore,  if  one  uses  time  step  and  nodal  spacing  that  are  derived 
above,  there  is  no  stability  problem. 


TIME  (sec) 


Figure  3-6  Sensitivity  analysis  of  the  horizontal  flame  spread  model-different 
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3.6  COMPARISON  WITH  AN  ANALYTICAL  SOLUTION 


In  order  to  validate  the  flame  spread  model,  the  comparison  between  numerical  results 
and  an  analytical  solution  using  Quintiere's  postulate  (see  Appendix  D)  is  carried  out.  For 
thermally  thick  solid  with  imposed  heat  flux,  q",  at  the  fuel  surface,  the  surface  temperature 
rise  of  the  solid  can  be  estimated  (Carslaw  and  Jaeger,  1959); 


T-T  = 


/- 


,// 


-dx 


(3-7) 


Therefore  for  a  point,  x,  on  the  fuel  surface  to  reach  pyrolysis  temperature,  Tp,  the 
time,  t,  in  eqn.  (3-7),  should  be  replaced  with  time  to  pyrolysis,  tp.  Eqn.  (3-7)  becomes: 


T-T  = 
-'p  ■'o 


Mp^o  ^|^ 


f— 


(3-8) 


It  is  assumed  that  the  imposed  heat  flux,  q",  consists  of  two  parts;  1)  an  external  heat 
flux,  q",  and  2)  a  flame  heat  flux  which  has  a  exponential  type  profile: 


.  // 


=Ae 


(3-9) 


One  solves  eqn.  (3-8)  together  with  the  heat  flux  profile  (eqn.  (3-9))  assuming 
constant  flame  spread  velocity,  V.  The  solution  of  eqn.  (3-8)  is  (eqn.  (D-7)  of  Appendix  D): 


T-T 

''o 


Ajblj? 

sjTikpc 


(3-10) 


where  Tj  is  surfece  temperature  rise  owing  to  external  heat  flux,  q".  It  can  be  shown  that  the 
time  to  pyrolysis  (Delichatsios  et  al.,  1991;  Quintiere,  1981): 


^  ^  kpc(T^-Tf 

^4 


(3-11) 
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Equation  (3-11)  can  be  substituted  into  eqn.  (3-10)  and  the  flame  spread  velocity,  V, 
can  be  solved  as: 


dt  4/^ 


(3-12) 


Equation  (3-12)  is  an  analytical  asymptotic  solution  for  constant  flame  spread  velocity 
which  is  also  referred  to  as  the  Quintiere's  postulate.  If  surface  reradiation  loss  is  taken  into 
account,  eqn.  (3-11)  can  be  modified  by  (Delichatsios  et  al.,  1991): 
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4  (A-CoT^'‘ 


where  C  is  the  coefficient  that  is  less  than  1 .  The  flame  spread  speed  can  still  be  calculated 
using  eqn.  (3-12). 


In  the  foflowing  comparisons,  it  has  been  assumed  that  the  fuel  surface  receives  two 
heat  fluxes:  1)  a  constant  external  heat  flux  over  the  entire  length  of  the  fuel  and,  2)  a  flame 
heat  flux  with  an  exponential  profile,  i.e. 


.  //  .  // 
<lconv=%  ^ 


(3-14) 


In  the  analytical  solution,  the  magnitude  of  this  heat  flux,  Qq,  is  used  in  eqns.  (3-1 1) 
or  (3-13)  to  calculate  the  time  to  pyrolysis  while  the  length  scale,  6,  is  used  in  eqn.  (3-12)  to 
quantify  flame  spread  speed.  It  also  has  to  be  pointed  out  that  the  instantaneous  surface 
temperature,  T,,  is  obtained  by  using  the  constant  external  heat  flux  together  with  the  one 
dimensional  heat-up  model. 

Figure  3-7  shows  the  comparison  of  flame  fi'ont  location  between  Quintiere's  postulate 
and  the  numerical  results  for  constant  heat  flux: 

q'^=5  kW/m^ 

and  the  flame  heat  flux: 

q'eomr^  70  exp[-150(x-Xp)]  kW/m^ 

Three  sets  of  numerical  results  are  illustrated  in  Figure  3-7  with  different  time  steps  and  nodal 
spacings.  Surface  reradiation  is  accounted  for  in  this  comparison.  It  can  be  seen  fi-om  the 
figure  that  all  numerical  results  foUow  the  analytical  solution  quite  well  and  the  numerical 
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result  with  Ax=2  mm  and  At=0. 1  seconds  (the  smallest  time  step  among  the  three)  is  fairly 
close  to  the  analytical  solution.  This  shows  the  numerical  model  predicts  flame  spread 
accurately,  at  least  for  this  case.  Therefore  in  the  following  comparisons,  only  numerical 
results  for  Ax=2  mm  and  At=0.1  .seconds  will  be  .shown. 


TIME  (sec) 


Figure  3-7  Comparison  between  numerical  solutions  of  the  horizontal  flame  spread  model 
and  an  analytical  solution  with  surface  reradiation  and  qe=5  kW/m^ 


The  comparison  of  the  case  with  no  external  heat  flux  (qe=0)  is  illustrated  in  Figure 
3-8.  Two  sets  of  analytical  solutions  using  different  C  (see  eqn.  (3-13))  are  also  plotted.  The 
numerical  results  are  closer  to  that  of  the  analytical  solution  with  C=0.5.  The  close  behavior 
between  the  numerical  result  and  the  analytical  solution  shows  the  validity  of  the  numerical 
model.  It  should  be  pointed  out  that  aU  the  prediction  with  surface  reradiation  loss  use  C=0.5 
including  Figure  3-7.  The  comparisons  at  other  external  heat  fluxes  (qe=10  and  20  kW/m^) 
are  shown  in  Figures  3-9  and  3-10  that  also  manifest  the  accuracy  of  the  numerical  model 
except  more  difference  is  observed  for  external  flux=20  kW/m^.  In  that  case,  the  flame  spread 
speed  is  very  higher  and  from  the  sensitivity  analysis  presented  earlier,  the  time  step.  At, 
should  be  reduced.  An  in^ortant  comment  concerning  Figures  3-7,  3-9  and  3-10  is  that  if  the 
external  heat  flux  is  less  than  surface  reradiation  loss  at  pyrolysis  temperature,  i.e. 

4;<q;>oT; 
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the  flame  spread  will  be  steady  after  initial  transition.  The  numerical  solution  should  match 
with  the  anal)dical  solution  asymptotically  as  shown  in  Figure  3-7.  Otherwise,  the  flame 
spread  speed  will  increase  continuously,  eventually  towards  infinity  (surface  flashes)  if  the 
flame  front  does  not  reach  the  end  of  sample  as  illustrated  in  Figures  3-9  and  3-10.  The 
surface  reradiation  loss  in  this  case  is  about  10  kW/m^. 

In  Figures  3-1 1  through  3-14,  the  same  sets  of  results  for  dq  surface  reradiation  losses 
are  shown.  The  plots  show  good  agreement  between  analytical  solutions  and  numerical 
results.  The  results  from  numerical  solution  at  4^=0  is  essentially  the  same  as  analytical 
solution  because  in  this  case  the  analytical  solution  is  exact.  Discrepancies  between  numerical 
results  and  the  analytical  solutions  have  been  observed  for  external  heat  fluxes  of  10  and  20 
kW/m^  and  in  these  cases,  the  time  steps  may  have  to  be  reduced. 
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Figure  3>8  Comparison  between  numerical  solutions  of  the  horizontal  flame  spread  model 
and  an  analytical  solution  with  surface  reradiation  and  <5^=0  kW/m^.  C  is  a  coefficient 
accounting  for  reradiation  loss 


Figure  3-9  Comparison  between  numerical  solutions  of  the  horizontal  flame  spread  model 
and  an  analytical  solution  with  surface  reradiation  and  qe=10  kW/m^ 
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Figure  3-10  Comparison  between  numerical  solutions  of  the  horizontal  flame  spread 
model  and  an  analytical  solution  with  surface  reradiation  and  q'^=20  kW/m^ 


TIME  (sec) 


Figure  3-11  Comparison  between  numerical  solutions  of  the  horizontal  flame  spread 
model  and  an  analytical  solution  with  no  surface  reradiation  and  q'j=0  kW/m^ 
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Figure  3-12  Comparison  between  numerical  solutions  of  the  horizontal  flame  spread 
model  and  an  analytical  solution  with  eg  surface  reradiation  and  qe=10  kW/m^ 


TIME  (sec) 


Figure  3-13  Comparison  between  numerical  solutions  of  the  horizontal  flame  spread 
model  and  an  analytical  solution  with  no  surface  reradiation  and  qe=5  kW/m^ 
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Figure  3-14  Comparison  between  numerical  solutions  of  the  horizontal  flame  spread 
model  and  an  analytical  solution  with  no  surface  reradiation  and  q"=20  kW/m^ 
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Chapter  4  Description  of  Combustion  Model  and 
Comparison  with  Experimental  Data 


4.1  INTRODUCTION 

The  primary  motivation  for  the  development  of  the  present  integral  combustion  model 
presented  here  has  been  the  necessity  to  predict  flame  spread  behavior  from  burning  of  an 
arbitrary  fuel,  whose  detailed  chemical  composition  would  not  be  known,  as  a  liquid  pool  or 
a  solid  surface.  This  model  would  then  be  coupled  with  other  phenomena  such  as  flame 
spread  or  used  as  a  component  of  enclosure  fire  models. 

Numerical  simulation  of  combustion  processes  would  require  a  direct  numerical 
solution  of  the  full  turbulent  Navier-Stokes  equation  as  well  as  species  relations  and 
incorporation  of  chemical  kinetics.  Simulation  of  a  turbulent  buoyant  diffusion  flame  for  an 
arbitrary  burning  material,  often  in  enclosures,  poses  added  conplexities  where  fuel  geometry, 
air  entraimnent  into  the  flame  and  external  parameters  have  to  be  modeled.  In  such  flows, 
large  turbulent  eddies,  normally  induced  by  buoyancy  and  flame  pulsations,  make  k-e-g  type 
models  prone  to  errors.  Moreover,  application  of  simple  correlations  (which  are  currently 
used  in  compartment  fire  models)  can  not  adequately  capture  the  important  effects  of 
geometry,  flame  radiation  distribution  and  vitiation.  The  use  of  an  integral  approach  for 
modeling  such  phenomena  has  shown  to  be  promising  (Delichatsios  and  Mathews,  1989; 
Tamanini,  1981),  is  accurate  within  experimental  validation  errors  and  computationally  much 
less  intensive  than  full  field  models.  The  integral  models  can  include  physics  of  complex 
combustion  flows,  such  as  those  encountered  for  pool  fires  (i.e.  large-scale  turbulent 
fluctuations),  via  use  of  experimentally  obtained  parameters,  properties  and  correlations 
whereas  k-e-g  models,  which  have  been  calibrated  for  isothermal  flows  having  relatively  small 
fluctuations,  would  simply  be  unreliable.  Prediction  of  flame  temperature  for  a  pool  fire  by 
a  k-e-g  model  may  be  in  error  by  as  much  as  200°C,  as  reported  by  Nicolette  et  al.  (1993). 

Use  of  an  integral  combustion  model  for  a  turbulent  diffusion  flame,  described  in  this 
chapter,  has  several  advantages  and  characteristics.  This  approach  offers  a  robust  method  in 
order  to  predict  soot  formation,  flame  radiation  for  optically  thin  flames  (and  applicable  to 
optically  thick  cases),  includes  turbulent  fluctuation  effects,  major  species  production,  degree 
of  completion  of  combustion  for  arbitrary  flammable  material  by  using  flammabihty 
properties.  For  soot  radiation,  the  laminar  smoke-point  heat  release  rate,  can  be  utilized 
to  predict  soot  concentrations  (Delichatsios  et  al.,  1992).  The  model  could  be  applicable  to 
wsdl,  pool  and  jet  fires  and  incorporate  vitiation  effects.  The  use  of  experimentally  determined 
parameters,  properties  and  correlations  (e.g.  entrainment  rates)  allows  the  model  to  include 
the  above  listed  physics  and  have  such  a  wide  range  of  applicabilities.  The  major  species  yield 
(i.e.  CO,  CO2,  etc)  are  determined  by  using  state  relationship  for  conserved  scalar  and  the 
flamelet  concept.  In  addition,  the  turbulent  fluctuations  of  the  conserved  scalar  are  modeled 
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using  a  direct  empirical  relationship,  which  has  been  tested  carefully  for  buoyant  jet  difiusion 
flames  with  very  promising  results. 

4.2  MODEL  DESCRIPTION  AND  ENERGY  BALANCE 


The  mtegral  combustion  model  presented  earlier  (Chen  et  al.,  1993)  has  been 
developed  for  turbulent  jet  difiusion  flames  and  includes,  as  its  most  important  component, 
flame  radiation  and  its  relation  to  flame  turbulence  and  the  sooting  tendency  of  the  fliel.  The 
buoyant  jet  is  modelled  using  integral  conservation  equations  and  an  ensemble  offlamelets 
to  model  reactions  and  combustion.  The 
energy  equation  is  simultaneously  solved 


with  the  flow  equations,  accounting  for 
radiation  losses  and  employs  a  simple 
soot  formation  model.  The  model 
computes  the  radiant  losses,  the 
incompleteness  of  combustion  and  the 
CO-yield  by  using  the  following  physics: 
1)  the  lower  part  of  the  flame  is  modelled 
by  difiusion  flamelets  using  a  probability 
distribution  function  and  2)  the  upper 
part,  as  the  flame  cools  due  to  radiant 
losses,  is  modelled  as  a  premixed  mixture 
sustaining  CO+soot  oxidation  (Figure  4- 
1).  The  main  characteristics  of  the 
integral  model  are  the  following: 

1)  Instead  of  an  entrainment  equation,  a 
dimensionless  correlation  for  entrainment 
spanning  the  range  fi’om  momentum  to 
buoyant  turbulent  jet  flames  has  been 
employed  and  validated  (Chen  et  al., 
1993).  2)  A  direct  relationship  for  the 
fluctuation  of  a  conserved  scalar, 


(kW/m) 

Figure  4-1  A  schematic  diagram  of  the  integral 
combustion  model 


inspired  by  recent  work  (Chatwin  and 

Sullivan,  1990)  and  validated  by  recent  experiments  is  also  used  (see  Appendbc  B  for  further 


detail): 


_  1 

af(Pl,-f)'  a(PV!) 

1+Y - ^ 

l-K 


(4-1) 
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3)  In  the  present  model,  radial  profiles  for  the  mixture  fi'action  and  velocity  are  used  in  order 
to  preserve  the  strong  nonlinearity  of  the  combustion  process  (Delichatsios  and 
Mathews,  1989),  i.e.: 


(4-la) 

(4-lb) 


where  u  and  Uc  are  the  mean  local  and  centerline  velocities  while  T  and  are  the  stream 
function  and  entrainment  rate,  respectively  (Delichatsios  and  Mathews,  1989). 

4)  Most  importantly,  a  flame  radiation  model  applicable  for  luminous  (soot  radiation)  and 
optically  thin  flame  is  presented  whose  sooting  tendency  is  characterized  by  the  laminar 
smoke-point  height. 


The  soot  radiation  model  is  founded  on  a  phenomenological  description  for  soot 
formation  which  has  been  established  fi-om  data  in  laminar  difiusion  flames  (Delichatsios, 
1993);  the  soot  concentration  can  be  determined  by  the  following  rate  equation 
(Delichatsisos,  1993): 


1-^ 


-E/RT 


(4-2a) 


wherein  ^  are  the  fuel  mixture  fi-action  and  its  stoichiometric  value,  E  is  the  activation 
energy  for  soot  formation  and  T  is  the  absolute  temperature.  As  an  approximation  the 
destruction  of  soot  can  be  ignored.  The  proportionality  factor  Aj  which  includes  the  fuel 
structure  chemistry  effects,  is  proportional  to  the  inverse  of  the  laminar  smoke-point  heat 
release  rate  (Delichatsisos,  1993): 


(4-2b) 


The  proportionality  constant  is  properly  defined  in  literature  (Delichatsisos,  1993). 

Prior  to  applying  eqn.  (4-2a)  to  turbulent  flows,  it  is  useful  to  know  what  is  required 
for  modelling  radiation  losses  in  an  optically  thin  medium.  The  flame  radiation  loss  due  to 
soot  can  be  treated  in  this  case  as  an  volumetric  emission  described  by  an  expression  which 
can  be  reduced  to  the  following  form  (Delichatsisos,  1993): 

^f=3.26-10-^— 7/^  W/m^K^  (4.3) 


where  P  is  the  pressure,  R  is  the  gas  constant  and  is  the  soot  density. 
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For  a  turbulent  flow,  one  has  to  average  the  volumetric  radiant  loss  (eqn.  (4-3))  over 
the  turbulent  fluctuating  field.  Local  instantaneous  values  of  the  temperature,  T,  and  the  soot 
concentration,  Yj,  needed  for  this  purpose,  are  obtained  by  using  the  representation  of 
turbulent  combustion  by  an  ensemble  of  flamelets  that  may  also  be  subjected  to  a  spectmm 
of  turbulent  straining  rates. 

Following  such  modelling,  the  instantaneous  laminar  flamelet  temperature,  Tl,  would 
be  a  function  of  the  mixture  fraction,  as  it  is  customarily  done  in  flamelet  models  (a 
modification  to  this  relationship  will  be  presented  later  in  this  work  to  account  for  the  radiant 
losses).  The  soot  concentration  can  be  obtained  by  integration  of  the  soot  formation  rate 
(eqn.  (4-2a))  over  the  local  flow  field  of  a  strained  diffusion  flame  (Delichatsisos,  1993).  The 
soot  concentration  becomes: 

Qsp 

wMe  Y^-0, 

where  y  is  the  inverse  straining  rate  and  F  is  a  function  obtained  by  integrating  eqn.  (4-2a) 
over  an  appropriate  streamline  (from  0  to  Yj).  Its  exact  form  is  not  required  for  the  present 
work,  which  focuses  on  optically  thin  flames.  Note  that  eqn.  (4-4b)  has  been  derived  from 
a  sinplification  of  the  oxidation  process  according  to  which  the  oxidation  rate  in  the  fuel  lean 
side  (^<Q  is  so  fast  that  soot  oxidizes  instantaneously  there.  This  situation  is  applied  as  long 
as  the  (maximum)  laminar  flamelet  temperature  (at  stoichiometric  condition)  is  greater  than 
a  certain  value,  about  (Delichatsisos  et  al.,  1992)  1500  K;  otherwise,  soot  can  not  be  further 
oxidized  and  is  emitted  as  smoke. 

By  using  eqn.  (4-4)  and  the  function  of  tenperature  in  terms  of  mixture  fraction,  T(0» 
one  obtains  from  eqn.  (4-3)  an  average  volumetric  radiation  loss: 

jf=3.26-10^— 

=3.26-lO-"^|/-X.F(5.5^r^)r\5)P(5,Y)^5rfY  (4-5a) 

wherein  P(?,y)  is  the  joint  pdf  between  the  mixture  fraction,  and  the  straining  rate,  y. 
Next,  we  use  the  common  approximation  that  mixture  fraction  and  straining  rate  are 
statistically  independent  so  that  eqn.  (4-5a)  becomes: 

_  1 

4"'^3.26-10-^-^^fFT*Pa)dl  (4.5b) 


(4-4a) 

(4-4b) 
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It  is  this  relationship  that  is  applied  in  the  energy  balance  of  the  present  integral 
combustion  model  along  the  jet  axis: 

■^jhdiff^-jq'J'lurdr  (4-6a) 

0 


where  h  is  the  total  enthalpy  of  the  product  of  combustion  (i.e.  chemical+convective)  and  T 
is  the  stream  function  for  the  mean  axial  jet  velocity.  In  this  form,  eqn.  (4-6a)  is  difficult  to 
solve  owing  to  the  iterations  needed  to  obtain  the  flame  temperature  for  the  laminar  flamelets 
or  for  the  turbulent  flow. 


In  the  spirit  of  the  integral  model  (and  to  save  computational  time),  we  make  the 
approximation  that  the  laminar  flamelet  temperature  at  a  given  height  from  the  source  is 
uniformly  reduced  by  a  constant  fraction  a  of  its  adiabatic  value,  namely: 

TS)-K‘a(Tj^iyTj  (4-6b) 

while  the  dependence  of  species  concentrations  on  ^  remains  the  same  independent  of  the 
radiation  losses  as  long  as  the  maximum  flamelet  temperature  is  greater  than  1,500  K(this 
approximation  is  supported  by  experiments  (Koylu  et  al.,  1992)).  With  these  approximations, 
the  energy  conservation  equation  (4-6a)  becomes  (see  Appendix  A): 


(4-6c) 


where  2^(z)  is  the  chemical  (heat  of  formation  flux)  heat  release  at  height  z  and  (1-a)  can 
also  be  identified  to  be  equal  to  the  radiant  fraction  up  to  height  z: 


Qd^) 


(4-6d) 


where  Qr(z)  is  the  radiation  losses  from  the  flame  over  a  length  extending  from  the  source 
to  a  given  height  z.  Even  in  the  sin5)lified  form,  eqn.  (4-6c)  requires  iterations  over  the  whole 
field  in  order  to  obtain  the  value  of  a  owing  to  the  complicating  dependence  of  radiant  losses 
on  temperature  as  indicated  in  eqn.  (4-5b).  A  further  simplification  of  eqn.  (4-6c)  has  been 
obtained  based  on  overall  radiation  correlation  (Delichatsios  et  al.,  1992)  and  the  specific 
form  of  radiations  losses  (eqn.  (4-5b)): 


jq'"lnrdr^ 

0 


J  P" 


(4-7a) 


The  following  remarks  are  made  concerning  the  derivation  of  eqn.  (4-7a): 
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1 .  Most  of  the  contributions  to  radiation  losses  originates  for  values  of  ^  ~  near  the 
stoichiometric  value  so  that: 

FT  delta  function  ~6(5-5,)  (4-7b) 


2.  The  average  strain  rate  y  (which  can  be  identified  with  the  Kolmogorov  strain  rate) 
is  nearly  independent  of  location  or  total  heat  release  rate  for  buoyant  flames 
(Delichatsisos  et  al.,  1992);  therefore  it  is  assumed  to  have  a  constant  value. 

3.  The  reduction  of  the  term  FT*  due  to  radiation  losses  is  proportional  to  1/Xr^-  This 
dependence  has  been  deduced  from  global  correlations  on  radiation  (Delichatsisos  et 
al.,  1992). 

4.  The  constant  C  incorporates  all  the  physical  constants  associated  with  the 
assumptions  1, 2, 3.  It  has  to  be  determined  from  experiments. 

5.  Comparison  with  experiments  demonstrates  the  applicability  of  the  approximations 
used  to  derive  eqn.  (4-6a).  More  importaiitly,  the  value  of  the  constant  C  should  be 
independent  of  the  fuel,  assuming  that  the  present  approximations  are  adequate. 


Finally,  with  these  approximations,  the  energy  balance  equation  (4-6c)  becomes: 


4  If, 


(4-8) 


We  niay  also  characterize  the  radiation  modelling  in  eqn.  (4-8)  in  the  following  way: 

1) .  All  radiation  originates  from  soot  that  is  near  the  instantaneous  maximum  flamelet 

temperature; 

2) .  This  flamelet  temperature,  however,  decreases  due  to  radiant  losses  as  the  radiation 

fraction  increases. 


Equations  (4-8)  is  singular  near  the  flame  base  where  (l-a)-O.  We  can  not  model  the 
flow  conditions  near  the  flame  base  at  this  moment.  Fortunately,  because  we  have  used 
integrated  value  over  the  entire  length  of  the  flame,  this  effect  has  been  eliminated. 


4.3  THE  STRUCTURE  OF  THE  COMPUTER  CODE 

A  computer  code  has  been  developed  based  upon  the  above  approach.  The  flow 
chart  of  this  code  is  shown  in  Figure  4-2. 
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PROGRAM  RADH.F 


MAJN 

DATA  INPUT  (FILE’  comhb.dat') 

CALHTO  Enthalpies  of  Amb.  and  Source 

- ZBRENTl  Find  Laminar  Temperature 

^ - CALTE  Enthalpy  Calculation 

-  QGAUS  Calling  Comb.-Radiation  Calculat 
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von 

-  ENTHA 

Species  Concentration  Cal. 

' - SPECl  or  SPEC 2 

Define  Species  Cone. 

I -  CALCON 

State  Relationship 

-  QGAUS3 

-  QRR 

pdf  Calculation 

RHS  of  the  Integrals 
RESULT  OUTPUT 


STOP 


Figure  4-2  Subroutines  and  interconnections  of  the  integral  combustion  model 
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4.4  COMPARISON  OF 
PREDICTION  WITH 
EXPERIMENTS 

Eqn.  (4-8)  together  with  the 
remaining  conservation  equations 
(i.e.,  momentum,  mixture  fraction), 
as  well  as  with  the  flamelet 
combustion  model,  form  a  well- 
behaved  system  of  equations  which 
have  been  applied  to  several  cases. 
Results  and  comparison  for  the 
radiation  distribution  are  presented 
in  this  section  for  turbulent  buoyant 
jets  and  pool  fires  of  various  fuels. 


.  Figure  4-3  Chemical  heat  release  rate  normalized 

re  iction  o  eat  release  by  the  total  heat  release  rate  versus  distance  from 
rate  from  the  integral  combustion  the  nozzle  exit,  normalized  by  twice  the  length,!,,, 
model  and  companson  with  over  which  half  of  the  combustion  has  been 
expenmental  data  (Tamanim,  1982)  completed 


Figure  4-4a  Radiation  fraction  at  a  given  height 
plotted  versus  the  height  normalized  by  twice  the 
length,  Lo  5,  over  which  half  of  the  combustion  has 
been  completed 


for  a  propane  (CjHg)  buoyant  jet 
flame  are  illustrated  in  Figure  4-3. 
The  abscissa  is  the  ratio  between 
the  height  from  the  nozzle  to  two 
times  Lq  j  is  the  distance  at 
which  one-half  of  total  heat  release 
occurs)  and  the  ordinate  is  the 
scaled  heat  release  rate,  i.e.  the 
ratio  of  the  heat  release  rate  at  that 
height,  2ft(z)  to  the  total  heat  rate, 
Qff,.  The  burner  diameter  used  in 
the  experiment  was  12.7  mm  and 
the  jet  flame  theoretical  heat  release 
rate  was  30.8  kW.  The  comparison 
in  Figure  4-3  shows  very  good 
agreement  between  prediction  from 
the  model  and  the  experimental 
data  for  heat  release  rate. 


The  constant  C,  in  the 

energy  balance  eqn.  (4-8),  was  determined  by  requiring  that  at  the  end  of  combustion  the 


42 


radiation  fraction  (1-a)  is  equal  to  the 
total  radiation  fraction  for  the  fuel, 
which  is  known  if  the  smoke-point 
heat  release  rate,  Qsp,  is  known 
(Delichatsisoset  al.,  1992)  (for  propane, 

<2^=303.4  W).  The  value  of  this 
constant  was  found  to  be,  C=0.36  W/s 
if  is  expressed  in  W.  By  using  this 
value,  the  calculated  radiation  fraction, 

Xr  vs.  the  scaled  height  from  the 
nozzle  is  plotted  in  Figure  4-4a  while 
the  scaled  radiation  per  unit  length  of 
the  flame  is  illustrated  in  Figures  4-4b 
and  4-4c.  AH  the  experimental  data 
were  obtained  by  Tamanini  (1981). 

The  close  behavior  between  the 
prediction  and  the  experiment,  shown 
in  Rgure  4-4a  reinforces  the  validity  of 
this  approach.  The  prediction  does 
give  the  right  trend  of  radiation  loss 
distribution  (Figures  4-4b  and  4-4c) 

though  the  difference  with  experiments  is  fairly  large.  It  can  be  noted  (may  not  be  obvious) 
that  the  area  under  the  curves  from  the  prediction  (Figure  4-4b)  is  the  same  as  that  from  the 
experiment. 


Figure  4-4b  Radiation  distribution  per  unit  height 
versus  normalized  distance  from  the  nozzle  exit  (cf. 
Fig.  4-4a).  The  experimental  data  are  from  Tamanini 
(1982) 


Finally,  Figure  4-5  illustrates 
the  maximum  laminar  flamelet 
temperature  variation  as  a  function  of 
the  height  from  the  flame  base.  The 
result  shows  decreasing  laminar  flame 
temperature  as  radiation  losses 
increase.  Near  the  flame  tip  (~  1  m), 
the  laminar  flamelet  temperature  is 
about  1680  K  indicating  that  no 
extinguishment  of  the  flame  due  to 
radiation  loss  occurs.  Therefore,  soot 
is  probably  completely  oxidized  inside 
the  flames,  which  agrees  with 
experiments  (Tamanini,  1981)  that  little 
smoke  is  emitted  from  propane  flames. 


Figure  4-4c  Same  data  as  in  Figure  4-4b,  except 
that  the  radiation  per  unit  height  is  normalized  by 
its  maximum  value.  The  experimental  data  are  from 
Tamanini  (1982) 


43 


The  same  constant  C  (=0.36  W/s).  can  be  used  to  calculate  the  combustion  and 
radiation  of  another  fuel,  propylene,  if  the  smoke-point  heat  release  rate,  is  known  (see 

eqn.  (4-8)).  We  have  used  this 


0.00  0.20  0.40  0.60  0.80  1.00 


HEIGHT  FROM  THE  FLAME  BASE  (m) 


approach  to  predict  the  radiation  of 
propylene  (C^Hg,  6^=463  W).  The 
results  are  shown  in  Figures4-6(a- 
c).  In  Figure  4-6a,  the  prediction 
of  radiation  fraction  is  plotted 
against  the  height  from  the  flame; 
the  total  radiation  fraction  is  about 
0.46  which  is  close  to  0.44  as 
measured  in  an  experiment 
(Delichatsisos  et  al.,  1992).  Figure 
4-6b  shows  a  comparison  of  non- 
dimensional  radiation  distribution 
between  the  prediction  and 
experiments  (Delichatsisos  et  al., 
1988).  These  results  show  that  one 
can  use  the  smoke-point  heat 


Figure  4-5  Reduction  of  laminar  flamelet  release  rate  to  predict  total 

temperature  with  height  owing  to  radiation  losses  radiation  and  radiation  distribution 
for  the  use  of  Figure  4-3  from  buoyant  jet  flames.  For 

conpleteness,  the  maximum  laminar 


flamelet  temperature  for  propylene 
flames  is  shown  in  Figure  4-6c.  In 
contrast  to  the  propane  flames. 
Figure  4-5,  the  temperature  drops 
below  1,500  K  at  a  height  at  0.60 
m;  such  a  low  temperature  will 
result  in  flame  quenching,  copious 
smoke  yield  and  incompleteness  of 
combustion,  which  might  explain 
the  higher  calculated  value  of  total 
radiation  fraction,  0.46,  relative  to 
the  numerical  value  of  0.44. 
Incompleteness  of  combustion  will 
be  addressed  as  the  present  model  is 
carefully  evolved. 

The  above  discussions  and 
results  presented  in  Figures  4-6a 
and  4-6b  have  shown  that  use  of  a 
single  constant  (representing 
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Figure  4-6a  Prediction  of  radiation  fraction  versus 
height  for  propylene  turbulent  buoyant  jet  flame.  No 
experiments  are  available  except  for  the  total 
radiation  fraction 
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collective  physical  parameters)  is  feasible  for  radiation  modeling.  This  also  points  out  that 
the  integral  model  has  accurately  included  the  essential  physics  for  buoyant  jet  flames. 

Another  ability  of  the  combustion  mode  is  to  predict  the  products  of  combustion. 
Figure  4-7  shows  the  comparison  of  prediction  with  experiment  of  CO  production  rate  for 
propane  jet  flame. 

A  final  demonstration  that  the  present  model  could  be  extended  to  combustion  and 
radiation  in  pool  fires  is  shown  in  Figures  4-8a,  and  4-8b,  for  a  propane  fire  0.38  m  in 
diameter.  The  same  conservation  equations  were  used  for  a  turbulent  buoyant  jet  flame 
including  the  constant,  C  (=0.36  W/s)  for  radiation  of  C3Hg.  The  following  modifications 
were  made: 


1)  an  entrainment  rate  correlation  for  pool  fires  was  used  (Delichatsisos,  1987)  and 

2)  the  fluctuations  were  represented  by  the  relationship: 


5(1-5) 


0.11 


(4-9) 


instead  of  eqn  (4-1)  because  eqn. 

(4-9)  was  shown  (Delichatsisos  and 
Mathews,  1989)  to  be  the  one 
applicable  for  pool  fires.  This 
expression  represents  a  level  of 
fluctuations,  (1’'^)^/T,  which  for 
values  of  ^  close  to  1  is  smaller  than 
what  eqn.  (la)  represents;  this 
behavior  is  consistent  with 
measurements  of  fluctuations  in 
pool  fires  (Gengebre,  1984). 

The  agreement  between 
predictions  and  experiments  as 
shown  in  Figures  4-8a  and  5-8b  is 
more  encouraging  than  should  be 
expected,  because  the  elliptic  flow 
field  near  the  pool  base  is  not 
modeled  in  the  present  approach.  In 
fairness  to  the  present  work,  it 

should  be  pointed  out  that  measurements  of  Xr  in  pool  fires  in  Tamanini's  (1982)  experiments 
(wherein  he  used  a  "catcher-quencher"  to  measure  heat  release  rates)  have  large  errors 
(±50%)  near  the  pool  base  (<.20  m)  because  a)  the  width  of  the  slit  radiometer  used  to 
measure  the  radiation  per  unit  height  was  54  mm  and  b)  the  "catcher-quencher"  might  have 


Figure  4-6b  Radiation  distribution  per  unit  height 
versus  height  for  a  propylene  turbulent  buoyant  jet 
flame 
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LAMINAR  FLAMELET  TEMEPRATURE  (K) 


influenced  the  combustion.  A  more  detailed  examination  of  the  shortcoming  of  pool  fire 
modeling  is  currently  being  undertaken. 


HEIGHT  FROM  THE  FLAME  BASE  (m) 


Figure  4-6c  Reduction  of  laminar  flamelet  temperature  vs.  height  owing  to 
radiation  loss  for  the  case  of  Figure  4-6a 
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Figure  4-7  CO  production  rate  vs.  scaled  height  from  the  flame  for  a  propane 
jet  flame 


HEIGHT  FROM  THE  FLAME  BASE  (m) 


Figure  4-8a  Chemical  heat  release  rate  versus  distance  from  the  base  of  a  propane  pool 
fire  0.38  in  diameter 
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THE  NOZZLE  DIAMETER  IS  38  cm 
THE  EXPERIMENTAL  DATA  ARE  FROM 
Tamanini  (1982) 


Tamanini's  POOL  FIRE  DATA 
-  PREDIOTION 


I  I  I  I  I  I  I  I  I  "I  I  I  I""!  I  T-  I  I  I  I  I  I  I  I  I  I  I  1  i'  I  I  "i'l-'T  I  I  I  I  I  1  I  I  I  ■]  T'  I  1  I  I  I  i  I  i' 1  I  1 

.OO  0.20  0.4-0  0.60  0.80  1.00  1.20 

HEIGHT  FROM  THE  FLAME  BASE  (m) 


Figure  4-8b  Radiation  friction  at  a  given  height  plotted  versus  the  height  from  the  base  of 
a  propane  pool  fire  0.38  m  in  diameter 


Chapter  5 


Description  of  Constant  Speed  Horizontal 
Flame  Spread  (CSHFS)  Experiment 


Following  the  completion  of  the  previous  analysis,  the  horizontal  flame  spread  code 
is  complemented  with  measurements  of  gas  phase  conduction  heat  flux  by  conducting  a 
constant  speed  horizontal  flame  spread  experiment  (CSHFS).  This  includes  testing  with  non¬ 
charring  (black  PMMA,  both  19  mm  (3/4")  and  6.35  mm  (1/4")  thick,  manufactured  by  the 
Polycast  Technology  Corporation)  and  charring  (particle  board,  19  mm  thick,  manufactured 
by  the  Potlatch,  Inc.)  materials. 

As  discussed  previously,  the  quantification  of  the  total  flame  heat  flux  immediately 
ahead  of  the  flame  front  is  very  important.  There  is  a  gas-phase  conductive  flux  which 
applies  over  a  very  short  distance  (in  the  order  of  1  mm)  just  ahead  of  flame  front.  To  try  to 
quantify  this  flux  and  verify  the  flame  spread  model,  a  particular  scheme  is  accomplished.  In 
this  approach,  the  fuel  sample  is  placed  on  a  movable  platform.  The  platform  speed  can  be 
controlled  using  a  stepping  motor. 

The  burning  sample  was  moved  by  a  slide  (manufactured  by  the  Unislide,  Inc.)  at  a 
constant  speed  in  the  experiments  to  obtain  a  constant  flame  burning  area  (Figure  5-1).  The 
transient  effects  can  therefore  be  eliminated.  The  flame  radiation  can  also  be  easily  deduced 
from  the  experiments  (flux  gage  measurements)  or  predicted  from  the  combustion-radiation 
model  as  described  in  Chapter  4.  The  only  remaining  unknown  is  the  gas-phase  conduction 
ahead  of  the  flame  front.  This  conductive  heat  flux  has  to  be  matched  with  the  flame  spread 
speed  (also  the  sample  moving  speed)  and  the  temperature  measurements  from  the 
thermocouples  as  described  later. 

To  create  a  condition  whereby  the  flame  spread  speed,  Vf,  is  equal  to  the  sample 
speed,  Vj,  the  flame  is  being  quenched  by  a  water-cooled  plate  placed  above  the  fuel  sample. 
The  quenching  serves  two  purposes.  First,  the  flame  impingement  on  the  radiant  panel  is 
limited  or  eliminated.  This  impingement  creates  an  additional  flux  which  can  not  be 
accounted  for.  Second,  the  quenching  helps  to  keep  the  flame  spread  speed  and  the  burning 
area  constant.  The  heat  release  rate  would  also  remain  constant  for  a  fixed  burning  area  and 
a  constant  external  flux  distribution  by  satisfying  the  condition  of  =  Vj.  The  external  heat 
flux  at  point  A  (see  Fig.  5-1)  wfll  always  be  the  same  under  these  conditions  since  the  external 
heat  flux  supplied  by  the  radiant  panel  is  decaying  along  the  sample  surface  as  seen  in  Figure 
5-2.  We  have  used  the  LIFT  apparatus  to  accomplish  this  The  challenge  in  this  scheme 
would  be  to  have  a  very  short  transient  time  after  ignition  and  establishing  the  steady  flame, 
positioning  of  the  pilot  flame,  instrumentation  and  measurements  and  finally  the  control  of 
the  movement  of  the  fuel  sample.  Thus,  it  is  possible  to  better  quantify  the  gas-phase 
conductive  heat  flux  from  the  flame  when  the  external  heat  flux  and  flame  radiation  are 
known. 
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Radiant  Panel 


Given  the  constraint  of  the  LIFT  apparatus,  a  slide  with  a  moving  distance  of  600  mm 
was  used  to  move  the  fuel  samples.  The  fuel  samples  are  instrumented  with  thermocouples 
and  heat  flux  gages.  Three  temperature  measuring  stations  situated  roughly  200  mm  (The 
exact  distances  between  two  stations  varies  depending  on  a  particular  experiment,  the 
variations  could  be  as  large  as  ±50  mm  and  the  distances  have  been  listed  in  Appendix  H) 
apart  are  used.  The  first  station  is  about  200  mm  from  the  leading  edge  of  the  sample.  At 
the  second  station,  a  1.6  mm  (1/16")  diameter  small  heat  flux  gage  is  also  used.  At  the  third 
station,  which  is  about  600  mm  from  the  leading  edge  of  the  sample,  a  25.4  mm  (1 ")  diameter 
(probe  size  9. 1  mm  (23/64")+copper  plate)  large  heat  flux  gage  is  instrumented.  The  flux 
gages  are  supported  in  the  bottom  of  the  marinite  board  insulation  to  insure  their  sensors 
flush  with  the  sanple  surface.  The  small  flux  gage  is  employed  since  the  flame  front  can  not 
go  over  the  large  sensor  (actually,  the  flame  will  go  around  the  sensor)  due  to  the  cool  sensor 
surface  (cooled  by  water).  K-type  thermocouples  (0.127  nmi  (0.005")  in  diameter)  are  used 
to  measure  surface  and  gas-phase  temperatures.  The  surface  thermocouples  are  imbedded 
by  running  two  thermocouple  wires  through  two  holes  and  righting  the  two  wires  at  the  back 
of  the  fuel  sample  to  maintain  good  contacts  between  thermocouple  beads  and  sample 
surface.  The  gas  phase  thermocouples  are  placed  just  above  the  sample  surface  (above  1 
mm).  It  is  expected  that  we  can  use  the  differences  between  the  surface  and  gas  temperatures 
to  quantify  the  gas-phase  conduction.  In  each  test,  video  recording  was  also  made  to  obtain 
the  flame  front  histories.  A  constant  burning  area  can  be  observed  from  the  video  tape. 

Figures  5-3  through  5-10  are  pictures  of  the  experimental  set-up  and  the  experiments. 
Figure  5-3  shows  the  long  view  of  the  LIFT/CSHFS  apparatus.  Figure  5-4  demonstrates 
early  stage  of  the  experiment.  Figures  5-5  and  5-6  illustrate  close-up  pictures  of  the 
experiment.  One  can  see  that  a  metal  bracket  with  scale  longitudinally  holds  a  piece  of 
burning  material  (in  this  case,  it  is  a  19  mm  particle  board).  On  the  left  of  the  flame,  there  is 
a  water  cooling  plate  which  is  very  close  to  the  sample  surface  so  that  there  is  no  burning 
underneath  the  cooling  plate  (the  sample  moves  from  right  to  left).  At  the  top  left  comer  of 
the  pictures,  the  radiant  panel  is  shown  (located  at  a  cute  angle  relevant  to  the  fuel  surface). 
Underneath  the  sample  holder  one  can  see  an  attachment  that  has  all  the  thermocouples  and 
heat  flux  connectors  so  that  the  extension  wires  can  be  moved  more  freely  with  the  sample. 
The  flame  width  is  about  7  cm  for  these  tests.  The  middle  stage  of  the  experiment  is  shown 
in  Figure  5-7  while  side  view  of  the  experiment  is  manifested  in  Figure  5-8.  A  large  heat  flux 
gage  probe  (circular  shape  in  the  middle  of  the  sample)  is  clearly  shown  in  Figure  5-8.  The 
flux  gage  probe  is  flush  with  the  san^le  surface  and  it  is  mounted  by  drilling  a  25.4  diameter 
hole  through  the  san5)le,  supported  at  the  back  of  the  sample  insulation  (25.4  thick  marinite 
board).  Figure  5-9  shows  the  experiment  is  just  stopped,  and  the  PMMA  samples  after  the 
experiments  are  shown  in  Figure  5-10. 

In  order  to  determine  the  accuracies  of  the  large  and  small  heat  flux  gages,  we  have 
compared  the  two  heat  flux  gage  readings  at  the  same  location  for  an  inert  test  (no  flame, 
only  radiant  panel  was  on).  The  small  heat  flux  gage  gave  higher  values  than  the  large  heat 
flux  gage  as  shown  in  Appendix  H.  The  large  flux  gage  was  painted  by  using  a  known 
emissivity  paint  and  then  calibrated  using  an  Factory  Mutual  apparatus  and  the  result  was  not 
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much  different  from  the  calibration  newly  done  by  the  manufacturer.  The  large  flux  gage 
reading  was  also  compared  with  a  flux  gage  calibrated  at  NIST  and  the  difference  was 
negligible.  Finally,  during  the  inert  test,  after  exposed  both  heat  gages  to  the  radiant  panel, 
the  flux  gages  were  covered.  The  small  flux  gage  went  to  a  non-zero  reading  (it  varied  from 
time  to  time!).  Therefore  it  has  been  determined  that  small  flux  gage  readings  were  not 
accurate  and  that  the  measurements  of  large  heat  flux  gage  is  more  accurate  (away  from 
flame  front)  (see  also  Figure  6-3  in  next  Chapter).  The  cause  of  non-accurate  readings  from 
small  flux  gage  is  not  fully  understood,  it  may  be  due  to  the  ineffective  cooling  of  the  sensor 
because  of  the  small  sensor/cooling  area.  The  other  problem  is  that  because  the  small  heat 
flux  gage  nose  is  fragile,  a  slightly  larger  hole  was  drilled  to  mount  the  small  flux  gage.  There 
was  some  open  space  between  the  gage  nose  and  the  wall  of  the  hole.  During  the 
experiment,  the  insulation  wall  as  well  as  the  fuel  itself  was  heated,  some  heating  from  the 
wall  to  the  gage  may  exist  and  that  would  cause  some  problems  for  cooling  the  sensor. 


52 


Figure  5-6  Same  as  Figure  5-5.  but  tor  dilTereiiL  tesi 


'igure  5-7  Flame  spread  at  the  middle  stage  of  the  experiment 


Figure  5-9  The  slide  slops,  the  experiment  finishes 


Figure  5- 10  The  PMMA  samples  alter  expcriinenls  (6.35  mm  ( 1/4" )  thick 
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Chapter  6  Experimental  Results  and  Comparison  of 
Simulation  of  Horizontal  Flame  Spread 
Model 

6.1  EXPERIMENTAL  APPROACH  AND  PROCEDURE 

Two  fuels,  PMMA  (non-charring)  and  particle  board  (charring  materials),  for  which 
flammability  properties  have  been  known  (Chen  et  al.,  1995),  were  used  in  the  CSHFS 
experiments  described  in  Chapter  5.  Two  thicknesses  of  PMMA  (19  mm  (3/4")  and  6.35  mm 
(1/4"))  and  a  single  thickness  of  particle  board  (19  mm)  particle  board  samples  were  tested. 
The  dimensions  of  the  san^jles  were  152.4x600  (WxL)  mm.  In  the  tests,  the  samples  were 
preheated  for  a  specified  time  and  this  time  was  determined  before  an  experiment  using  the 
pre-heat  model  of  the  flame  spread  model  discussed  in  Chapter  3.  This  time  depends  on  the 
external  heat  flux  level;  then  an  ignition  flame  was  applied  at  the  leading  edge  of  the  sample 
respectively.  The  sample  (slide)  was  moved  at  a  pre-determined  constant  speed  after  flame 
has  grown  to  certain  size  (established  flame).  This  speed  was  determined  by  using  the  flame 
spread  model  and  assuming  flame  heat  fluxes.  This  speed  was  verified  or  altered  based  on 
a  preliminary  experiment.  The  criterion  for  choosing  the  speed  was  that  the  burning  area 
should  not  be  too  large  (if  the  speed  was  too  small)  such  that  the  flame  would  not  impinge 
on  the  radiant  panel  and  create  an  additional  heat  flux  that  can  not  be  quantified.  The  burning 
area  could  not  be  too  small  either  (if  the  speed  was  too  large)  because  relative  errors 
introduced  by  fluctuations  of  flame  front  could  become  too  large  or  all  the  flame  would  be 
quenched.  The  proper  burning  width  is  10  cm  from  experience.  There  were  three  measuring 
stations  along  the  burning  sample  and  they  were  approximately  200, 400  and  600  mm  away 
from  the  leading  edge  of  the  sample.  The  surface  temperatures  at  three  stations,  gas 
temperatures  at  first  two  stations  and  total  heat  fluxes  at  last  two  stations  were  measured. 
At  the  second  station,  a  small  heat  flux  gage  (1.6  mm  (1/16")  probe  diameter)  was  used  and 
a  9.1  mm  (23/64")  sensor  diameter  heat  flux  gage  was  used  at  the  third  station.  All  heat  flux 
gages  were  Gardon  type  manufactured  by  Medtherm,  Inc. 

Data  obtained  from  the  experiments  described  in  Chapter  5  include:  1)  temperature¬ 
time  histories  from  thermocouple  measurements,  2)  heat  flux-time  histories  from  the  heat  flux 
gage  readings  and,  3)  flame  front  location  relative  to  the  original  leading  edge  of  sample 
which  was  obtained  from  both  visual  observations  and  video  tape  recordings.  The  difference 
between  the  total  heat  flux  measurement  and  the  external  heat  flux  at  the  measuring  stations 
would  be  the  total  flame  heat  flux  incident.  The  distance  from  the  cooling  plate  edge  to  the 
measuring  station  at  any  instant  of  time  is  also  known  because  the  time  to  start  moving  the 
sartple  and  the  slide  speed  were  also  recorded.  Therefore,  the  flame  heat  flux  can  be  plotted 
versus  the  distance  from  the  cooling  plate  to  the  measuring  station. 
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6.2  ANALYSIS  OF  EXPERIMENTAL  DATA 


Figure  6-1  shows  measured  temperatures  for  a  particle  board.  The  external  heat  flux 
at  425  mm  away  from  the  lower  edge  of  the  radiant  panel  was  20  kW/m^  (external  heat  flux 
has  a  distribution  as  described  in  Chapter  5).  The  preheating  time  before  applying  the  pilot 
ignition  flame  was  120  seconds.  The  time  to  start  moving  the  slide  was  147.  seconds.  GAS 
TCI  and  SUR  TCI  are  gas  and  surface  temperatures  at  the  first  measuring  station.  Rest  of 
the  curves  can  be  inferred.  It  can  be  seen  that  at  the  same  measuring  station,  surface 
temperatures  are  higher  than  gas-phase  temperatures  (though  it  is  somewhat  hard  to 
distinguish)  when  the  flame  is  not  close  to  the  measuring  station.  Gas-phase  temperature 
starts  to  exceed  surface  temperature  when  flame  approaches  the  measuring  station  and  it 
remains  so  when  the  measuring  station  is  immersed  in  the  flame.  The  reason  is  that  when  the 
flame  was  far  away  from  the  measuring  station,  the  fuel  surface  received  radiation  both  from 
the  radiant  panel  and  the  flame.  The  hot  surface  heats  the  gases  above  it.  Therefore,  a 
thermal  layer  developed  wherein  the  gas  temperature  was  therefore  lower  than  the  surface 
temperature;  when  flame  approached  the  measuring  station,  however,  the  gases  above  the 
fuel  surface  below  involved  with  the  flame,  the  gas  phase  temperature  became  higher.  One 
can  observe  that  there  is  always  a  sudden  increase  of  temperature  when  flame  approaches 
thermocouples.  We  can  employ  these  periods  of  sudden  increase  to  obtain  flame  conductive 
heat  flux  profiles  as  we  shall  discuss  in  Section  6.4.  We  can  also  observe  in  Figure  6-1  that 
the  temperature  always  starts  the  sudden  increase  at  approximately  same  temperature  (in  this 
case  it  is  about  270°C)  at  all  measuring  stations.  This  manifests  that  the  experiment  was 
indeed  steady  state  and  this  characteristic  temperature  (e.g.  Tg,j,=270°C,  note  that  this  is  not 
ignition  or  pyrolysis  temperature,  but  rather  it  is  the  temperature  rise  caused  by  external  and 
flame  radiative  heat  fluxes)  together  with  temperature  histories  during  the  sudden  increase 
period  can  be  used  to  derive  both  the  magnitude  and  length  scale  of  the  gas  phase  conduction 
(see  Section  6.4). 
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Figure  6-1  Temperature  histories  of  a  particle  board  test,  external  heat  flux  level:  20 
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Figure  6-2  Total  heat  flux  histories  of  a  particle  board  test,  external  heat  flux  level:  20 


The  total  heat  flux  measurement  is  illustrated  in  Figure  6-2  where  a  1.6  mm  diameter 
small  probe  flux  gage  was  used  at  station  2  and  a  9.1  mm  diameter  probe  was  used  in  station 
3.  The  total  heat  flux  measured  at  station  2  shows  similar  behavior  as  the  temperature 
measurement  demonstrated  in  Figure  6-1.  The  peak  heat  flux  is  about  70  kW/m^  from  the 
small  flux  gage  reading  when  the  flame  front  reaches  the  measuring  station.  This  heat  flux 
is  the  summation  of  external,  flame  radiation  and  gas  phase  conduction.  The  total  heat  flux 
measurement  at  station  3  does  not  show  a  large  size  peak  for  two  reasons:  1)  the  flame 
quenches  and  moves  around  the  25.4  mm  diameter  gage  surface  which  includes  the  detector 
surface  and  copper  cylinder  (sink)  (see  Chapter  5),  and,  2)  the  experiment  is  stopped  once 
the  flame  front  reaches  station  3  because  of  the  length  constraint  of  the  LIFT  (the  sample  can 
not  move  more  than  600  mm  because  the  leading  edge  of  the  sample  is  blocked  by  the 
structure  of  the  LEFT).  Flame  heat  flux  is  determined  by  subtracting  external  heat  flux  from 
the  radiant  panel  from  the  total  heat  flux  measurements  in  Figure  6-2.  The  result  is  plotted 
in  Figure  6-3  where  the  abscissa  has  been  changed  from  time  to  distance  from  the  edge  of  the 
cooling  plate  by  using  the  initial  location  of  measuring  station  and  the  slider  speed.  It  can  be 
seen  from  the  small  flux  gage  measurement  at  station  2  that  the  flame  heat  flux  (flame 
radiation+gas  phase  conduction)  is  about  50  kW/m^.  The  distance  at  which  the  peak  flame 
heat  flux  occurs  is  in  fact  the  width  of  the  burning  area  (if  there  is  no  burnout)  and  in  this  case 
is  about  100  mm.  This  length  is  consistent  with  video  recording. 

It  is  seen  in  Figure  6-3  that  flame  heat  fluxes  are  not  the  same  at  same  distance  from 
the  cooling  plate  edge,  namely  the  reading  of  small  flux  gage  is  much  larger  than  that  of  large 
flux  gage.  The  same  problem  exists  for  aU  tests.  The  possible  explanations  are  offered  in 
Chapter  5. 
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6.3  EXAMINATION  OF  THE  FLAMMABILITY  PROPERTIES  OF  PMMA 
SHEETS  AND  PARTICLE  BOARDS  USED  IN  THE  EXPERIMENTS 


Before  we  conduct  detailed  calculation  of  gas-phase  conductive  heat  flux  and 
compare  results  from  the  numerical  model  and  experimental  data,  the  flammability  properties 
of  PMMA  and  particle  boards  have  to  be  obtained.  The  properties  derived  in  earlier  work 
(Chen  et  al.,  1993;  1995)  are  used  as  inputs  for  the  heat-up  portion  of  the  horizontal  flame 
spread  model,  using  the  total  heat  flux  station  3  measured  by  the  large  heat  flux  gage  (as 
stated  before,  this  measurement  is  reliable).  The  calculated  surface  temperature  histories  are 
con:q)ared  with  experimental  measurement  of  station  3  surface  temperatures.  The  accuracies 
of  these  properties  are  therefore  checked  as  described  below. 

Figure  6-4  illustrates  the  surface  temperature  histories  of  a  19  mm  particle  board 
experiment  and  predictions  from  the  heat-up  submodel  of  the  HFSCS  code.  The  properties 
used  in  the  prediction  are:  k=0.25  kW/m.K,  p=765  kg/m^  and  c=2800  J/kg.K.  The  circular 
line  shows  the  experimental  data.  The  triangle  line  indicates  the  prediction  without 
convection  loss  while  the  square  line  illustrates  the  prediction  with  convection  loss  from  the 
fuel  surface.  The  convection  loss  is  evaluated  as  natural  convection  on  a  horizontal  plate 
(q^n=h(Ts-TJ).  The  latter  one  obviously  agrees  better  with  the  experimental  measurement. 
The  comparison  shows  that  the  prediction  using  convection  loss  as  agrees  both  with  the 
experiment  indicating  that  the  properties  used  in  the  prediction  are  correct.  We  have  used 
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Figure  6-4  Comparison  of  surface  temperature  histories  for  a  particle  board  test.  This 
comparison  is  used  to  check  the  values  of  material  properties  (k,  p,  c) 
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these  properties  to  predict  surface  temperature  histories  of  other  experiments.  The 
predictions  all  agrees  with  the  experiments.  These  results  are  presented  in  Appendix  G  where 
the  comparisons  of  surface  temperature  histories  at  the  3rd  station  between  predictions  and 
experiments  for  19  mm  PMMA  sheet  are  also  shown.  The  PMMA  properties  used  in  the 
predictions  are;  k=0.2798.  kW/mK,  p=12(X)  kg/m^  and  c=2200  J/kg.K.  It  should  be  pointed 
out  that  all  predictions  account  for  convection  losses  from  the  surface. 

6.4  ESTIMATION  OF  GAS-PHASE  CONDUCTIVE  HEAT  FLUX 


It  is  the  intention  that  the  validation  of  the  horizontal  flame  spread  model  (HFSCS) 
using  the  data  from  the  constant  speed  horizontal  flame  spread  experiments  (CSHFS)  would 
be  carried  out.  Before  doing  that,  the  gas-phase  conductive  heat  flux  profiles  (ref.  discussion 
of  Chapter  3),  has  to  be  quantified.  Gas-phase  conductive  heat  flux  profiles  have  been  shown 
both  experimentally  (Babrauskas  and  Watterlund,  1995)  and  numerically  (di  Blasi  et  al., 
1989)  to  follow  an  exponential  profile: 
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where  the  magnitude,  and  the  length  scale,  6g  have  to  be  obtained.  These  parameters  can 
be  derived  from  the  experimental  data  of  surface  tenperature  histories.  The  heat  flux  profiles 
are  essentially  the  same  shape  as  that  of  the  temperature  profiles,  if  solid  conduction  is 
assumed  to  be  one-dimensional  into  the  solid,  i.e.  conduction  parallel  to  the  fuel  surface  is  not 
important.  For  justification,  readers  can  refer  to  Appendix  E.  The  surface  temperature  rise 
of  the  fuel  at  any  location  for  thermally  thick  solid  can  be  expressed  (eqn.  (E-5)  of  Appendix 
E): 
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where  tp  is  the  time  to  pyrolysis  and  V  is  the  velocity  at  which  the  flame  front  travels.  T^^t  is 
surface  terrperature  rise  owing  to  other  heat  fluxes  (external  and  flame  radiative  heat  fluxes). 
It  can  be  easily  seen  by  comparing  equations  (6-1)  and  (6-2)  that  the  profile  of  surface 
temperature  rise  also  exponentially  decays  at  the  same  rate  as  the  gas-phase  conductive  heat 
flux  distribution. 


If  a  point  at  a  fuel  surface  starts  to  pyrolyze,  i.e.  Ts=Tp,  Eqn.  (6-2)  changes  to: 
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One  can  carry  out  sirrple  manipulation  by  inspecting  Eqn.  (6-2),  assuming  two  surface 
temperatures,  Tj,  and  Ts2,  at  two  different  times,  tj  and  t^: 
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dividing  eqn.  (6-4a)  by  eqn.  (6-4b),  one  has 
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If  we  assume  the  left  hand  side  of  eqn.  (6-5)  equals  e'*,  then  V(t2-ti)=6g.  Therefore 
we  can  chose  a  surface  temperature,  Ts2  during  a  sudden  increase  period  and  read  off  the 
time,  t2,  that  corresponds  to  temperature,  Ts2,  from  a  detailed  plot  of  surface  temperature 
during  a  sudden  increase  period  as  illustrated  in  Figure  6-5.  By  using  eqn.  (6-5)  with  RHS=e‘ 
*,  Tj,  can  be  evaluated.  One  then  reads  off  the  time,  t„  from  the  same  plot.  The  length  scale, 
6g,  can  finally  be  obtained  as  6g=V(t2-ti).  Obviously,  this  can  be  done  only  if  T^  is  a 
continuous  and  smooth  function  between  tj  and  t2.  Finally  one  can  solve  for  the  magnitude 
of  the  heat  flux,  (jg,  by  rearranging  eqn.  (6-3): 


(6-6) 


where  the  pyrolysis  temperature,  Tp,  is  obtained  from  literatures  (Chen  et  al.,  1995),  i.e. 
Tp=370°C  for  PMMA  and  Tp=320°C  for  particle  board.  The  pyrolysis  temperature  can  not 
be  determined  from  current  experiments,  but  rather  should  be  determined  from  a  separate 
pyrolysis  test  (Delichatsios  et  al.,  1991).  Again  it  should  be  pointed  out  that  this  solution  is 
only  applicable  for  thermally  thick  solid.  The  materials  tested  here  can  all  be  treated  as 
thermally  thick  as  far  as  gas-phase  conductive  heat  flux  is  concerned. 
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The  above  analysis  of  the  CSHFS  data  has  been  carried  out  and  the  results  are 
illustrated  in  Tables  6-1  through  6-3  (for  details,  see  Appendix  G).  In  Table  6-1,  8  Potlatch 
particle  board  (19  mm  thick)  tests  are  shown.  There  are  three  external  heat  flux  levels  at  425 
mm  away  from  the  radiant  panel  lower  edge  (see  Figure  5-2):  10, 20  and  30  kW/m^  (Column 
5  in  Table  6-1).  Three  slide  (flame  spread)  speeds  (0.32,  1.27  and  2.56  mm/s,  Column  4) 
were  used  corresponding  to  the  three  heat  flux  levels.  The  surface  temperatures  before  a 
sudden  increase,  Tg*,  are  listed  in  Column  6.  The  temperature,  T  ^2  is  320°C.  The  time 
differences,  'r=t2-ti  (see  eqn  (6-5))  are  presented  in  Columns  8  and  9  for  surface 
thermocouples  at  the  first  and  second  stations,  respectively.  For  each  test,  only  surface 
thermocouple  measurements  at  stations  1  and  2  are  used  since  the  experiment  was  stopped 
when  the  flame  front  approached  station  3.  Columns  10  and  1 1  show  the  magnitudes  of  the 
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Figure  6-5  The  detailed  surface  temperature  histories  during  the  sudden  increase  period 
for  a  particle  board  test 

heat  fluxes.  The  final  three  colunms  of  Table  6-1  list  average  characteristic  time,  length  scale, 
6g,  and  magnitude,  qg,  of  the  gas  phase  conduction  heat  flux.  There  are  considerable 
differences  among  6g  and  for  different  flux  levels  and  even  for  same  flux  levels. 

The  analysis  of  PMMA  test  data  is  shown  in  Tables  6-2  and  6-3.  In  6.35  mm  PMMA 
test  (Table  6-2),  0.5  mm/s  slide/flame  spread  speeds  were  used  at  the  nominal  heat  flux  level 
of  10  kW/m^.  The  heat  flux  levels  shown  in  Column  5  are  the  values  at  400  mm  from  the 
lower  edge  of  the  radiant  panel.  The  19  mm  PMMA  test  results  are  demonstrated  in  Table 
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6-3  where  the  slide  speeds  at  10  kW/m^  level  were  not  constants  (0.5, 0.2  and  0.35  mm/s) 
because  of  difficulties  of  obtaining  constant  burning  areas  in  the  experiments. 

It  was  postulated  that  the  gas  phase  conductive  heat  flux  profile  is  independent  of  fire 
size  (Atreya,  1984) ,  i.e.  for  the  same  material  and  ambient  conditions,  6g  and  (jg  should  be 
constants.  It  can  be  clearly  seen  that  the  results  are  inconsistent  and  greatly  scattered.  The 
wide  differences  led  to  the  speculation  that  the  derived  6g' s  are  really  the  length  scales  of  the 
gas-phase  conductive  heat  flux  profiles.  Careful  examination  has  revealed  that  the  6g  are 
essentially  the  length  scales  of  the  surface  temperature  distribution  along  the  fuel  surface! 
The  more  explanation  of  this  abnormal  phenomenon  will  be  discussed  in  the  next  section  and 
the  experimental  data  will  be  reprocessed  using  a  new  formulation. 


66 


— 

00 

ON 

<N 

m 

t—H 

VO 

Tt 

ON 

O 

Cv| 

cs 

o 

vd 

vd 

uS 

csi 

<pM 

tM 

CN 

1-H 

o 

VO 

00 

00 

cn 

00 

r-' 

00 

o 

00 

ON 

m 

m 

o“S 

ri 

cs 

'*d- 

cd 

VO 

r- 

t> 

ON 

1-H 

« s? 

tT 

00 

<N 

Tt 

cq 

m 

00 

iq 

00 

00 

r-H 

cvi 

d 

— 

r- 

o 

VO 

VO 

o 

cs 

in 

CNJ 

r- 

1— «< 

m 

o 

Os 

iri 

rd 

id 

d 

d 

cd 

vd 

f-H 

T>^ 

ro 

CN 

■■ 

VO 

VO 

ON 

ON 

in 

f— 4 

o 

ON 

00 

o 

o 

VO 

t> 

00 

1-H 

vd 

<d 

csi 

d 

id 

K 

rM 

Cn| 

cs 

cn 

n 

cn 

■ 

T— I 

in 

1? 

00 

ON 

00 

VO 

cn 

o 

D 

csi 

CM 

<N 

cq 

tq 

B 

1^ 

l-H 

cd 

d 

'cT 

VO 

r- 

in 

u 

ON 

»M 

vq 

rq 

rq 

cs 

cs 

c^’ 

r— ( 

csi 

d 

T3 

o 

o 

o 

o 

o 

o 

o 

o 

cd 

cs 

<N 

cs 

<N 

CM 

CN 

cs 

<N 

O 

H  w 

m 

m 

m 

m 

ro 

m 

cn 

cn 

nQ 

t-* 

00 

tq 

in 

cs  ^ 

B  O 

o 

00 

id 

ON 

Tt 

ON 

m 

ca 

t-- 

in 

00  p-H 

ed 

H  w 

cs 

(M 

cs 

<s 

CN 

CNl 

cs  in 

®  i"s 

VO 

VO 

rM 

00 

rf 

00 

X 

= 

l> 

oi 

t> 

<S 

vq 

c<i 

CN 

id 

d 

o 

m 

d 

00 

csi 

1 

<N 

cs 

CN 

CN 

c^ 

ON 

Cm 

CO 

O 

c 

cs 

cs 

<N 

oo 

00 

00 

00 

VO 

CO 

p 

ro 

m 

<N 

cs 

cq 

<s 

in 

«M| 

>  O 

o 

d 

d 

rM 

l-H 

csi 

CO 

o 

<N 

VO 

00 

cn 

1-H 

CN| 

(N 

(N 

1— < 

1—4 

CA 

VO 

VO 

VO 

VO 

vq 

00 

00 

'S 

o 

(U  1 

W 

w 

; 

M 

.1 

o 

o 

o 

o 

o 

o 

o 

o 

r5  i5 

r-H 

1-H 

cs 

<N 

(S 

cs 

cn 

'Hh  IZ 

H 

H 

H 

H 

H 

H 

H 

H 

*M 

P- 

X 

W 

o 

CO  X 

r"  ►§ 

cs 

Tf 

in 

cn 

1  t-' 

H 

H  ^ 

CN 

CN 

<N 

cs 

(N 

I 

VO 

■ 

■ 

o 

d 

H 

2 

B 

C>1 

B 

m 

VO 

00 

cd 


O 

c 


The  test  data  show  that  two  surface  temperatures  start  to  rise  quickly  at  two  different  values  indicating  that  flame  spread  was 


Table  6-2  The  Experimental  Results  of  6.35  mm  (1/4")  PMMA  Sheet 


^  N/A  here  means  that  either  surface  thermocouple  failed  to  work  or  the  data  fluctuate  too  much  to  read 


Table  6-3  The  Experimental  Results  of  19  mm  (3/4")  PMMA  Sheet 
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6.5  ENERGY  BALANCE  AT  THE  FLAME  LEADING  EDGE  FOR 
HORIZONTAL  SPREAD;  A  NEW  INTERPRETATION  OF 
EXPERIMENTAL  DATA 


In  an  attempt  to  explain  the  inconsistences  of  gas-phase  conduction  profiles  for  the 
same  materials  (cf  last  two  columns  of  Tables  6-2  and  6-3),  we  had  to  reexamine  the  physics 
of  this  process  in  detail.  A  schematic  diagram  of  this  process  is  shown  in  Figure  6-6  where 
different  heat  fluxes  near  the  flame  fi'ont  are  illustrated.  We  have  derived  gas-phase 
conduction  profiles  fi'om  surfece  temperature  histories  based  on  the  postulate  that  gas-phase 
conduction  profile  is  the  same  as  the  surface  temperature  profile  or  in  other  words,  only  solid 
conduction  normal  to  the  fuel  surface  was  considered  in  those  calculations  while  solid 
conduction  upstream  of  the  flame  front  (parallel  to  the  fuel  surface)  was  not  accounted  for. 
This  explains  why  the  length  scale  of  gas  phase  conduction  for  6.35  mm  PMMA  was 
generally  shorter  than  that  for  19  mm  PMMA,  since  the  path  of  conduction  is  limited  for 
thirmer  PMMA  (Wiliams,  1976).  There  is  also  a  uniform  heat  flux  applied  in  the  pyrolvzing 
region.  qS  which  also  contributes  to  the  heating  upstream  of  the  pyrolysis  front  because  it 
frees  som^ of  the  flame  heat  bv  assisting  in  the  vaporization  process.  The  new  energy  balance 
for  the  control  volume  of  Figure  6-6  is: 


or 


(6-7) 


where  6^  is  the  vertical  thermal  length: 


a. 


and  6ij  is  the  vertical  thermal  scale  for  the  surface  temperature  distribution: 

5  _  SX  S 


The  coefficient  for  each  term  in  equation  (6-7)  may  not  be  exactly  equal  to  one.  But  there 
are  not  much  different  from  one  (de  Ris,  1969;  Delichatsios,  1986).  e  ^  is  the  surface 
emissivity.  The  LHS  of  equation  (6-7)  represents  rate  of  energy  gain  by  the  control  volume; 
the  first  term  on  the  RHS  is  the  conducted  heat  from  the  flame  near  the  pyrolysis  front;  the 
second  term  contributes  to  the  heating  upstream  of  flame  front  and  q'^  is  defined  as: 

q^=netheatflux=q"^-aTp  (6-8) 
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the  third  term  is  the  contribution  from  the  preheating  of  external  heat  flux+flame  radiation. 
If  an  exponential  decaying  profile  of  this  heat  flux  is  assumed,  i.e.  q'l'exp  (|x|/fii),  the  surface 
temperature  rise  due  to  this  heat  flux  can  be  estimated  as: 


9sCfT-TW 


K  .  //» 

S 


(6-9) 


Tj  is  surfece  temperature  due  to  external  heat  flux;  the  last  term  of  the  RHS  of  the  equation 
(6-7)  is  the  reradiation  loss  for  the  control  volume.  One  can  combine  equations  (6-9)  and  (6- 
7),  substitute  arrive  at  a  new  energy  equation  concerning  the  control  volume  in  Figure  6-6 
can  be  derived  as: 


N 


6^  //  // 

cc  -e — —6„ 


(6-10) 


By  examining  equation  (6-10),  one  can  see  that  the  last  term  of  right  hand  side  is  small 
compared  to  other  terms  and  therefore  it  can  be  ignored  for  simplicity. 

Eqn.  (6-10)  can  be  used  to  derive  gas  phase  conduction  profile  from  the  experimental 
data.  The  approach  is  described  as  follows: 

1)  The  flame  spread  speed,  ¥„  of  the  experiment  is  known  exactly,  the  surface 
temperature,  Tj,  is  measured  in  the  experiment  (see  Fig.  (6-5)).  The  material 
properties  values  are  known  as  described  in  Section  6.3 .  qj  is  also  known  from  an 
experiment. 

2)  By  applying  eqn.  (6-10)  for  two  different  experiments  (i.e.  different  external  heat  flux 
levels  at  425  mm  away  from  the  radiant  panel  edge  and  different  flame  spread  speeds), 
two  equations  with  qg  and  6g  are  obtained.  Subtract  the  first  equation  from  second 
equation,  one  equation  with  only  6g  is  derived  since  the  first  term  on  the  RHS  of  eqn. 
(6-10)  is  eliminated.  The  remaining  algebraic  equation  can  be  simply  solved  for  8g,  q^ 
can  be  consequently  derived. 

We  have  applied  above  approach  to  the  experimental  data  discussed  above.  The 
results  in  Tables  6-1  to  6-3  are  reprocessed  and  shown  in  Tables  6-4  to  6-6.  Note  that  we 
have  used  averaged  values  of  surface  temperatures,  T,,  for  the  same  external  heat  flux  level 
and  flame  spread  speed.  The  derived  gas-phase  conduction  heat  flux  profiles  parameters  (q^ 
and  6^  are  consistent  (cf  the  values  of  8g  and  q^  in  Column  4  and  6  of  Tables  6-5  and  6-6): 
(^'“70  kW/m^and  6g“1.4mmforPMMA.  For  particle  board,  qg“30  kW/m^  and  6g“1.0  mm. 
The  consistency  of  the  results  for  different  thicknesses  of  PMMA  and  different  tests  for 
particle  board  shows  the  validity  of  this  analytical  approach.  We  will  use  above  values  of  gas- 
phase  conduction  in  the  flame  spread  model  to  simulate  flame  spread  and  numerical  results 
will  be  compared  with  the  experimental  data  which  are  discussed  in  the  following  section. 
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Finally,  an  error  analysis  is  carried  out  for  the  20  kW/m^  particle  board  test  which  is 
representative  of  e)q)erimental  errors.  These  errors  arose  mainly  from  two  variables,  surface 
temperature  and  flame  front  location  (which  would  cause  errors  in  determinng  the  heat  flux 
behind  the  flame  front,  q^).  There  are  many  factors  that  can  cause  errors  in  surface 
temperature  measurements,  but  two  factors  are  particularly  important:  1)  location  of 
thermocouple,  (I.e.  it  is  important  to  make  sure  that  the  thermocouple  does  measure  surface 
temperature  instead  of  gas-phase  or  solid  temperature  inside  the  solid)  2)  the  other  factor  is 
not  really  an  experimental  error,  rather  it  is  from  the  use  of  the  surface  temperature  at  the 
point  of  the  sudden  increase  in  the  above  analysis.  In  the  analysis  shown  above,  the  surface 
temperature  measurement  at  the  time  of  the  sudden  increase  is  used.  In  fact,  there  is  still 
temperature  rise  due  to  external  heat  flux  and  flame  radiation  during  the  sudden  increase 
period  which  was  not  accounted  for  in  the  above  analysis.  This  uncertainty  of  surface 
temperature  can  be  estimated  from  examining  the  detailed  plot  of  surface  temperature  during 
period  (e.g.  Figure  6-3).  For  20  kW/m^  particle  board  test,  the  uncertainty  is  estimated  to  be 
±7.5°.  Using  this  uncertainty,  the  magnitude  of  the  gas  phase  conduction  of  qg=27.63  kW/m^ 
is  determined  with  the  external  imposed  flux  of  20  kW/m^  (see  Table  6-4).  The  uncertainty 
of  flame  front  location  is  ±1  cm  from  experimental  observation  and  this  could  be  resulted 
from  fluctuations  of  flame  front  locations.  Accounting  for  this  uncertainty,  a  value  of 
q^-32.62  kW/m^  is  obtained.  Therefore,  we  can  determine  that  the  overall  uncertainty  in  qg 
is:  [(27.63-28.69)2+(32.62-28.69)2]‘^  kW/m^  or  about  14%.  Similarly,  the 

uncertainty  in  6g  is  about  21%. 


2  “external  heat  f/ux  to  pyroiyztng  solid 


Energy  balance  In  solid: 

pc  (T-T)V6mk 
s  s  p  »  S  V  a  Cp 

(solid  convoction)  (ftamo  conducffon) 


07^ 


(oxtomsi  Iluxos)  (rondiadon) 


Figure  6-6  A  diagram  of  heat  fluxes  at  the  flame  leading  edge 
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6.6  NUMERICAL  SIMULATION  OF  FLAME  SPREAD  AND  COMPARISON 
WITH  THE  EXPERIMENTAL  DATA 


After  determining  the  gas-phase  conduction  heat  flux  profile  (i.e.  q"g  and  6g)  from 
the  experimental  data,  the  following  steps  have  been  taken  in  order  to  determine  flame 
spread  using  the  horizontal  flame  spread  model  (HFSCS)  and  to  compare  with  experi¬ 
mental  data:  1)  the  gas-phase  conductive  heat  flux  profile  was  used  in  the  code  as  an 
imposed  heat  flux,  2)  the  flammability  properties  examined  in  Sec.  4.2  were  the  code 
input,  3)  because  the  large  heat  flux  measurements  were  reliable  when  flame  was  away 
from  the  large  heat  flux  gage,  the  external  and  flame  radiative  heat  fluxes  at  any  point  on 
the  fuel  surface  at  any  instant  were  inferred  from  the  large  heat  flux  measurement  for  each 
experiment.  This  was  possible  since  the  distance  between  any  point  on  the  fuel  surface  and 
the  cooling  plate  edge  at  any  instant  of  time  was  known,  and  assuming  that  the  external 
heat  flux+flame  radiation  were  the  same  for  the  same  distance  from  the  cooling  plate  edge 
(see  discussion  of  Fig.  6-3).  The  inferred  heat  flux  histories  are  used  in  the  predictions. 

Figures  6.7-9  illustrate  the  comparisons  of  the  measured  (via  video  recording)  and 
predicted  flame  front  locations  from  the  sample  leading  edge  for  the  19  mm  particle  board 
experiments  for  three  different  external  heat  flux  levels.  The  video  recording  was  used  to 
obtain  flame  front  position  at  discrete  time  intervals  using  frame  by  frame  analysis.  One 
can  see  that  all  predictions  agree  with  the  experimental  data.  There  are  large  differences 
shown  in  Figure  6-10.  This  would  probably  due  to  1)  the  flame  spread  was  not  steady  at 
the  beginning,  2)  no  attempt  was  made  to  model  the  exact  initial  burning  area  and  ignition 
length.  This  could  result  in  discrepancies  at  matching  the  initial  flame  front  location  as 
seen  in  Figure  6-7.  The  main  aim  is  here  is  to  compare  the  slope  of  the  flame  leading  edge 
position  vs.  time  in  Figures  6.7-9,  which  represents  the  flame  spread  velocity. 

The  same  approach  was  used  to  predict  experimental  results  at  three  different 
external  heat  flux  levels  for  the  19  mm  thick  PMMA  specimens,  shown  in  Figures  6. 10-12. 
Good  agreement  was  observed  except  for  the  case  M20K809  where  there  were  some 
problems  obtaining  a  steady  flame  spread  during  the  experiment.  Similarly  good 
agreements  can  be  seen  for  the  comparisons  for  the  6.35  mm  PMMA  thick  samples  shown 
in  Figures  6.13-15. 

Finally,  we  show  in  Figure  6-16  the  comparison  of  detailed  surface  temperature 
histories  during  the  sudden  increase  period  between  the  prediction  and  a  particle  board 
test  at  the  second  measuring  station.  The  temperature  history  profiles  of  prediction  and 
experiment  are  the  same.  The  surface  temperatures  increase  suddenly  at  almost  same 
value  (-270°)  for  both  the  predication  and  the  experiment  indicating  that  the  model 
follows  the  surface  temperature  histories  during  the  flame  spread  experiment  quite  well. 
The  flat  part  of  the  prediction  at  320°C  simply  shows  the  pyrolysis  temperature  because 
after  pyrolysis  starts,  the  pre-heating  model  is  no  longer  valid.  Therefore,  the  flat  part  has 
no  physical  meaning.  The  surface  temperature  comparison  for  a  19  mm  PMMA  is  shown 
in  Figure  6-17.  Similar  results  are  shown  in  Appendix  L. 
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Some  comments  concerning  above  prediction  and  the  flame  spread  code  should  be 
made  to  help  the  reader  to  better  understand  the  current  approach;  1)  the  use  of  heat  flux 
measured  at  the  station  3  arrives  from  the  fact  that  the  calculation  of  the  heat  flux  for  the 
experiments  would  be  difficult  due  to  the  presence  of  the  radiant  panels  (results  in 
additional  unknown  flux).  A  radiation  model  can  be  used  to  calculate  the  heat  flux  ahead 
of  the  flame  for  an  arbitrary  condition,  2)  for  any  other  materials  being  modeled,  the 
constant  speed  horizontal  flame  spread  experiments  have  to  be  carried  out  to  derive 
parameters  of  gas-phase  conductive  heat  flux  profiles,  q"g  and  8g.  These  parameters  then 
can  be  used  in  the  flame  spread  code,  3)  material  flammability  properties  have  to  be 
determined  apriori  as  input  to  the  model  and,  4)  the  completed  horizontal  flame  spread 
model  can  only  predict  the  flame  radiative  heat  flux  based  on  a  radiation  model,  5)  a 
lumped  model  (with  combustion)  can  include  vitiation  in  a  limited  manner. 

6.7  LIMITATION  OF  THE  HORIZONTAL  FLAME  SPREAD  MODEL 

The  horizontal  flame  spread  model  developed  in  this  work  can  handle  sudden  jump 
of  heat  flux  as  described  in  Section  3.4.  The  model  can  also  use  transient  external  heat 
flux  with  no  problems,  at  least  for  the  cases  that  have  been  tested.  The  transient  effect  of 
the  flame  growth  on  variation  of  flame  heat  fluxes  can  also  be  dealt  with  easily.  The 
model  has  been  designed  to  accommodate  the  effects  of  vitiation  (i.e.  changes  in 
combustion  and  flame  heat  fluxes).  The  model  has  not  been  tested  extensively  for  sudden 
decrease  of  heat  fluxes,  such  as  heat  flux  drops  due  to  suppression. 
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Figure  6-7  Comparison  of  predicted  and  measured  flame  spread  for  a  particle  board 
experiment.  External  heat  flux  level;  10  kW/m^ 
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Figure  6-8  Comparison  of  predicted  and  measured  flame  spread  for  a  particle  board 
experiment.  External  heat  flux  level:  20  kW/m^ 
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Figure  6-9  Comparison  of  predicted  and  measured  flame  spread  for  a  particle  board 
experiment.  External  heat  flux  level:  30  kW/m^ 
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Figure  6-10  Comparison  of  predicted  and  measured  flame  spread  for  a  19  mm  thick 
PMMA  sample  experiment.  External  heat  flux  level:  lOkW/m^ 
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Figure  6-11  Comparison  of  predicted  and  measured  flame  spread  for  a  19  mm  thick 
PMMA  sample  experiment.  External  heat  flux  level:  20  kwW 
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Figure  6-12  Comparison  of  predicted  and  measured  flame  spread  for  a  19  mm  thick 
PMMA  sample  experiment.  External  heat  flux  level;  30  kW/m^ 
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Figure  6-13  Comparison  of  predicted  and  measured  flame  spread  for  a  6.35  mm  thick 
PMMA  sample  experiment.  External  heat  flux  level;  10  kWW 


TIME  (sec) 


Figure  6-14  Comparison  of  predicted  and  measured  flame  spread  for  a  6.35  mm  thick 
PMMA  sample  experiment.  External  heat  flux  level:  20  kWW 


Figure  6-15  Comparison  of  predicted  and  measured  flame  spread  for  a  6.35  mm  thick 
PMMA  sample  experiment.  External  heat  flux  level:  30  kW/m^ 
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Figure  6-16  Detailed  comparison  of  surface  temperature  at  the  second  measuring  station 
during  the  sudden  increase  period  with  the  model  prediction  for  the  particle  board  test. 
External  heat  flux  level:  20  kW/m^ 


M20K809 

3/4"  PMMA  Sheet,  V=1  .27  mm/s 


300  310  320  330  340  350  360  370  380  390 


TIME  (sec) 

Figure  6-17  Detailed  comparison  of  surface  temperature  at  the  first  measuring  station 
during  the  sudden  increase  period  with  the  model  prediction  for  the  19  mm  thick  PMMA. 
External  heat  flux  level:  20  kW/m^ 
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6.8  CHAR  DEPTH  MEASUREMENTS  AND  ADDITIONAL  DISCUSSION 


The  material  charring  is  an  important  parameter  for  determining  the  material 
flammabihty  properties  and  consequently  the  flame  spread  prediction.  The  pyrolysis 
model  that  we  have  developed  as  part  of  the  horizontal  flame  spread  model  can  predict  the 
char  depth  for  the  material  being  exposed  to  external  flux.  Although  the  focus  of  our 
research  was  not  on  in-situ  validation  of  the  charring  model,  the  char  depth  was  measured 
for  all  the  particle  board  test  specimens.  The  char  depth  was  measured  every  5  cm  at  3 
positions  across  the  material  width  by  cutting  the  specimen  along  the  width  and  the  three 
char  depth  measurements  were  averaged.  The  data  was  collected  and  tabulated  using 
three  specimens  for  each  of  the  external  flux  levels  (10,  20  and  30  kWW)  and  are 
included  in  the  report  for  completeness  in  Appendix  M.  A  plot  of  the  data  is  shown  in 
Figure  6-18.  Tests  1-3  for  each  of  the  external  flux  levels  are  shown  here.  Clearly, 
variations  in  the  material  homogeneity,  initial  material  conditions,  test  conditions  and 
burning  area  size  would  affect  the  char  depth.  Never-the-less,  the  three  distinct  levels  of 
charring  is  obtained  for  each  external  flux  levels.  Furthermore,  the  largest  variation  in  the 
char  depth  is  about  20%  fi-om  the  mean.  The  scatter  seems  to  be  more  for  the  10  kWW 
case,  however,  the  percent  variation  is  no  greater.  The  char  depth,  as  expected,  is  highest 
for  the  lowest  flux  since  the  flame  spread  is  slowest  in  that  case  and  thus  the  burning  zone 
is  being  moved  under  the  cooling  plate  (i.e.  the  flame  is  quenched)  at  a  slower  rate. 


Particle  Board 


Figure  6-18  Char  depth  measurements  along  the  particle  board  specimens  exposed 
to  3  different  external  flux  conditions.  Three  tests  for  each  condition. 


We  have  shown  (Chen,  et  al.,  1995)  that  the  char  depth  can  be  characterized  by  a 
characteristic  char  conductance  length,  dc,  shown  below: 

k^T* 

dc  =  -  (6-11) 

Aq" 

o 

where  kc  is  the  thermal  conductivity  of  the  char  and  T*  is  the  maximum  surface 
temperature  attained  if  the  imposed  external  flux,  q” ,  is  equal  to  the  surface  reradiation 

o 

losses. 


The  data  presented  in  Figure  6-18  was  normalized  using  equation  (6-1 1).  The  char 
depth,  Xc,  was  averaged  for  the  three  tests  and  normalized  by  the  char  conductance  length, 
dc,  and  is  shown  in  Figure  6-19.  This  plot  shows  that  for  a  steady  state  flame  spread 
condition,  the  char  formation  is  a  strong  function  of  the  externally  imposed  flux  and  does 
not  depend  on  the  flame  fluxes.  Furthermore,  we  can  conclude  that  this  char  conductance 
length  is  a  reasonable  parameter  in  that  all  the  data  collapse  approximately  on  a  single 
curve. 
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Figure  6-19  Normalized  char  depth  along  the  specimen  for  3  externally  imposed 
conditions,  constant  flame  spread  experiments. 
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TR  lma£ing 


We  have  shown  temperature  data  in  this  report  obtained  from  thermocouples. 
Although,  the  measurements  were  relatively  reliable,  thermocouples  only  provide  local 
measurements.  An  alternative  measurement  approach  would  be  to  use  IR  imaging  for 
surface  temperature  measurements.  The  advantage  is  the  field  measurements  and  ability  to 
have  evaluate  the  overall  experiment.  The  disadvantages  are  the  flame  interference  with 
the  surface  measurements,  reflections  causing  errors  and  the  relative  uncertainty  in  the 
measurement  which  must  be  carefully  calibrated  and  dependency  on  the  surface  emissivity. 
In  addition,  the  cost  of  such  system  would  be  much  higher  than  using  thermocouples.  We 
used  a  demonstration  system  to  obtain  some  preliminary  data  for  a  flame  spread  over  a 
PMMA  sample.  Figures  6-20  through  6-23  show  the  thermal  images  of  the  sample.  The 
burning  zone  measurements  are  totally  unreliable  since  a  special  filter  would  be  needed  to 
reduce  the  flame  interference.  An  emissivity  of  0.9  was  assume  for  the  painted  surface  of 
the  PMMA.  The  burning  zone  is  clearly  distinguishable.  The  lateral  uniformity  of  the 
burning  zone  is  not  typical  and  was  improved  in  actual  tests.  The  last  figure  is  an  actual 
contour  plot  of  the  surface  temperature  generated  by  the  image  processing  software.  The 
PMMA  p5^olysis  temperature  of  330"  is  quite  closely  measured  here  despite  the  lack  of 
very  careful  calibration  and  relying  on  the  internal  calibration  of  the  system. 
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Figure  6-20  -  IR  Temperature  map  ("C),  PMMA  Test,  19  mm  thick,  10  kW/m^ 
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Figure  6-21  -  IR  Temperature  map  fC),  PMMA  Test,  19  mm  thick,  10  kW/m^ 
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Figure  6-22  -  IR  Image,  Surface  temperature  ("C)  contour  map  of  Figure  6-21, 
PMjMA,  19  mm  thick,  10  kW/m^ 


Figure  6-23  -  Surface  temperature  profile  of  Fig.  6-21  along  the  designated  line. 
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Chapter  7  Conclusions  and  Future  Work 


A  comprehensive  opposed  flow  horizontal  flame  spread  model  has  been  developed 
and  validated  by  comparison  with  an  analytical  solution  and  experimental  data.  The  model 
uses  measured  gas-phase  convective  heat  flux  distribution,  which  has  a  characteristic 
exponential  decaying  profile  with  a  length  scale  5g  and  a  maximum  flux,  q"^.  The  model 
includes  all  phenomenon  involved  in  horizontal  flame  spread  as  follows; 

1)  Solid  conduction  and  pyrolysis  for  both  charring  and  non-charring  materials  is 
calculated  using  a  one  dimensional  integral  method  validated  by  comparison  with 
experiments. 

2)  Combustion  of  pyrolysis  gases  and  radiation  from  turbulent  difliision  flame  is 
simulated  using  a  combustion  model.  A  phenomenolo^cal  soot  formation  model 
based  on  laminar  smoke-point  height  was  applied  in  turbulent  combustion  to  predict 
soot  concentration  and  flame  radiation.  As  a  first  step,  this  soot  and  radiation 
model  was  used  for  optically  thin  flames  in  order  to  derive  an  energy  equation. 
Comparison  of  predictions  with  radiation  measurements  in  turbulent  buoyant  jet 
flames  (See  Chapter  4)  for  different  firels  provided  a  partial  validation  to  the  present 
model.  For  pool  fires,  application  of  the  model  provides  encouraging  results,  but  it 
also  points  to  a  need  for  a)  fiirther  examination  of  the  flow  near  the  base  of  the  pool 
fire  b)  relaxation  of  the  optically  thin  flame  approximation  as  the  pool  size  increases 
and  c)  state  relationships  for  real  fuels. 

3)  Derivation  of  gas-phase  convective  heat  flux  from  measurements  of  a  constant  speed 
flame  spread  experiment.  A  quasi-steady  horizontal  flame  spread  experiment 
(constant  speed  and  burning  area)  was  designed  and  conducted  for  two  materials; 
PMMA  and  particle  board.  Sur&ce  tenq>erature  and  heat  flux  histories  were 
measured.  The  sur&ce  temperature  histories  were  analyzed  using  an  analytical 
solution,  to  obtain  the  profiles  of  gas-phase  convective  heat  flux.  It  was  deduced 
that  the  solid  conduction  from  flame  front  to  upstream  of  the  flame  front  has  to  be 
considered  (conduction  parallel  to  fiiel  sur&ce)  to  correctly  derive  the  gas-phase 
conductive  heat  flux  proffle.  A  new  analytical  approach  of  energy  balance  has  been 
developed  that  considers  contribution  to  flame  spread  from  external  heat  flux  in  the 
pyrolyzing  re^on  behind  the  flame  front  to  flame  spread  (See  Chapter  6).  Profiles 
of  gas-phase  convection  was  finally  determined  using  this  approximation.  The 
uncertainty  of  the  magnitude  of  the  gas-phase  convective  heat  flux  from  the  analysis 
of  the  experimental  data  for  a  representative  particle  board  test  is  about  14%  for 
q  g  and  21%  for  5g.  This  profile  has  itKX>rporated  effects  on  gas  phase  conductive 
heat  flux,  namely,  gas  phase  finite  chemistry,  complex  flow  near  the  flame  front  as 
well  as  ambient  condition.  The  d^ved  profile  was  finaUy  used  in  the  horizontal 
flame  spread  model  to  calculate  flame  spread.  The  model  has  been  validated 
through  experimental  data  of  flame  spread  and  surface  temperature  histories. 
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The  future  work  can  include  the  following  items: 


1)  Complete  coupling  of  the  horizontal  flame  spread  model  with  the  integral 
combustion  model  should  be  implmented.  Such  work  would  involve: 

a)  validated  model  of  heat  flux  to  the  burning  surface,  b)  measurements  of  state 
relationships  for  arbitrary  ftiel,  c)  improving  numeric  of  the  combustion  model,  and 
d)  extension  of  the  radiation  sub-model  to  optically  intermediate  and  thick  flames. 

2)  Experiments  could  be  performed  for  other  materials  (e.g.  wood)  to  derive  gas 
phase  conduction  heat  fluxes. 

3)  Insertion  of  the  horizontal  flame  spread  model  into  a  compartment  models  (e.g. 
CFAST)  should  also  be  carried  out. 

4)  The  flame  spread  should  be  tested  extensively  for  sudden  decrease  of  heat  fluxes. 
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Appendix  A  An  Integral  Model  for  Prediction  of 
Radiation  and  Combustion  of  Turbulent 
Diffusion  Flames 

1.  The  mean  species  concentration  and  temperature  calculations 
Fluctuation  correlation^*^: 

_  1 

1+Y - = — 

1-^ 


where 

5 :  local  mean  mixture  fraction  and 
mean  centerline  mixture  fraction 
The  parameter  X  used  in  T-function  therefore  is: 

la-I)  ^  l-^3x0.1M5l^-O  ^ 
"  0.16(5l^-I) 


(A-1) 


for  a=0.16,  P=5  and  y=3 
^  is  defined  as: 


h-h 

CO 


h 


(A-2) 


The  method  or  approach: 

At  each  height  of  flame,  z,  one  can  calculate  the  entrainment  rate  by  using 
appropriate  entrainment  equations  (e.g.  whether  it  is  momentum  or  buoyancy  driven  flow). 
Then  is  calculated  as: 

(A-3) 


where  ihf  is  the  flow  rate  at  the  nozzle.  The  assumed  mixture  fraction  profile  is  used; 


A-1 


where  i|;  is  the  so-called  stream  function  and  is  given  by: 


r 


0 


Having  known  one  uses  eqn.  (1)  to  obtain  X.  The  probability  function  P(^,  ?  '^)  can 

therefore  be  calculated  by: 


r(X) 


r{n)r{m) 


m-1 


where  n=A  ^  and  m=A,(l-  ^). 

Given  we  can  use  Sivathanu  and  Faeth's  generalized  state  relationship^^^  to  calculate  the 
laminar  species  concentration  (see  the  relevant  part  of  the  paper). 

The  average  species  concentration  can  be  calculated  by: 

_  _  1  _  _ 

0 


while  the  laminal  temperature  can  be  easily  evaluated  by  using  equation  (A-2),  noting  that: 
hi-:  enthalpy  at  the  source 
h„:  enthalpy  of  the  ambient 
h:  enthalpy  in  the  flame 

(A-4) 


h  can  be  obtained  if  the  laminar  species  concentration  and  temperature  are  known,  which  is 
given  by: 


A-2 


(A-5) 


N  W  r 

■^n 


Y.C  .(T)dT 

i  pi  ' 


where 

hfii  heat  of  formation  of  species  i  at  temperature  Tq  (usually  298  K) 

T:  laminar  temperature 
Cp  (T):  heat  capacity  of  species  i 

On  the  other  hand,  knowing  h  (i.e.  from  RHS  of  eqn.  (A-4)),  one  can  evaluate  the 
laminar  temperature  by  trial  and  error.  It  is  straightforward  that  the  average  temperature 
T  ( O  can  be  calculated  as: 

_ 1  _  _ 

(A-6) 

0 


2.  Centerline  velocity  calculations: 

Using  the  momentum  equation  derived  by  Delichatsios,  et.  al^^^: 


dz 


(A-7) 


this  equation  can  also  be  written  as: 


2  ^  dz  2  ^  dz 


AT 
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Or 


dz 


i|/^  dz 


2  RHS 


A- 3 


Or  in  terms  of  u^: 


2  ARHS 

dz  1|J^  dz  ^  llJj. 


This  equation  takes  the  form  of 

y'^Pix)  y=Q(x) 


The  solution  of  this  type  of  equation  is: 


Or 


y=AjQ{x)  e^^^''^‘^dx^CA=Af{x)  ^CA 


where 


A=e 


-/P(x) 


dx 


f(x)=j'Q(x) 


Hence 


0  x-xl 

y^.A^fU)\^,;CA^  @  x-x^ 
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(A-9) 


(9)xAi-(8)xA2: 


^2^1  ^l'^2  -^1-^2  ^  lx=X2  Ix-xl^ 


e/"‘^’^dx 


Therefore: 


In  our  case: 


i|;^=0.018p„z 


5/2 


AH 


\  (5.1)  C  T 


hence: 
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•  dz 


And 


A- 5 


(Noting  @ 


The  first  term  can  be  calculated  by  using  normal  integration  procedure.  The  second  term  has 
to  be  handled  specifically  as  follows: 
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Assuming  R=0.999,  then 
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Therefore: 
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For  buoyancy  driven  flow: 

i|jt=az+ihf 
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3.  Generalized  state  relationship  from  Sivathanu  and  Faeth^^^: 
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f  CO  co^s 


(A-15) 


9m-MY^  Y 

ijj  =  ( - ^ ^ ) 

^  \os 


(A-16) 


One  can  do  some  simple  algebra  transformation  to  obtain  the  species  concentration  Y  as  a 
function  of  the  state  function  il/j; 

Y-fOl/i) 

N^: 
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i.e. 


Y  -Y  Y 

Y  IF*’  Y 

N^S  N^ 


■Y.,  J\J^ 


Y  llJ 

N^S'^N^ 


Therefore: 


Y  -til  Y 

rj  hj 


^2'“ 


(A-17) 


^2- 


32n+8/n+M,Y-  oo  q,v, 

f  0,,“  3Z+Bin 

- - +M, 


T  ih 

02'“^02 


{32n*Sm)  T  tIj 

^Of  **  ' 


'2'“  ^2 


32 n^Sm-^M^Y^  -M^Y^  }\!^ 

f  o^,-  f  o,,»^a 


2' 
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CH^ 


(A-19) 


COj: 


44n-M  Y  AAn-M  Y 

f  C02S_  ^^f^C02 

r  Ilf  ~  Y 

C02S^C02  COj 
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(A-20) 


H,0: 


(A-21) 


CO: 


H,: 


44/5y„  i|;^^ 

CO^s^  CO 


'CO 


44n-M^Y:,,  MX,, 

f  COjS  f  co^s^co 


(A-22) 


Y  = 


f  E2OS  f  S^OS^E^O 


{A-23) 


For  hydrocarbon  reaction: 

CH^^a  ( O^^kN^)  -^hCOy^  {m/2 )  H^O^ia-n-m/A )  0^-^akN^ 


V.M. 

r  _  11 

M^-^a  {32*28 k) 


where 

V;:  the  number  of  moles  of  a  product  species,  i,  per  mole  of  fuel. 
At  stoichiometry: 
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(n+m/A)  kM 

A)  {32^28 k) 


nM, 


co^ 


'CO^S 


M^+  ( n+m/  A ) 


(32+28A:) 


^ 

( n-m/ A)  (32^28 k) 


For  CH4  fuel,  n=l,  m=4. 

Y„^f0.12il  i'^^,.0.1514  S'^<,.0.1239 


The  following  is  the  correlation  of  mass  fraction  with  fiiel-air  equivalence  ratio. 
Species  mass  fraction:  Yj 
Fuel-air  equivalence  ratio,  (j) 

From  S  and  Fs  paper^^^  assuming  x=log<|),  y—logiljj 

1). 

yN2=1.446xl0-^-1.0x 
@  ^=0,  (i)=0,  Yn2=0.767 
^=1,  (j)=co,  Yn2=0 

2) .  Y02: 

(j)<l,  i|io2=0.9898-0.9062(j) 

4»1.  yo2=-0-79481-3.33465x+7.74823  x^-8.58911  xM.31337  x^-0.810195  x^ 
@  ?=0,  Yo2=0.233 
^=1,  Yo2=0.0 

3) .  Fuel-CH4 

4)^l,Yp„„=0 

4»1.  yFuei=-2.47524+3.0703  x-2.04325  x40.406869  x^ 

@  ^=0.  ypuePO- 
ypuel-l-O 

4) .  Yco2: 

4><l>yco2=-0.0332629+0.977861  x 

yco2=-0.092255+0. 104665  x-0.512132  x^-hO.l  10464  x^ 

®  ^=0,  yco2=0 
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C=l,  Yco2=0 

5) -  Yh20- 

yH2o=-0.0332629+0.977861  x 
4»1,  yH2o=-0.032632+0.171055  x-0.301843  x^ 

@  ^=0,  Yh2o=0 
^=1,  Yh2o=0 

6) .  Yco 

yco=-l. 00413+2. 103 14  x-2.80715  xVl. 19247  x^-0.180829  x^ 

@  ^=0,1,  Yco=0 

7) .  Yh^: 

@  ^=0, 1,  yH2~0. 

yH2=-2.08 179+2.7 159  x-5.80981  x45.59833  x^-2.57435  x'‘+0.445662  x^ 


4.  The  entrainment  correlations: 
1).  Buoyancy  driven  flow: 


f./p 

— = — =0.018  (z/len) 

{MjpJ 


3/2^ 


(M/p.)  ""2 


(a-24) 


M;,:  initial  momentum 


00 


dj:  the  nozzle  diameter 

u^:  initial  velocity  of  gas  flow 

therefore: 


A- 12 


Ar,  ah  (x,-x«) 

T  '  C  (S-.1)  T 

oo  p  oo 


if  z/len>5  use  correlation  in  buoyancy  driven  flow  region  as  above 
if  z/len<5,  momentum  driven  flow 


yp. 


=  .2012. 


yp.. 


5.  Radiation  losses  in  the  combustion  model 
Energy  balance: 


dz 


(A-25) 


h :  enthalpy  per  unit  mass  at  certain  location  (height) 
ho:  enthalpy  per  unit  mass  of  the  ambient 


Stream  function 


Q':  radiation  losses  per  unit  height  and  it  is  given  by: 


^R,  T^ch 


) 


{A-26) 


Xrt-  total  turbulent  buoyant  radiation  fraction  (for  methane,  Xrt=0.18;  for  propylene, 
Xr;,=0.4OI‘'I). 

Flame  height  calculations: 


A-13 


13.5F2/5 


^  s 

1 

where 

S:  Stoichiometric  ratio  of  air  to  fuel 


P’2  ^ 

1 

^rf 

(S.l)3  AHJl-X,) 

1 

Uj:  velocity  at  the  nozzle 
d^:  nozzle  diameter 
pji  density  of  fuel  at  the  nozzle 
subscript  0:  ambient  condition 
Xr=^Xr,t 

Qch-  tot^l  heat  release  rate=mf  AH^. 

The  f(z/Lf)  in  eqn  (A-26)  can  be  found  from  GRI  report.  By  using  eqn  (A-25),  one  gets: 

_  _  z 

j' (h-h^)  d[]f-ih-h^)ih^-j‘Q^dz  {a-27) 

0  0 


where  (¥-¥o)ihf  is  enthalpy  at  the  source. 

From  this  equation,  one  can  determine  the  temperature  of  the  laminar  flamelet. 
In  reference  [4],  one  has: 


q' 

(O') 

'  ^r'  max 


(A-28) 


Or  by  integration,  one  gets: 
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Or 


Q^-jQ'dz.{Q'^)^ff{z/L,)dz 

if 


Or 


Q  =  {Q')  Lr 

max  f 


where 

if 

ojf{z/L^)d(z/L^)  a  constant 
0 


But 


q/xa,=(q;)_-l,c 


So  that: 


(O') 

r  max  q 


By  inserting  this  eqn  into  eqn  (A-26),  one  has: 


f[z/L^) 


Thus: 
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QL^  f{z/L^) 

Y  Q  C 

^R^ch 


where  f(z/Lf)  is  the  same  as  that  in  Fig.  17  of  reference  [4],  and 

1 

C-j  f{z/L^)  d{  z/L^  =jf{w)dw 

0 


The  idea  for  incorporating  this  into  the  complete  model  is  as  follows: 

a.  given  different  one  can  obtain  the  laminar  adiabatic  temperature,  Tl^j. 

b.  assuming  the  linear  relationship  between  laminar  radiative  temperature  Tl  and  Tt 


T-  T  =a  ( 

L  a  ' 


) 


a  is  assumed  to  be  a  constant  for  each  height  and  in  fact  is  what  we  are  looking  for. 

c.  using  ^  and  ^  the  PDF  is  obtained: 

d.  calculating  average  temperature: 

_ 1 

T{^)=fT{^)PDFd^ 

0 


e.  using  eqn  (A-27),  the  a's  can  be  evaluated  (by  trial  and  error). 

The  function  f(z/Lf)  is  obtained  as  follows: 
assuming  x=z/Lf 
when  X:^ 0.387637 

y=f(z/Lf)=-1.42573E- 17+1. 80933  x-21.1433  x^=141.767  x^-212.667  x'^ 
when  x>0.387637 

y=-3.92171+28.685  x-54.3652  x2+35.801 8^-5.92926  x'^ 

6.  Calculation  of  the  radiation  constant  Cj 
Radiation  through  elemental  volume  of  flame: 

27101^  K;  dr 
Energy  balance: 
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dz 


00 


One  also  knows  for  optically  thin  flame: 


GT^K  =>GT'‘r=noT 

s  s  FI 


i-is 


The  soot  formation  is  assumed  occurred  somewhere  between  ^5  and  ^=1. 
Taking  PDF: 


jfaT‘ 


1-5. 


FDFd^ 


Assuming  ^=Y  where  y  is  a  number  larger  than  1. 
Therefore: 


lift 

Ql(Z)^cJ  ifoT^ 

0 


j^PDFd^) 

S 


Up 


And  at  the  flame  height: 


e,(L,)-cJo;(z)dz.xA 

0 


(A-29) 


Knowing  Xr,  and  Iq  Q'(z)dz,  one  is  able  to  solve  for  Cj. 

7.  The  simplified  combustion  model 

If  there  is  no  radiation  loss,  the  energy  balance  for  combustion  process  is: 


A-17 


(A-30) 


with  radiation  taking  into  account: 


dQJ.z) 

dz 


{A-31) 


where  T  is  the  actual  temperature  after  radiation. 
Assuming  Cpi  is  constant,  eqn.  (A-31)-eqn.  (A-30): 


dz 


dQ^iz) 

dz 


where 

Tj  is  adiabatic  temperature  and  T  is  real  temperature: 

T-T„=a(T,-TJ 

Therefore: 


d 

dz 


Vt _ 

Jsr,Cp.(l-a)  r  dilf. 

0 


dQ^iz) 

dz 


As  a  is  constant  for  certain  height  z. 


d 

dz 


_ 

(l-a)  Jsr.Cp.r  dii;. 
0 


dQ^(z) 

dz 


Performing  integration: 

Vc _ 

(l-a)  j'SY.C^.T^dllr.Qjz) 

0 


(A-32) 


On  the  other  hand,  we  know: 


A-18 


_  'ft  _ 

0  0 

where  Q^ht  is  the  enthalpy  at  the  source. 

Therefore: 

'I't _  _ 

0  0 


Qjh(z)  is  the  chemical  heat  release  rate  at  height  z. 

One  knows  and 

Vt _ 

|sy,c„r^diu 

0 

Qch(z)  can  be  calculated  easily: 

Vt _ 

/r,c„rdilr.d,,(2) 

0 


Thus  from  eqn.  (A-32),  one  has: 

QAz) 


1-a- 


(A-33) 


8.  The  combustion  model  using  radiation  fraction  Xr: 

It  can  be  guessed  that  the  left  hand  side  of  eqn.  (A-29)  can  be  simplified.  The  LHS  of  eqn. 
(A-29)  can  also  be  replaced  using  eqn.  (A-33).  Eqn.  (A-29)  can  be  changed  to: 


d 

dz 


[(l-a)c^^(z)]. 


dij; 

pE 


Suppose  one  is  able  to  obtain  C  from  total  radiation. 
Assuming 
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Then: 


d 

dz 


[(l-a)!3^^(z)]. 


CRjz) 

(1-a)" 


Assuming:  Y=(l-a)Q(.h(z),  then: 

^  CQ^^(z)Riz) 
dz 


Or: 


1  dY^ 
4  dz 


-CC^^(z)i;(z) 


(A-34) 


Or: 


r/-r/-4  cjQ^^{z)R{z)dz 


{A-35j 


At  flame  height: 

{1-Ot)  Q^^{z) -total  radiation  loss-m^H^Xg^ 


Thus  one  can  obtain  C  from  eqn.  (A-35). 

From  eqn.  (A-35),  one  has: 

22 

r/.r^Mc|c^^(z)i?(z)dz 


Or: 

^2 

1-a-  [  Y^*4  cfQ^l(z)R(z)dz]^'^/Q^Jz) 
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One  can  transform  R(z)  to  more  simple  form  as  follows  using: 


and 


Therefore: 


di); 

pu 


cfi); 


CO  S 


pu 


where 


r(X) 

—  £  ^  ^ (IE  ) 

r(n)r(m) 


References 

1.  Chatwin,  P.C.  and  Sullivan,  PJ.,  J.  Fluid  Mech.,  212,  p.  533, 1990. 

2.  Delichatsios,  M.A.  and  Mathews,  M.K.,  Second  International  Symposium  on  Fire  Safety 
Science,  1989. 

3.  Sivathanu,  Y.R.  and  Faeth,  G.M.,  Combust.  Flame,  82:21 1-230(1990). 

4.  Delichatsios,  M.A.,  Markstein,  G.H.,  Orloff  L.  and  de  Ris,  J.,  Gas  Research  Institute,  GRI- 
88/0100,  1988. 


A-21 


Appendix  B  Scalar  Fluctuations  in  a  Turbulent  Buoyant 
Jet  Flow  and  Comparison  with 
Experimental  Data 

The  calculation  of  scalar  fluctuations  is  essential  in  modeling  turbulent  flows  and 
combustion.  A  direct  correlation  for  the  fluctuation  of  a  conserved  scalar  inspired  by  Chatwin 
and  Sullivan  (1990)  was  used  and  incorporated  in  our  recent  work  for  a  turbulent  combustion 
model  (Motevalli  et  al.,  1992).  The  direct  correlation  of  root-mean-squared  fluctuations  of 
mixture  fraction  is  (Chen  et  al.,  1993); 

_ 1 

a?(pf^-0  a(p!^-I)  (B-1) 

1+Y - - 


where  a,  P,  y  are  parameters  determined  from  the  experiments  by  Koylu,  et  al.  (1990)  and 
and  ^  are  the  centerline  and  local  mean  mixture  fractions  respectively. 

A  physical  interpretation  of  eqn.  (B-1)  is  as  follows:  Consider  turbulent  mixing  at  a 
given  distance  from  the  origin  of  the  jet  where  the  centerline  value  for  the  mixture  fraction 
equals  One  may  visualize  that  mixing  at  this  height  occurs  by  large  scale  eddies  produced 
due  to  concentration  gradients  between  a  concentration  multiple  of  the  centerline  value,  i.e., 
Pl^  and  the  concentration  of  ambient  air  (where  ^=0).  If  molecular  mixing  is  absent,  then 
the  fluctuations  of  the  scalar  would  be  (Chatwin  and  Sullivan,  1990): 

(B-2a) 


To  account  for  the  molecular  mixing,  we  introduce  another  coefficient  a,  so  that  eqn.  (B-2a) 
becomes: 

_  _  ^RHS^Fcnd)  (B-2b) 

We  have  chosen  the  right  hand  side  of  eqn.  (B-2b)  to  be  a  function  of  the  scalar,  (and  not 
a  constant  as  it  was  done  in  Chatwin  and  Sullivan  (1990))  for  two  reasons:  a)  to  match  the 
mixing  conditions  at  initial  stages  near  the  nozzle  and  b)  to  preserve  the  condition:  ^  ^  (1- 
^),  which  must  always  be  satisfied.  A  simple  form  was  selected  as  shown  in  eqn.  (B-1)  where 
Y  must  always  be  greater  than  one.  Note  that  near  the  nozzle  where  ^  =  1  eqn.(B- 1)  becomes: 
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(B-3) 


l'\  1 

5(1-5)  Y 


For  mixing  in  non-combusting  jets,  Chatwin  and  Sullivan  (1990),  have  suggested,  after 
careful  corrparison  with  experiments,  the  following  values:  a=0.16,  P=1.3  while  the  RHS  of 
eqn.  (B-2b)  is  equal  to  one. 

In  this  work,  we  compare  assumed  scalar  fluctuations  calculated  from  different 
combinations  of  parameters  a,  p  and  y  (no  virtual  origin  was  used)  with  the  experimental 
data  given  by  Koylu  et  al.  (1990).  Koylu  et  al.  (1990)  measured  root-mean-squared- 
fluctuations  of  mixture  fraction  from  a  turbulent  diffusion  flame  generated  from  a  mixture  of 
Carbon  Monoxide  (CO)  and  Hydrogen  (Hj).  The  experimental  data  discussed  hereafter  are 
for  a  nozzle  diameter  of  7.04  mm,  CO  flow  rate  of  691  mg/s  and  H2  of  14.3  mg/s.  This 
conditions  represent  a  turbulent  buoyant  jet  flame  except  very  close  to  the  nozzle  where  the 
flow  is  momentum  dominated  (Delichatsios  et  al.,  1988). 


The  stream  function,  Y  is  related  to  radial  distance,  r  as: 


(B-4) 


where  u  is  mean  velocity,  p  is  mean  density  and  is  taken  as  a  constant  ( p=(  pg+  p„)/2,  average 
density  in  this  height)  in  the  following  calculation  and  u<,  is  centerline  velocity.  In  the 
calculation  presented  here,  the  velocities  from  the  experiment  (Koylu,  1990)  at  corresponding 
heights  are  used. 

Figures  B-la,  B-lb  illustrate  comparisons  of  the  experimental  data  and  the  results 
from  eqn.  (1)  with  a=4.5,  P=1.4  and  y=1  at  different  heights  from  the  flame  base  (x/d=40 
and  x/d=50).  The  abscissa  is  the  normalized  radial  distance  and  the  ordinate  is  the  scaled 
root-mean-squared  fluctuation  normalized  by  the  corresponding  centerline  mean  value.  Two 
numerical  results  are  shown  in  Figures  B-la,  B-lb:  1)  using  calculated  centerline  mean 
mixture  fraction,  (Motevalli  et  al,  1992)  and  2)  using  the  measurement  of  given  by 
Koylu  et  al.  (1990)  (therefore,  to  be  consistent,  Yt=2rnf/^c  is  incorporated  and  ihf  is  fuel  flow 
rate).  The  value  calculated  for  the  centerline  mean  conserved  scalar,  does  not  agree  with 
the  experimental  values  at  the  same  distance  from  the  nozzle,  because  our  integral  model  does 
not  represent  the  flow  conditions  near  the  nozzle  well:  no  adjustment  for  virtual  origin  was 
used  in  the  present  case. 

AU  calculated  results  show  higher  fluctuations  than  the  experimental  data  indicating 
that  adjustment  of  parameters  a,  P  or  y  in  eqn.  (B-1)  is  necessary  in  order  to  match  with  the 
experimental  data.  Therefore,  we  have  changed  y  to  be  2,  a  and  P  being  the  same  as  above. 
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to  calculate  fluctuations  using  eqn.  (B-1).  The  results  are  shown  in  Figure  B-3  where  the 
experimental  data  and  calculated  results  of  fluctuations  using  experimental  value  of  Tc  are 
plotted.  Good  agreements  are  observed  in  Figure  B-2.  The  same  conclusion  can  be  obtained 
from  Figure  B-3  where  the  comparisons  of  the  experimental  data  and  calculated  results  for 
a=3,  p  and  y  being  the  same  as  in  Figures  B-la,  B-lb,  are  shown. 


Figure  la  Comparison  of  fluctuations  for 
a=4.5,  P=1.4  and  y=1 


a=4.5,  p=1.4  and  y=1 


Figure  2  Comparison  of  fluctuations  for 
a=4.5,  P=1.4  and  y=2 


a=3,  P=1.4  and  y=1 


The  conparisons  shown  above  are  promising  in  that  one  can  select  a  set  of  parameters 
a,  p  and  y  in  order  to  use  a  direct  general  correlation  to  represent  experimental  reality. 
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It  is  decided  to  check  the  sensitivity  of 
selected  conserved  scalar  fluctuation 

relationships  on  predications  from  our 

integral  model.  Two  of  the  relationships 
examined  in  Figures  1, 2, 3  were  used  for  this 
purpose:  a)  a=4.5,  P=1.4,  y=1  (Figure  1) 
and  b)  a=4.5,  P=1.4  and  y=2  (Figure  2). 
The  correlation  for  case  b  (Figure  3)  agrees 
better  with  the  experiments  than  the 

correlation  for  case  a  (Figure  1).  It  should  be 
pointed  that  the  overall  predictions  for 

combustion,  radiation,  and  products 
generation  rates  based  upon  the  fluctuation 
relationships  we  have  proposed  above  do  not 
vary  significantly.  Figure  4  illustrates  the 
con^arison  between  the  data  from  Tamanini 
(1982)  and  numerical  calculations  of  cases  a 


1.00  -q - 
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Figure  4  Comparison  of  radiation  fraction, 
Xr  for  different  fluctuations 


and  b.  Little  difference  between  predictions  for  the  radiation  fraction,  Xr  was  found  between 


Y=1  and  y=2  indicating  the  validity  of  use  of  the  fluctuation  correlation  shown  in  eqn.  (B-1). 


We  can  use  any  one  of  the  above  fluctuation  relationships  in  the  integral  model  for  the 


combustion  process. 
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Appendix  C  Application  of  the  Combustion  Model  to 

Pool  Fire 


Delichatsios  (1988)  has  derived  the  entrainment  correlations  for  a  pool  fires  as: 
1).  For  z/D<  1.0,  close  to  the  pool  surface: 


m 


-F^^.0.086(  — 


(C-1) 


2).  For  1.0<z/D<4.0,  around  and  beyond  neck-in  area: 

— .-0.093  (  — 

(S*l)jn^  D 


{C-2) 


3).  Further  downstream  but  before  the  flame  tip: 


jn 


|i^F^^-0.018  (4)^"'^ 


{C-3) 


where  the  Froude  Number  is  defined  as: 


■  rf 


AH  , 
P. - -D^ 


S*1  \ 


AH 

(S*l)  C  T  ^  ^ 


and  the  theoretic  heat  release  rate: 

Q=ihfAH, 

The  momentum  equation  for  the  flame  is  (Delichatsios  and  Mathews,  1989): 


(^) 


d  r  1  ,  n  r 

—  [-rU^ri  ]./g ^ 

dz  2  J  u 


-dl|; 


This  equation  can  be  transformed  to  (Motevalli  el  al.,  1992): 
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ARHS 


du^  2  , 

C  ^  ^  t  „  2 

- + .  .  U  - - 

dz  i|;^  dz 


If  we  specify; 
y=Uc 


P{x) 


2 

ft  dz 


and  Q(x)=4  RHS/f, 

The  solution  for  the  momentum  equation  is: 


where; 


A-e 


-j'P{x)dx 


For  z/D<L0 


llj  -in  +in  -m  _+0 . 08  6  {5+1 )  in.  (  — ) 

^  t  ent  f  f  '  ^  D 


(C-4) 


where; 


0.086  (S+l)2n, 


a,.- 


1/2 


For  1.0<z/D<4.0: 


f 


3/2 


where; 


C-2 


0.093  (S*l)jn, 


3/2 


For  z/D>4.0: 


where; 


0.018  (S*l)mj, 


For  z<  1.0: 


p(x). 


2  _ 


2  11 
■3, 


1/2 


and 


fp(x)dx-f- 


1/2  \  .,1/2 


■dz 


(m^*a^z^'^)  z 


Assuming  x=z‘^,  then  x^=z.  Therefore: 


j*P(x)  dx-j'- 


a^2x 


(m^*a^x)  X 


-dx 


■/ 


2a, 


m^^+a^x 


-dx 


-21n  (m^+a^x)  -21n  (ni^+a^z 


and 


A~e 


-fp(x)dx 


1/2  \  2 
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Therefore: 


For  1.0<z/D<4.0: 


P(x). 


o  3a 

3^1/2  "^^2^ 


m.^*a^z 


3/2  2  2 


ia.^*a^z 


3/2 


and 


^ P  (x)  dx-21n  {m^*a^z 


3/2  > 


A- 


1 


(A^*a^z 


3/2  J  2 


For  z/D>4.0: 


P(x). 


5/2  3  2 


5/2 


j* P  (x)  dx-21n  (m^+a^z 


Appendix  D  Analytical  Solution  for  Constant  Horizontal 
Flame  Spread  Velocity 

For  a  semi-infinite  solid  subject  to  time  vaiying  heat  flux,  q"  on  the  surface,  the 
surface  temperature  can  be  evaluated  as  (Carslaw  and  Jeager,  1959); 


=/- 


.  // 

q 


■dx 


]/nJcpci 

where  is  the  surface  and  Tq  the  initial  temperatures. 

If  the  time  for  a  point  on  the  surface  to  reach  pyrolysis  temperature  is  tp,  then 


(D-1) 


T  -r^ - - —  f  -  dx 


(D-2) 


where  Tp  is  obviously  the  pyrolysis  temperature. 

In  the  case  that  q"  consists  of  two  heat  fluxes:  1)  a  constant  external  heat  flux,  qg,^,, 
and  2)  a  gas-phase  conductive  heat  flux  with  exponential  profile:  q"onv=A  exp[-(x-Xp)/6], 
equation  (D-2)  becomes: 


T  -r„- 

P  0 


•1  ^ 

=f 


-dx 


(D-3) 


Assuming  constant  flame  spread  velocity,  V,  then  x=Vtp  and  Xp(t)=Vt.  Therefore  equation 
(D-3)  can  be  transformed  as 


ynjcpc  Q 

For  the  time-being  we  only  deal  with  the  exponential  heat  flux  or  the  first  integral.  We 
assume: 

tp-T=a) 
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therefore  -dT=d<o 

for  T  from  0  to  tp,  co  from  tp  to  0,  we  obtain: 


P 

/ 


Ae 


-dx 


p 

■/ 


■da 


(assuming  ci)=(t)^,  ^  from  0  to  tp‘^^) 


where  P=(V/6)*^^cf) 

Equation  (D-4)  can  therefore  be  simplified  (for  constant  external  heat  flux,  qe,^): 


T  -r„- 

P  0 


2A 


^/nkpc\ 


y'nArpC 


using  Vtp=x: 


T  -r„-2 

P  0 


^n/cpc 


^nArpc  ^ 


(D-5) 


The  integral  is  exactly  an  error  function  erf[(x/6)‘'^]  (including  2ln''^). 

If  we  X  is  large  relative  to  6  (this  is  usually  case  since  6  is  small  (order  of  1  mm)),  then 
(x/6)*^  ->«>  and  erf  («)=!.  Therefore  equation  (D-5)  is  changed  to: 


T  -T  • 

p  0 


A^/yfV  2 


\/kpc  ^Jnkpc  “‘'N 


X 

V 


(D-6) 


where  the  second  term  of  the  RHS  is  essentially  the  surface  temperature  owing  to  external 
heat  flux. 

Therefore  equation  (D-6)  can  be  written  as: 
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y  y  A^/yfV 


(D-7) 


We  can  notice  that  time  to  reach  pyrolysis  from  the  surface  temperature,  Tj  is: 

n  kQC{T^-Ty 

^ps'  4 


Combining  this  with  equation  (7),  one  has: 

dt  4  t 

ps 


(D-8) 


Equation  (D-8)  is  the  analytical  asynqttotic  solution  for  constant  flame  spread  velocity. 
On  the  other  hand,  if  one  can  obtain  constant  flame  spread  velocity,  V,  the  combination  of  A 
and  6  (A6  can  be  deduced  (cf.  equation  (D-7)). 
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Appendix  £  Justification  for  Gas-Phase  Conductive 
Heat  Flux  Profile  Being  the  Same  as 
Surface  Temperature  History 


The  sur&ce  temperature  rise  from  Carslaw's  eqn.: 


T,(t)To^ 


yl^pc  {  -Jt-T 


dr+T^-To 


(E-1) 


where  4'  ('t)  is  the  flame  heat  flux,  assuming  no  surfrce  reradiation  loss.  Therefore: 


r(r) 


4o"e 


x.j^ 

s  ■■ 


Substituting  this  into  eqn.  (E-1)  and  ignoring  4'  for  a  moment,  one  has; 

O 


Ts(t)-To  = 


lees 
•Jjikpc  \  t-T 


dr+T^-To 


(E-2) 


Ignoring  Text-To  on  the  right  hand  side  for  moment  without  losing  generality: 
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assuming; 


hence: 


6 


dr  =  -2G)S /vdct) 


therefore: 
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b/v  2co  da-ly/bjv  J  e 


dco 


substituting  this  into  eqn.  (E-3),  the  result  is; 

Tjt)T^ - - - e  5  2  f  e-“'dC0 

i/nkpc  i 


(E-4) 
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2  2 

—  f  e^dcO-1 

i 


Then  putting  back  all  the  variables  into  eqn.  (E-4),  finally  it  becomes: 


TAt)-T 
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giyt) 


^nkpc 


yfbjv 


(E-5) 


One  can  see  from  above  equation  that  the  profile  of  surface  temperature  is  exactly  the  same 
as  that  of  the  gas  phase  conductive  heat  flux. 
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Appendix  F  Algorithm  of  Heat-up  Calculation  for  a 
Sudden  Jump  (Drop)  in  Heat  Fluxes  at 
the  Surface 

When  the  &el  surface  experiences  sudden  jump  (drop)  of  the  heat  fluxes  (e.g. 
sharp  increase  of  flame  heat  flux  when  flame  approaches;  sharp  decrease  of  total  heat  flux 
when  water  applies  to  the  fuel  surflice),  the  integral  heat-up  model  has  problem  in 
calculating  sur&ce  temperatures  (Delichatsios,  M.M.  et  al.,  1990).  Therefore,  the  results 
of  flame  spread  calculation  are  not  reliable.  This  error  arises  fi'om  the  flict  that  there 
exists,  in  essence,  two  temperature  profiles  within  the  solid:  one  being  profile  due  to 
previous  heat  flux  without  sudden  jump  (drop)  and  the  other  being  profile  caused  by  the 
additional  sudden  jump  (drop)  heat  flux  (es).  The  following  description  will  show  the 
reader  that  one  can  separate  the  effects  of  two  heat  fluxes  and  superimpose  the  two 
surface  temperatures  calculated  independently  fi'om  these  two  heat  fluxes. 

We  begin  with  one  dimensional  heat  conduction  equation,  boundary  and  initial 
conditions: 


dt  ^ 

(F-1) 

_  de. 

II 

l-e 

1 

(F-2) 

II 

1 

1 

(F-3) 

0(2,0)  =0 

(F-4) 

where  0=T-To  and  To  is  the  initial  temperature  of  the  solid;  q\  is  the  total  net  heat  flux  at 
the  fuel  surface  including  external,  flame  as  well  as  sur&ce  reradiation  loss.  We  will 
ignore  surfikce  reradiation  loss  for  simplidty  of  current  discussion. 

When  a  sudden  jump  (drop)  of  heat  flux  is  {q)plied  to  the  fuel  surface  at  a  time,  ti. 

The  total  iiet  heat  flux  at  the  fuel  surface  is  qi ,  the  net  increase  of  the  heat  flux  is  qi-qx. 
One  can  assume  there  are  two  temperature  rises  (drops)  (0i  and  02)  associated  with  two 

heat  fluxes  {qi  and  q2-q\ .)  and  set  timilar  equations  as  Equation  (F-1  to  F-3): 
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(F-Sa) 
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dz 


l*=0 


ddj  ^02 
dt  ^  ^z 


-  k 


dz 


\2=0  ^2  ^2 


0 


(F-5b) 

(F-5c) 

(F-6a) 

(F-6b) 

(F.6c) 


Combining  Eqns.  (F-5a)-(F-5c)  withEqn.  (F-6a)-(F-6c)  one  can  obtain; 

(F-7a) 

(F-7b) 

(F-7c) 

where  0=6i+02 .  It  can  be  seen  that  equation  (F-7a)-(F-7c)  are  exactly  same  exactly  same 
except  q  2  replaces  q  1 ,  which  is  expected  because  qt  is  the  n^  heat  flux  after  time,  ti. 

F-2 


dt  ^  &z 


^  dz 


One  can  solve  equations  (F-5a)-(F-5c)  and  equations  (F-6a)-(F-6c)  using  the  same  way 
that  solves  equations  (F-l)-(F-4)  with  the  initial  conditions  changes  to; 


-  0i(x,11)  (F-7d) 

and 

02(x,tl)  =  O  (F-7e) 

In  the  case  that  reradiation  loss  is  taken  into  consideration,  we  can  change  Eqn.  (F- 
5b)  and  Eqn.  (F-6b)  to: 

80 

-k-^l^  =  q,e<TTt,  (F-8a) 

80 

-k-^U  =  q,qjea(TtTti)  (F-8b) 

The  above  equation  can  be  easily  implemented  into  the  integral  preheat  model  to 
calculate  flame  spread. 
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Appendix  G  The  input  of  the  Horizontal  Flame  Spread 
Code 

There  are  numerous  variables  that  are  used  as  input  to  the  model.  In  general,  the 
external  heat  flux,  flame  gas-phase  conductive  heat  flux  and  flame  radiation  (at  least  at  this 
moment)  have  to  be  prescribed.  The  material  properties  have  to  be  known.  The  input 
parameters  along  with  their  explanations  listed  below  can  help  interested  users  to  run  the 
model.  Two  input  files  are  needed  to  run  the  flame  spread  code  hfsspv70.f.  They  are 
named  hflinp  and  fluxinp.  The  first  input  file,  hflinp,  contains  all  the  necessary  input  and 
the  second  file,  fluxinp,  is  used  to  specify  the  external  flux. 

hflinp  input  file  listing 

19.93, 1000, 1 .27, 1 92,371,592, 

2.20.. 33.5, 

250. 1200.. 0.,  1 100.,2.2, 1 .,  1 .,  1 .06E-4,2.7984E-4, 1 .,  1000.,/N,Rhov,RhoC,Htpyrl,Cpv,Cpc, 

Cpg,tkc,tkv,dt,Timend 

0.1,1.,  /E1,XN 

.30,l.,23400,8.5,296.15,643.15,5.67e-ll,.01/CHIR,ChiA,Htcomb,sc,Tinf,Tpyr,epsig,T% 

/(x(i),i=l,n) 

30.. 30..0.0.0,  ,.0,00,.00,.00,.000,.000,.000,.000,.000,.000/(QE(i),i=l,n) 


Nomenclature 
Line  1 

Line  1  parameters  are  only  used  to  predict  Horizontal  Flame  Spread  experiment,  users  can 
ignore  this  line  (but  keep  it  here) 

Q42;  The  nominal  heat  flux  at  the  425  mm  location  from  the  lower  edge  of  the  radiant 
panel,  kW/m^2 

TMV;  Time  to  start  moving  the  slide/sample,  seconds 
SMV:  The  slide/sample  moving  speed,  mm/s 

Nl,  N2,  N3:  the  surface  temperature  output  at  desired  points  which  corresponds 
to  the  experimental  conditions 


Line  2 

kchar:  The  flag  for  indicating  if  the  fuel  is  charring  or  non-charring  material,  if  charring, 
kchar=l;  noncharring,  kchar=2 

crm:  Critical  mass  flux  that  can  sustain  burning  (do  not  use  currently) 

Qig:  The  ignition  heat  used  for  first  two  nodes  (kW/m'^2) 
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Line  3 

N:  Number  of  Nodes  used  for  the  fuel  surface  in  x  direction 

Rhov:  The  density  of  the  virgin  fuel,  kg/cu  m 

Rhoc:  The  density  of  the  charred  fuel,  kg/cu  m 

Htpyrl:  The  heat  of  pyrolysis,  kJ/kg 

Cpv:  The  specific  heat  of  the  virgin  fuel,  kJ/kg.K 

Cpc:  The  specific  heat  of  the  charred  fuel,  kJ/kg.K 

Cpg:  The  specific  of  the  ambient  gas,  kJ/kg.K 

tkc:  The  thermal  conductivity  of  the  charred  fuel,  kW/m,K 

tkv;  The  thermal  conductivity  of  the  virgin  fuel,  kW/m.K 

dt:  The  initial  time  step  used  for  the  calculation,  seconds.  If  the  code  is  converged 

within  dt,  dt  is  also  the  output  interval 
Timend:The  final  time  of  the  calculation,  seconds 

Line  4 

El:  The  fuel  thickness,  m 

XN:  The  fuel  length,  m  (not  used  currently) 

Xl^owr:  power  exponent  used  for  Xf  computation,  if  no  input,  Xfpowr=0.6667 
Xfcof:  The  coefficient  used  for  Xf  computation 

Line  5 

CHAIR:  The  radiation  fraction  of  the  flame 
Chi:  The  combustion  efficiency  of  the  flame 
Htcomb:The  heat  of  combustion  of  the  fuel,  kJ/kg 
sc:  The  Mass  air/fuel  stoichiometric  ratio 

Tinf:  The  initial  temperature  of  the  fuel,  K 
Tpyl:  The  p3Tolysis  temperature  of  the  fuel,  K 

epsig:  The  fuel  surface  emissivity  x  Stafan-Boltzmann  constant,  kW/sq  m.  KM 
T%:  A  percentage  of  Tp:  (Tp-T)/Tp  used  for  calculation,  normally  0.01 

Line  6 

x(i):  Nodal  specification,  e.g.  node  1,  0,  node  2,  0.1,  m,  if  no  input  in  this  file,  the 
program  calculates  x(i)  using  N 

Line  7 

Qe(i):  External  heat  fluxes  for  the  nodes,  kW/sq  m 

Comments 

When  a  non-charring  fuel  is  used  in  the  flame  spread  calculation  (i.e.  kchar=2),  the 
density  of  char  (rhoc)  is  set  rhoc=0  inside  the  program. 

When  a  charring  fuel  is  used  (kchar=l),  the  code  sets  rhoc=0  only  in  the  pyrolysis 
subroutine  to  speed  up  convergence.  Justification  for  using  rhoc=0  is  presented  in  Chen's 
thesis  (Chen,  1991). 
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fluxinp  input  file  listing 


3,0.,20.,0.3/KEFLUX,QCEXT,QE0,XEXT 

10.,0.1,/QFRO,DFRAD 

70.0,0.0015,/QFC0,DFC 

Line  1 

KEFLUX:  external  heat  flux  option;  1 -experiment,  2-constant,  3 -exponential 
QCEXT:  the  value  of  constant  external  heat  flux,  kW/m^2  ifKEFLUX=2 
QEO;  the  magnitude  of  exponential  external  heat  flux,  kW/m^2  if  KEFLUX=3 
XEXT;  the  length  scale  of  exponential  external  heat  flux,  m  ifKEFLUX=3 

LINE  2 

QFRO;  the  magnitude  of  exponential  flame  radiative  heat  flux,  kW/m'^2 
DFRAD:  the  length  scale  of  exponential  flame  radiative  heat  flux,m 

LINE  3 

QFCO:  the  magnitude  of  exponential  flame  condcutive  heat  flux,  kW/m'^2 
DFC:  the  length  scale  of  exponential  flame  conductive  heat  flux,m 
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Appendix  H  Experimental  Data  and  Comparisons  with 

Numerical  Predictions 

In  this  section,  processed  experimental  data  and  comparisons  with  the  numerical 
predictions  of  the  horizontal  flame  spread  model  are  listed.  The  first  page  shows  the 
comparison  of  readings  of  large  and  small  heat  flux  gages  during  a  inert  test  where  86862  is 
a  small  flux  gages  and  others  are  large  flux  gages  (two).  The  rest  of  the  materials  are 
arranged  in  the  following  order:  1)  A  sheet  recording  experimental  conditions,  2)  the 
temperature  history  measurements,  3)  heat  flux  history  measurement  from  the  two  heat  flux 
gages,  4)  flame  heat  fluxes  as  function  of  the  distance  from  the  cooling  plate  edge,  5)  detailed 
surface  temperature  plots  for  first  two  stations  (if  there  is  any),  6)  the  comparison  of 
experimental  data  of  surface  temperature  at  the  station  3  with  the  numerical  prediction  using 
heat  fluxes  measured  by  the  large  flux  gage  (not  for  every  case),  and  7)  the  comparison  of 
experimental  data  of  detailed  surface  temperatures  at  the  stations  1  and  2  (if  there  is  any)  with 
the  numerical  prediction  (not  for  every  case). 
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TP — The  difference  between  pyrolysis  temperature  and  the 

ibient  temperature 
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ENDIF 

CALL  ENTHA2(COSM,VLAMDA,TCE) 
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IF(VLAMDA.LT.O.)  VLAMDA=0. 
VN=VLAMDA*COSM 
VM= VLAMDA*(1  .-COSM) 
GLAMDA=GAMMLN(VLAMDA) 
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VLAMDA=8,09 

ENDIF 

CALL  ENTHA1(C0SM,VLAMDA,TC) 
VOLK = (TC-TA)/(TA*(SIT-SCC)) 
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9/3/93 

IF(COSM.LE.0.068)  THEN 
VLAMDA=(l.-COSM)/(1.5*COSM)-l. 
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ENDIF 

CALL  ENTHA1(C0SM,VLAMDA,TC) 
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THE  FLAME  ENTHALPY  AT  HIGHT  Z,  MIXTURE  FRACTION  COSM 
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FUNCTION  VOLI3(SCC) 
IMPLICIT  REAL*8(A-H,0-Z) 
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Now  the  unit  of  CPG(I)  is  in  KJ/Kmol.K.  This  should  be  changed 
as  KJ/kg.K 
CPG(l)=CPG(l)/28. 
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EXTERNAL  SPEC1,SPEC2 
EXTERNAL  MIDSQUl 
EXTERNAL  MroSQU2 
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FUNCTION  SPECl(CO) 

PARAMETER  (NPAR=64) 

IMPLICIT  REAL*8(A-H,0-Z) 

REAL*8  YS(NPAR),DELH(NPAR),A(NPAR,5),B(NPAR,5) 
COMMON/BASIS/A, B,DELH 
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IF(POl.GE.0..AND.POl.NE.l.)  THEN 
BB1=0.5**VN 
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IF  (C0.LE.EXP((l.-P01)*(-86.)))  THEN 
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COMMON/COMP/YS 
common/pw/pol  ,P02 
IF  (PO2.GE.0.)  THEN 
COl  =  1  .-CO*t'(l  ./(I  .-P02)) 
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IF(P01.LT.O.)  THEN 
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ELSE  IF  (POl.GE.O)  THEN 
SPEC6=TEMP*DIV*(1.-C01)**(-P02)/(1.-P01) 
ELSE 
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ELSE  IF  (PO2.GE.0)  THEN 
IF(.P01*DL0G(C01).LE.(-86.))  THEN 
SPEC7=0. 
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TEMPI = ZBRENT2(CALTE,T0,TMAX,EPS) 
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Y(I)=0. 

ENDIF 

CONTINUE 
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P2=(-0.79481-3.33465*XM+7.74823*XM**2- 

8.58911*XM**3+4.31337*XM**4-0.810195*XM**5) 
IF  (P2.LE.(-43.))  THEN 
YM(2)=0. 
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comi=Y(2) 

DO  120  1=1, N 

IF  (Y(I),GT.  1.0)  THEN 

Y(I)=1. 
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SUBROUTINE  QGAUS3(FUNC,A,B,SS) 
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EXTERNAL  FUNC 
The  abscissa  and  weights. 

DATA  X/.  1488743389, ,4333953941, .6794095682, .8650633666, 
.9739065285/ 


DATA  W/.2955242247,. 2692667193, .2190863625, .1494513491, 
.0666713443/ 

XM=0.5*(B+A) 

XR=0.5*(B-A) 
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SS = SS + W(J)*(FUNC(XM + DX) + FUNC(XM-DX)) 
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by  ten-point  Gaxiss-Legendre  intergration:  the  function  is 
evaluated  exactly  t&a  times  at  interior  points  in  the  range  of 
integration 
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FUNCTION  ZBRENT1(FUNC,X1,X2,T0L) 
IMPLICIT  REAL*8(A-H,0-Z) 
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IF(ABS(XM).LE.TOLl.OR.FB.EQ.0.)THEN 
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Attempt  inverse  quadratic  interpolation 
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Appendix  K  Description  of  Execution  of  the  Program 

and  Output 

The  main  code  of  the  HFSCS  is  named  hfsspv70.f  This  program  is  written  in  FORTRAN 
77  and  consist  of  a  large  number  of  subroutines.  The  input  files  for  the  program  are  hflinp 
and  fluxinp  which  are  described  in  details  in  Appendix  G.  A  small  number  of  input  may  be 
required  to  change  the  parameters  within  the  code  and  if  the  output  needs  to  be  modified. 
The  nodal  spacing,  Ax ,  can  be  either  fixed  or  variable  and  is  set  in  the  code  or  specified  in 
the  iput  file  hflinp.  The  external  radiative  flux  distribution  is  specified  by  ihe  fluxinp  file 
where  3  choices,  experimental  (HTFT  distribution  profile),  constant  or  exponential  are 
available  as  described  in  Appendix  G.  The  flame  radiative  and  gas-phase  conductive  flux 
profiles  are  also  specified  using  peak  and  length  scale  values  for  exponentially  decaying 
profiles  by  fluxinp. 

There  are  two  output  files  currently  active  in  the  program,  although  a  number  of  other 
output  files  have  been  set-up  (and  are  documented  in  the  code)  for  debugging  and  special 
data  outputs  (e.g.  to  make  contour  plots  of  temperature)  are  commented  out  and  can  be 
readily  used.  The  primary  output  files  are: 

hfssprf  is  the  main  output  file  reporting  the  temperature  and  flux  distribution  on  the  fuel 
surface,  respectively,  at  specified  spatial  distribution  by  the  user.  The  time  interval  of  the 
output  is  also  user  specified.  These  specifications  are  made  in  the  code  at  the  output 
subroutine. 

hfssout.  is  another  output  file  containing  some  of  the  key  parameters  that  are  selected  by 
the  user  when  the  code  is  initially  executed.  These  are  the  burn-out  front,  Xb,  pyrolysis  or 
flame  fi'ont,  Xp,  flame  height,  Xf,  rate  of  mass  burning,  Emtotl,  the  chemical  heat  release 
rate,  Qch  and  the  incident  flux,  Qinc.  The  parameters  Xf  and  Qch  are  based  on  a 
calculated  mass  loss  rate,  Emtotl  and  are  calculated  using  an  empirical  relation  for  Xf  and 
the  heat  of  combustion  for  Qch.  However,  Emtotl  would  need  an  accurate  value  for  the 
radiation  feed-back  from  the  flame  to  the  burning  surface.  This  can  be  accomplished  using 
a  combustion  model  and  a  radiation  model  that  could  account  for  the  soot  concentration 
near  the  flame  base.  Since  such  models  are  still  under  development,  a  true  coupling 
between  the  combustion  model  and  the  flame  spread  can  not  be  accomplished  without 
further  work. 

The  combustion  model,  /7oo/2./ (Appendix  J),  relies  on  the  combb.dat  (Appendix  I)  file  for 
input  and  a  modified  program  raJM/ (listing  not  provided  but  included  on  a  diskette)  can 
be  used  for  jet  flames  with  an  input  file  combbl.dat  (Appendix  la). 


The  combustion  model  can  be  pseudo  coupled  to  the  EDFSCS  by  using  the  radiation 
vertical  profile  output  of  the  combustion  model  to  estimate  the  flame  radiative  profile. 
The  HFSCS  can  be  simply  modified  to  vary  the  peak  radiative  flux  as  the  burning  area 
grows.  As  noted  before,  lack  of  basic  data  for  state  relationships  for  real  fuels,  time 
intensive  calculation  of  the  combustion  model  and  difficulty  in  developing  a  robust  and 
accurate  radiation  model  prevents  a  true  coupling  of  the  HFSCS  with  the  combustion 
model. 
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imem 

ivieasuremenis  oi  lemperature  ana 
Fluxes,  Data  Sheets  and  External  Flux 
Profile.  Plots  are  results  after  significant 
processing  of  the  raw  data 


The  first  plot  shows  the  comparison  of  the  external  flux  profile  measured  by  three  different 
flux  gages  as  the  sample  holder  was  moved  at  a  speed  of  0.5  m/s.  Flux  gages  27S44  and 
525842  are  the  large  gages.  They  were  placed  at  different  locations,  but  the  data  was 
processed  to  show  the  plot  as  a  fiinction  of  distance  from  the  leading  edge.  These  two 
gages  measured  approximately  the  same  flux  level  (there  may  be  some  small  variations  in 
time  of  the  radiant  panel).  The  small  gage  (86862)  clearly  measures  a  higher  flux  level.  It 
is  noisier  and  does  not  drop  to  a  near  zero  when  reaches  the  leading  edge  of  the  cooling 
plate  (as  other  two  gages  do).  This  indicated  that  the  small  gage  was  not  reliable.  We 
also  used  a  large  gage  (1  inch  in  diameter  and  the  same  type  as  27844  and  525842) 
calibrated  at  NIST  to  check  the  measurements. 

Arrangement  of  the  information 

For  each  test,  the  following  is  contained  in  this  appendbc; 

1 .  Data  sheet,  listing  all  the  pertinent  information 

2.  Plot  of  surface  and  gas  temperatures  measured  vs.  Time,  noting  the  time  of  slide 
movement  and  speed  of  movement. 

3 .  Total  heat  flux  measured  at  the  surface  as  a  fimction  of  time. 

4.  Kame  heat  flux  (total  -  externally  imposed  radiation)  profile  ahead  of  the  flame.  It  is 
obtained  by  transforming  the  flux  vs.  time  data  using  the  traversing  speed  of  the  slide. 

5.  Prediction  vs.  measurement  of  the  surface  temperature  during  heat-up  using  the  flame 
flux  measurement  (which  would  be  aU  flame  radiation  at  some  finite  distance  before 
the  flame  leading  edge  region).  This  procedure  was  used  to  determine  the  material 
thermal  inertia,  kpc. 

6.  Detailed  plots  of  the  surface  temperature  measurements  near  the  flame  fi'ont  to 
estimate  the  pyrolysis  temperature  and  perform  the  analysis  described  in  6.5  and 
.^pendixE. 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 


Experimental  No.:  21  Date:  6/17/94  Tu  Time:10:30pm 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  10 
Slider  Speed  (mm/s):  0.32  mm/s 

Burning  Sample  Data 

Material:  3/4"  Potlatch  Particle  Board  Dimension(LxWxH,  mm):  600x152x18.85 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

File  Column 

Note 

1 

TC 

175 

0.79 

B 

GAS  TCI 
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SURFACE 
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3 

TC 
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0 
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3 

HG 
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0 

i 

H 

SER.  525842 

Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  3.70-12.7576  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  270 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  20 

Speed  Used  (SlM"Steps"):  50 

Number  of  Steps  (IlM'’Steps"):  80000 

File  Names  (.PRN):  T10K6i7 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  402.2 
Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (°C):  23 
Observations  Burning  area  ~  10  cm 
Personnel  V.M.,  Y.C 

Processing  of  Data:  Tp=320X,  Surface  TCI,  1/t=0.14,T^=240°C,  kpc=0.17x765.x2800=364140, 
q'^==(kpc)''^(Tp-T«0(l/T)*'H364140)*'2x(320-240)(0,14)*'^18.06kW/ra'' 

Surface  TC2,  1/t=O.10,  T«,=240°C,  q’<;=(364140)‘'^x(320-240)(0.10)''^15.27  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 


Experimental  No.:  22  Date:  6/20/94  Tu  Time:8:30pm 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  10 
Slider  Speed  (mm/s):  0.32  mm/s 

Burning  Sample  Data 

Material:  3/4"  Potlatch  Particle  Board  Dimension(LxWxH,  mm):  600x152x18.85 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

t 

X  (mm) 

z  (mm) 

FUe  Column 

Note 

1 

TC 

170 

0.8 

B 

GAS  TCI 

1 

TC 
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0 

C 
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TC 
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0.41 

D 
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HG 
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TC 
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0 

F 
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3 

HG 
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0 

H 

SER.  525842 

Experiment 

Hux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  3.70-12.7576  kWtn? 

Preheating  Time  without  PUot  Flame  (seconds):  240 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  20 

Speed  Used  (SlM'Steps"):  50 

Number  of  Steps  (IlM"Steps"):  80000 

File  Names  (.PRN):  T10K620 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  die  Sample  (sec):  309.8 
Ignition  Heat  Flux  (Calculated,  kW/m*):  16.7 
Ambient  Conditions 
Temperature  ("C):  24 

Observations  Burning  area:  -10  cm,  -2  cm  burnout  area  so  the  burning  width  is  in  fact  8  cm. 
Personnel  J.C,  Y.C 

Processing  of  Data:  Tp=320’’C,  Surface  TCI,  l/t=0.20,T„=238.7°C,  kpc=0.17  x765.x2800=364140, 
qS  =(kpc)''^(Tp-T„^(l/T)‘"=(364140)''^x(320-238.7)(0.20)>'^=21.94  kW/m^ 

Surface  TC2,  1/t=0.078,T„,=238.7’’C,  qS=(364140)''^x(320-238.7)(0.078)‘'^=13.70  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 

Experimental  No.:  23  Date:  6/21/94  Wed  Time:  10:30  am 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  10 
Slider  Speed  (mm/s):  0.32  mm/s 

Burning  Sample  Data 

Material:  3/4"  Potlatch  Particle  Board  Dimension(LxWxH,  mm):  600x152x18.85 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

z  (mm) 

File  Column 

Note 

1 

TC 

167 

0.7 
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GAS  TCI 
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Experiment 


Rux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  3.70-12.61  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  252 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  18 

Speed  Used  (SlM’Steps’):  50 

Number  of  Steps  (IlM’Steps"):  80000 

File  Names  (.PRN):  T10K621 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  355.8 
Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  ("C):  24 

Observations  Burning  area:  14-4=10  cm,  -4  cm  burnout  area,  HGl  did  not  flush  with  the  surface  and 
SUR  TC3  had  some  problems  in  the  middle  of  the  test 
Personnel  J.C,  Y.C 

Processing  of  Data:  Tp=320‘>C,  SurfeceTCl,  1/t  =0.088,T„,=225.8®C,  kpc=0.17  x765.  x2800=364140, 
qS=(kpc)*“(Tp-T^(l/t)‘"=(364140)'®x(320-225.8)(0.088)'"=16.86kW/m^ 

Surface  TC2,  l/t=0.071,T„,=225.8‘’C,  qS=(364140)''2x(320-225.8)(0.071)’'^=15.15  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 

Experimental  No.:  24  Date:  6/21/94  Wed  Time:  8:30  pm 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  20 
Slider  Speed  (mm/s):  1.28  mm/s 

Burning  Sample  Data 

Material:  3/4"  Potlatch  Particle  Board  Dimension(LxWxH,  mm):  600x152x18.85 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

z  (mm) 

File  Column 

Note 

1 

TC 

174 

1 

B 

GAS  TCI 
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SURFACE 
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Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  7.45-25.38  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  120 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  0 

Speed  Used  (SlM'Steps*):  200 

Number  of  Steps  (IlM'Steps"):  80(X)0 

File  Names  (.PRN):  T20K621 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  126.8 
Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (®C):  24 

Observations  Burning  width:  12  cm,  no  burnout  area,SUR  TC3  is  not  good 
Personnel  J.C,  Y.C 

Processing  of  Data:  Tp=320'’C,  Surface  TCI,  lh=0.6l,T^=2n°C,  kpc=0,17  x765.x2800=364140, 
q;=(kpc)*'2(Tp-T„^(l/t)*'2=(364140)*'^x(320-271)(0.61‘°=23.09kW/m2 

Surface  TC2,  1/t=0.077,T„,=271°C,  qS=(364140)‘'^x(320-271)(0.77)‘'^=30.16  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 

Experimental  No.:  25  Date;  6/22/94  Wed  Time:8:30pm 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^);  20 
Slider  Speed  (mm/s):  1.28  mm/s 

Burning  Sample  Data 

Material:  3/4"  Potlatch  Particle  Board  Dimension(LxWxH,  mm):  600x152x18.85 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

imQ^yiiiii 

z  (mm) 

File  Colunm 

Note 

1 

TC 

176 
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GAS  TCI 
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SER.  525842 

Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  6.00  mV-  20.44  kW/m^ 

Preheating  Time  without  Pilot  Hame  (seconds):  120 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  0 

Speed  Used  (SlM"Steps’):  200 

Number  of  Steps  (IlM"Steps");  80000 

FUe  Names  (.PRN):  T20K622 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec);  2500 
Time  to  Start  Moving  the  Sample  (sec):  129.2 
Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (°C):  23 
Observations  Burning  width  -7  cm 
Personnel  J.C,  Y.C 

Processing  of  Data:  Tp=320®C,  Surface  TCI,  1/t=0.77,T„=259.4«C,  kpc=0.17  x765.x2800=364140, 
q;=(kpc)‘'2(Tp-T^(l/i:)"^=(364140)‘'2x(320-259.4)(0.77)'"=32.09kW/m^ 

Surface  TC2,  l/t=0.68,  T„,=259.4'’C,  qj=(364140)‘'^x(320-259.4)(0.68)''^=30.16  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 


Experimental  No.:  41  Date:  8/16/94  Tu  Time:4.00pm 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  20 
Slider  Speed  (mm/s):  1.28  mm/s 

Burning  Sample  Data 

Material:  3/4"  Potlatch  Particle  Board  Dimension(LxWxH,  mm):  600x152x18.85 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

z  (mm) 

File  Column 

Note 

1 

TC 

175 

1.47 

B 

GAS  TCI 

1 

TC 

175 

0 

C 

SURFACE 

2 

TC 

350 

1.09 

D 

GAS  TC2 

2 

TC 

350 

0 

E 

SURFACE 

2 

HG 

350 

0 

G 

SER.  86862 

3 

TC 

573 

i 

0 

F 

SURFACE 

3 

HG 

573 

0 

H 

SER.  525842 

Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  5.80-19.9984  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  120 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  7 

Speed  Used  (SlM'Steps*):  200 

Number  of  Steps  aiM"Steps'):  80629 

File  Names  (.PRN):  T20K816 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  147.0 
Ignition  Heat  Flux  (Calculated,  kW/m^):  29.4 
Ambient  Conditions 
Temperature  (°C):  22 
Observations  Burning  area  -10  cm 
Personnel  V.M.,  Y.C 

Processing  of  Data:  Tp=320‘>C,  Surface  TCI,  l/t=0.57,T„=274.7‘’C,  kpc=0.17  x765.x2800=364140, 
q;=(kpc)‘'^(Tp-T^(l/T)''*=(364140)‘'^x(320-274.7)(0.57)‘"=20.64kW/m^ 

Surface  TC2,  l/t=0.79,  T«=274.7°C,  <5;=(364140)"^x(320-274.7)(0.79)*'^=24.30  kW/m^ 


T20l<81  6 

Potlatch  Particle  Board,  V=l  ,28  rnm 


^  CM  rO 

O  O  CJ  O  O 


<  ~>  <  Z)  Z) 
O  (/)  O  oO  00 


(3o)  3dniVd3dl/^31 


o 

o 

o 

O 

o 

O 

O 

O 

o 

o 

o 

00 

CD 

CnJ 

O 

00 

CD 

CM 

o 

n 

rO 

rO 

rO 

ro 

Cs) 

CNJ 

OnI 

CM 

C\l 

(Oo)  3dnivd3diAi3i  dovjdns 


200 


o 

O 

o 

O 

o 

O 

o 

ro 

(Oo)  3dniVEI3dl^31  30VJdnS 

200  300  400 

TIME  (sec) 


(Oo)  3ymVd3dl/^3130VJdriS 


Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  43  Date:  8/18/94  Th.  Time:9:30pm 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  20 
Slider  Speed  (mm/s):  1.28  mm/s 

Burning  Sample  Data 

Material:  3/4"  Potlatch  Particle  Board  Dimension(LxWxH,  mm):  600x152x18.85 


Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

z  (mm) 

File  Column 

Note 

1 

Tr* 

174 

1  47 

B 

1 

•vn 

174 

0 

r 

O 

TT* 

1  27 

D 

HMMHMIi 

o 

•m 

'x^'y 

0 

E 

xjrz 

'X^O 

0 

G 

n 

T-r* 

0 

F 

3 

HG 

575 

0 

H 

SERIAL  525842 

Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  5.90  mV-  20.3432  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  120 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  0 

Speed  Used  (SlM"Steps"):  200 

Number  of  Steps  (IlM"Steps"):  80629 

File  Names  (.PRN):  T20K818 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  Ae  Sample  (sec):  141.2 
Ignition  Heat  Flux  (Calculated,  kW/m^):  23.4 
Ambient  Conditions 
Temperature  (®C):  22 

Observations  Burning  width  ~7.5  cm,  constant 
Personnel  J.C,  Y.C 

Processing  of  Data:  Tp=320^C,  Surface  TCI,  1/t=0.44,T^=280.6“C, 

kpc=0.17x765.x2800=364140,  (5S=(kpc)‘'^(Tp-Tea)(l/x)‘"=(364140)*"x(320-280.Q(0.44)‘'^=15.75 
kW/m^ 

SurfeceTC2,  1/t=0.30,  Te«=280.6“C,  (i;=(364140)‘'"x(320-280.6)(0.30)‘'^=13.02  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 


Experimental  No. :  27  Date:  6/23/94  Th.  Time:8:30pm 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  30 
Slider  Speed  (mm/s):  2.56  mm/s 

Burning  Sample  Data 

Material:  3/4"  Potlatch  Particle  Board  Dimension(LxWxH,  mm):  600x152x18.85 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station  Sensor  x  (mm)  z*  (mm)  File  Colunm  Note 


3  HG  575  0  H  SERIAL  525842 


Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  9.65  mV-  32.88  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  50 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  0 

Speed  Used  (SlM"Steps"):  400 

Number  of  Steps  (IlM"Steps"):  80000 

FUe  Names  (.PRN):  T30K623 

Set-up  File  Name; 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  55 
Ignition  Heat  Flux  (Calculated,  kW/m^):  40 
Ambient  Conditions 
Temperature  (°C):  23 

Observations  Burning  width  -10  cm,  flame  size  little  bit  large 
Personnel  J.C,  Y.C 

Processing  of  Data:  Tp=320°C,  Surface  TCI,  l/T=4.0,Te„=289®C,  kpc=0.17  x765.x2800=364140, 
qS  =(kpc)*"(Tp-Te«)(l/T)‘'2=(364140)‘'2x(320-289)(4.0)*^=37.71  kW/m^ 

Surface  TC2,  1/t  =  1.29,  Te«=251.5‘'C,  qS=(364140)’'2x(320-251.5)(1.29)‘'2=46.95  kW/m" 
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Potlatch  Partical  Board,  V=2.56  mm/s 
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T30K624 

Potlatch  Partical  Board,  V=2.56  mm/ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.;  2  Date:  4/18/94  Mon.  Time:  10:30  pm 

Pre-Experiment 

Heat  Flux  at  400  mm  Location  (Planned)  (kW/m^):  10 
Slider  Speed  (mm/s):  0.64  mm/s 

Burning  Sample  Data 

Material:  1/4"  PMMA  Sheet  Dimension(LxWxH,  mm);  800x152x6.35 

Instrumentation 

Number  of  Thermocouples:  6  Number  of  Heat  Flux  Gauges:  2 


Station  Sensor  x  (mm)  z  (mm)  File  Column  Note 


3  HG  599  0  I  SERIAL  525842 


Experiment 

Flux  Gauge  Reading  @  400  mm  (Serial  No.  27844,  mV):  3.10-10.56  kW/m^-10.16  kW/m^  @425  mm 

Preheating  Time  without  Pilot  Flame  (seconds):  150 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  20 

Speed  Used  (SlM"Steps"):  100 

Number  of  Steps  (IlM"Steps"):  60000 

File  Names  (.PRN):  P10K418a,  P10K418b 

Set-up  File  Name: 

Sampling  Rate  (Hz):  10  Duration  of  Sampling  (sec);  600 
Time  to  Start  Moving  the  Sample  (sec):  -170 
Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (°C):  18. 

Observations  Burning  width  remains  at  about  6  cm,  slide  stops  @60000  steps  (travels  over  381  mm) 
then  stops  the  data  acquisition  for  10,  seconds,  begin  data  file  P10K418b  which  is  the  data 
file  when  the  slide  stops 
Personnel  Jeff  Cote,  Y.C. 

Processing  of  Data:  SUR  TCI:  1/t=0.41,  Tp=370“C,  T„,=235.4°C, 

kpc=0.2798  X 1200  X2200=738672, 

qS = (7tkpc)i"(Tp-Te,J(V/6)''2=(738672)*'2  x  (370.-235.4)(0.41)''2=73 .91kW/m2 
SUR  TC2;  1/t=0.53,  Tp=370“C,  T^=235.4°C, 

q;;=(7tkpc)'"(Tp-TeJ(V/6)^'2=(738672)i/2x(370-235.4)(0.53)i^=84.48kW/m2 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No. :  11  Date:  5/16/94  Mon.  Time:2:30pm 

Pre-Experiment 

Heat  Flux  at  400  mm  Location  (Planned)  (kW/m^):  10 
Slider  Speed  (mm/s):  0.50  mm/s 

Burning  Sample  Data 

Material:  1/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  800x152x6.35 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


t 


Station  Sensor  x  (mm)  z  (mm)  File  Column  Note 


3  HG  602  0  H  SERIAL  525842 


Experiment 

Flux  Gauge  Reading  @  400  mm  (Serial  No.  27844,  mV):  3.20-10.90  kW/m^- 10.58  kW/m^  @425  mm 

Preheating  Time  without  Pilot  Flame  (seconds):  240 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  80 

Speed  Used  (SlM"Steps"):  79 

Number  of  Steps  (IlM"Steps"):  80000 

File  Names  (.PRN):  P10K516 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  1500 
Time  to  Start  Moving  the  Sample  (sec):  -376 
Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (°C):  16. 

Observations  Surface  &TC,  HG  aU  are  unpainted.  Area  is  constant  (10  cm  wide).  84501  water  outlet 
temperature:  16°C,  about  constant. 

Personnel  Jeff  Cote,  Y.C. 

Processing  of  Data:  SUR  TCI;  1/t=0.444,  Tp=370°C,  Te«=221.7°C, 

kpc=0.2798  X 1200  X2200  =  738672, 

q;;=(7tkpc)‘'^(Tp-Te«)(V/6)*'2=(738672)^"x(370.-221.7)(0.444)*'2=84.95kW/m2 
SURTC2:  1/t=0.333,  Tp=370°C,  T,«=221.7°C, 

q;;  =(iikpc)i^2(Tp.T^j(V/6)i^2=(738672)i^2 x(370.235.4)(0.333)1/2=73.57  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.;  5  Date:  4/29/94  Fr.  Time:  10:30  pm 

Pre-Experiment 

Heat  Hux  at  400  mm  Location  (Planned)  (kW/m^):  20 
Slider  Speed  (mm/s):  1.27  mm/s 

Burning  Sample  Data 

Material:  1/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  800x152x6.35 

Instrumentation 

Number  of  Thermocouples;  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

z  (mm) 

File  Column 

Note 

1 

TP 

1Q5 

1  57 

R 

RPISfSPHI^I 

_ 1 

TC 

.  195 

0 

r 

IlfflnfllBSBRnillii 

7 

TC 

377 

1  77 _ 

n 

_ 2 

TC 

377 

iRm^n 

F 

_ 2 

HG 

L377 

_ 0 _ 

a 

_ 3 

TC 

598 

IIRHHH 

F 

3 

HG 

598 

0 

H 

SERIAL  525842 

Experiment 

Flux  Gauge  Reading  @  400  mm  (Serial  No.  27844,  mV):  5.85-19.93  kW/m2^19.35  kW/m^  @425  mm 

Preheating  Time  without  Pilot  Flame  (seconds):  90 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  30 

Speed  Used  (SlM"Steps"):  200 

Number  of  Steps  (IlM"Steps"):  80000 

FUe  Names  (.PRN);  P20K429 

Set-up  File  Name: 

Sampling  Rate  (Hz);  5  Duration  of  Sampling  (sec):  1500 
Time  to  Start  Moving  the  Sample  (sec): 

Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (®C):  24.5 

Observations  Burning  area  is  good  (-6  cm).  Data  acqmdtion  starts  @35  seconds  after  starting  up  of 
preheating,  connection  to  TC  not  good.  Connection  to  HG2  discontinues  after  flame  firont  gets  to  second 
station 

Personnel  Y.  C. 

Processing  of  Data;  SUR  TCI:  kpc=0.2798x  1200  x2200=738672, 

SURTC2:  1/t=1.20,  Tp=370®C,  T^=245.7°C, 

qS = (itkpc)‘"(Tp-TeJ(V/6)‘'2=(738672)*'=‘ x (370-245. 7)(1 .20)^"= 1 17.0  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  6  Date:  5/04/94  Wed.  Time:  11 :30  am 

Pre-Experiment 

Heat  Flux  at  400  mm  Location  (Planned)  (kW/m^):  20 
Slider  Speed  (mm/s):  1.27  mm/s 


Burning  Sample  Data 

Material:  1/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  800x152x6.35 


Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

File  Colunm 

Note 

_ ] _ 

TC 

_ 201 _ 

0  77 _ 

r 

IPHISHHI 

_ 1 

TC 

R 

_ 2 

TC 

.  37Q 

F 

2 

TC 

379 

n 

_ 2 

!  HG 

_ 320 _ 

G 

_ 3 

TC 

F 

3 

HG 
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0 

H 

SERIAL  525842 

Experiment 

Flux  Gauge  Reading  @  400  mm  (Serial  No.  27844,  mV):  5.87-20.00  kW/m2-19.42  kW/m^  @425  mm 

Preheating  Time  without  Pilot  Flame  (seconds):  90 

Time  to  Sanq)le  Ignition  After  Applying  the  Pilot  Flame  (seconds):  20 

Speed  Used  (SlM"Steps"):  200 

Number  of  Steps  (IlM"Steps"):  80000 

File  Names  (.PRN):  P20K504 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  800 
Time  to  Start  Moving  the  Sample  (sec): 

Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (°C):  23. 

Observations  Surface  TC3  seems  fluctuating  at  the  end  of  the  test.  84501  reading  low.  High  heat 
spray  paint  was  applied  on  die  fuel  surface 
Personnel  Jeff  Cote,  Y.C. 

Processing  of  Data:  SUR  TCI:  1/t  =  1.20,  Tp=370°C,  Te«=260.4®C, 

kpc=0.2798  X 1200  X2200=738672, 

qo=(nkpc)^^(Tp-Text)(V/6)i'2=(738672)i'2x(37o..260.4)(1.2)‘^=103.4kW/m2 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  4  Date:  4/27/94  Wed.  Time:  11 :30  am 

Pre-Experiment 

Heat  Flux  at  400  mm  Location  (Planned)  (kW/m^):  20 
Slider  Speed  (mm/s):  1.27  mm/s 

Burning  Sample  Data 

Material:  1/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  800x152x6.35 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

t 

z  (mm) 

File  Colmnn 

Note 

_ 1 

TC 

0  85 

R 

_ 1 

TC 

HBIHiHii 

--  0 

r 

_ 2 

TC. 

378 

064 

n 

_ 2 

TC 

378 
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F 

_ 2 

hg 
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0 

G 

3 

TC 

.599 

_  0 

F 
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3 

hg 
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0 

h 

SERIAL  525842 

Experiment 

Hux  Gauge  Reading  @  400  mm  (Serial  No.  27844,  mV):  5.90-20.10  kW/m^- 19.68  kW/m^  @425  mm 

Preheating  Time  without  Pilot  Flame  (seconds):  90 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  60 

Speed  Used  (SlM"Steps''):  200 

Number  of  Steps  (IlM"Steps"):  80000 

File  Names  (.PRN):  P20K427 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  800 
Time  to  Start  Moving  the  Sample  (sec): 

Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (®C):  21. 

Observations  Stops  at  70481,  HGl  wire  disconnects,  burning  width  too  large  (8  cm),  constant  burning 
area  is  obtained 
Personnel  Jeff  Cote,  Y.C. 

Processing  of  Data:  SUR  TCI:  1/t=1.0,  T  =370‘’C,  Te«=232.4°C, 

kpc=0.2798  X 1200  X2200=738672, 

qS=(7ikpc)‘"(Tp-Tex,)(V/6)*"=(738672)‘'2x(370.-232.4)(0.2857)‘“=118.23kW/m2 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  7  Date:  5/06/94  Wed.  Time:  11 :30  am 

Pre-Experiment 

Heat  Flux  at  400  mm  Location  (Planned)  (kW/m^):  20 
Slider  Speed  (mm/s):  1.27  mm/s 


Burning  Sample  Data 

Material:  1/4"  PMMA  Sheet  unpainted 

Instrumentation 

Number  of  Thermocouples:  5 


Dimension(LxWxH,  mm):  800x152x6.35 
Number  of  Heat  Flux  Gauges:  2 


Flux  Gauge  Reading  @  400  mm  (Serial  No.  27844,  mV);  5.85-19.93  kW/m2-19.35  kW/m^  @425  mm 

Preheating  Time  without  Pilot  Flame  (seconds):  90 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  27 

Speed  Used  (SlM"Steps"):  200 

Number  of  Steps  (IlM"Steps"):  80000 

File  Names  (.PRN);  P20K506 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec);  800 
Time  to  Start  Moving  the  Sample  (sec): 

Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (“C);  23. 

Observations  HG2  polarity  wrong,  surface  TC2  short  circuit,  surface  TC2  was  painted  10  cm  up  and 
down  streams.  All  TC’s  are  painted.  Note  that  die  file  column  numbers  were  meissfyf  up  at  the 
beginning. 

Personnel  Jeff  Cote,  Y.C. 

Processing  of  Data:  SUR  TCI:  1/t=2.08,  T_=370®C  T  =272  0°C 

kpc=0.2798x  1200  x2200=738672,  >  ext  •  , 

qo=(7tkpc)‘“(Tp-TeJ(V/6)*'2=(73g672)i/2x(370.-272.0)(2.08)*'2=i21.6kW/m2 
SUR  TC2:  1/t  =  1.05,  Tp=370°C,  T,«=272.0®C,  kpc=0.2798x  1200  x2200=738672, 
qS=(7tkpc)^/2(-p^-TeJ(V/6)^"=(738672)‘'2x(370.-272.0)(1.05)‘'2=86.42  kW/m^ 


P20K506 

1  / 4"  PMMA  Sheet,  V=1  .27  mm 


Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  31  Date:  6/29/94  Wed.  Time:8:30pm 

Pre-Experiment 

Heat  Flux  at  400  mm  Location  (Planned)  (kW/m^):  10 
Slider  Speed  (mm/s):  0.50  mm/s 

Burning  Sample  Data 

Material:  1/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  800x152x6.35 

Instrumentation 

Number  of  Thermocouples:  3  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

t 

z  (mm) 

File  Column 

Note 

1 

TC 

N/A 

_ 1 

TC 

201 

IHHHHi 

R 

TP 

_ NZA _ 

7 

TC 

IHHIHH 

C 

7 

HG 

IHiRHHI 

0 

F. 

-3 

TC 

601 

0 

D 

3 

HG 

601 

0 

F 

SERIAL  525842 

Experiment 

Flux  Gauge  Reading  @  400  mm  (Serial  No.  27844,  mV):  3.20-10.22  kW/m2-9.92  kW/m^  @425  mm 

Preheating  Time  without  Pilot  Flame  (seconds):  240 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Hame  (seconds):  30 

Speed  Used  (SlM"Steps"):  79 

Number  of  Steps  (IlM"Steps"):  lOOtKX) 

File  Names  (.PRN):  P10K629 
Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  25(X) 

Time  to  Start  Moving  the  Sample  (sec):  -422.4 
Ignition  Heat  Flux  (Calculated,  kW/m^):  5  3.1 

Ambient  Conditions 
Temperature  (°C):  23. 

Observations  Burning  width- 10  cm  but  is  reduces  somewhat-7.5  cm,  eventually  to  5  cm 
Personnel  Jeff  Cote,  Y.G. 

Processing  of  Data:  SUR  TCI:  l/x  =0.2857,  Tp=370°C,  Te«=228.9“C, 

kpc=0.2798  X 1200  X2200=738672, 

qj;=(7tkpc)‘'2(Tp-Te,,)(V/6)^'2=(738672)i^x(370.-228.9)(0.2857)i"=64.82kW/m=' 

SURTC2:  1/t=0.2874,  Tp=370“C,  T„t=228.9“C, 

qS = (7tkpc)‘"(Tp-Te;tt)(V/6)‘"=(738672)‘^2  x  (370-228.9)(0.2874)*®=65.01  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  39  Date:  8/12/94  Fri  Time:6.30  pm 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  30 
Slider  Speed  (mm/s):  1.90  mm/s 

Burning  Sample  Data 

Material:  3/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  600x152x6.35 


Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Hux  Gauges:  2 


Station 

Sensor 

X  (mm) 

f  (mm) 

File  Column 

Note 

1 

TC 

1.38 

B 

GAS  TCI 

1 

TC 

0 

C 

2 

TC 

350 

1.43 

D 

GAS  TC2  II 

2 

TC 

HKSillHBi 

0 

E 

2 

HG 

350 

HSIHIH 

G 

_ 3_ 

TC 

mmm 

0 

F 

3 

HG 

573 

0 

H 

SERIAL  525842 

Experiment 


Hux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  8.7-29.9976  kW/m^ 

Preheating  Time  without  Pilot  Hame  (seconds):  60 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Hame  (seconds):  0 

Speed  Used  (SlM"Steps"):  300 

Number  of  Stq)s  (IlM"Steps"):  80629 

File  Names  (.PRN):  M30K812 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  121.4 
Ignition  Heat  Flux  (Calculated,  kW/m^:  40.1 
Ambient  Conditions 
Temperature  (*C):  23.5 
Observations  Burning  area  11  cm,  constant 
Personnel  Rich  Kashian,  Y.C. 

Processing  of  Data:  SUR  TCI,  l/r=0.26,  Tp=370*C,  Te^=269.2*C, 

269  2)(0  26)‘^-43  41  ^  o  =(kpc)‘'2(Tp-T,«)(V/5)‘'^=(738672)*'2  x(370- 

m3)Si?3)-=52;wkw5^^  T„=306.3-C,  «=(kpc)“(T,-T.J(V«)'«=C738672)-x(370- 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  30  Date:  6/28/94  Wed.  Time:  8:30  pm 

Pre-Experiment 

Heat  Flux  at  400  mm  Location  (Planned)  (kW/m^):  20 
Slider  Speed  (mm/s):  1.27  mm/s 

Burning  Sample  Data 

Material:  1/4"  PMMA  Sheet  unpainted  Dimension(LxWxH,  mm):  800x152x6.35 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

File  Column 

Note 

1 

TC 

_ NZA _ 

GAS  TCt _ 

1 

TG 

HHHHIIIi 

B 

2 

-JTC _ 

.  N/A 

BUSRHIillliiii 

2 

TC 

379 

C 

- . ^ 

HG 

_ 329 _ 

F 

3 

1 

TC 

0 

D 

3 

HG 

600 

0 

F 

SERIAL  525842 

Experiment 

Hux  Gauge  Reading  @  400  mm  (Serial  No.  27844,  mV):  5.85-19.93  kW/m2-.19.35  kW/m"  @425  mm 

Preheating  Time  without  Pilot  Flame  (seconds):  90 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  27  ( 

Speed  Used  (SlM"Steps"):  200 
Number  of  Steps  (IlM"Steps"):  80000 
File  Names  (.PRN):  P20K628 
Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  -138.80 
Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  ("C):  23. 

Observations  Burning  width  -10  cm 
Personnel  Jeff  Cote,  Y.C. 

Processing  of  Data:  SUR  TCI:  1/t=0.85,  Tp=370°C,  T„,=255.5°C, 

kpc=0.2798  X 1200  X2200=738672, 

qS=(Jtkpc)‘^(Tp-T^(V/6)‘^=(738672)‘'2x(370.-255.5)(0.85)i“=91.02kW/m2 
SURTC2:  1/t=0.86,  Tp=370°C,  T,«=272.0°C, 

qS=(itkpc)*'2(j^-TeJ(V/6)i'2=(738672)‘'2x(370.-255.5)(0.86)‘^=91.41kW/m2 
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Experimental  No.:  44 


Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 


Date:  8/22/94  Wed.  Time:  9:30  pm 


Pre-Experiment 

Heat  Flux  at  400  mm  Location  (Planned)  (kW/m^):  30 
Slider  Speed  (mm/s):  1.90  mm/s 

Burning  Sample  Data 

Material:  1/4"  PMMA  Sheet  urpainted  Dimension(LxWxH,  mm):  800x152x6.35 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

z  (mm) 

File  Colunm 

Note  1 

_ 1 

Tr 

101 

_ L37 _ 

R 

1 

TC 

101 

IRH 

r 

^  . 

Tr 

371 

n 

_ 2 

/cc 

371 

IRHHHIi 

__E 

2 

HG 

371 

_ 0 _ 

G 

_3.  . 

TC 

504 

F 

3 

HG 

594 

0 

H 

SERIAL  525842 

Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  8.70-29.9976  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  42 

Time  to  Sanq)le  Ignition  After  Applying  the  Pilot  Flame  (seconds):  0 

Speed  Used  (SlM"Steps"):  300 

Number  of  Steps  aiM"Steps"):  80629 

FUe  Names  (.PRN):  P30K822 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  -50 
Ignition  Heat  Flux  (Calculated,  kW/m^):  BV  - 
Ambient  Conditions 
Temperature  (°C):  23. 

Observations  Burning  width  -15  cm,  constant 
Personnel  Hui  Fang,  Y.C. 

Processing  of  Data:  SUR  TCI:  1/t=0.78,  Tp=370°C,  Te«=287.8°C, 

kpc=0.2798  X 1200  x2200=738672, 

qS=(7tkpc)‘'2(Tp-TeJ(V/6)‘"=(738672)^^X(370.-287.8)(0.78)*"=80.75kW/m2 
SURTC2:  1/t=0.61,  Tp=370°C,  T^=287.8“C, 

qS=(7ikpc)*"(Tp-TeJ(V/6)‘'2=(738672)i'2x(370.-287.8)(0.61)*"=71.4kW/m2 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  45  Date:  8/23/94  Wed.  Time:  10:30  pm 

Pre-Experiment 

Heat  Flux  at  400  mm  Location  (Planned)  (kW/m^):  30 
Slider  Speed  (mm/s):  1.90 

Burning  Sample  Data 

Material:  1/4"  PMMA  Sheet  unpainted  Dimension(LxWxH,  mm):  800x152x6.35 


Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

z  (mm) 

File  Column 

Note 

1 

TC 

199 

1  12 

JB 

1 

Tr 

199 

0 

r 

2 

Tr 

.  .  378 

1-70  . . 

D 

GA.S  TGI _ II 

--  2 _ 

Tr 

.  378 

_ 0 _ 

F. 

2 

HG 

378 

G 

3 

TC 

598 

iRmHii 

F 

3 

HG 

598 

0 

H 

SERIAL  525842 

Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  8.70-29.9976  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  30 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  0 

Speed  Used  (SlM"Steps"):  300 

Number  of  Steps  aiM"Steps"):  80629 

FUe  Names  (.PRN):  P30K823 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 

Time  to  Start  Moving  the  Sample  (sec):  -45 
Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (°C):  20. 

Observations  Burning  width  -7.5  cm,  constant 
Personnel  Hui  Fang,  Y.C. 

Processing  of  Data:  SUR  TCI:  1/t=0.51,  Tp=370°C,  T,„=273.8“C, 

kpc=0.2798  X 1200  X2200=738672, 

qS = (iikpc)‘"(Tp-T,  J(V/6)'" =(738672)*'^  x  (370.-273 .8)(0.51)‘" =59.09  kW/m" 

SUR  TC2:  1/t=0.8,  Tp=370“C,  Te«=273.8»C, 

qS = (iikpc)‘"(Tp-Te«)(V/6)''2=(738672)‘"  x  (370.-273.8)(0.8)>"=73.99  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  36  Date:  8/09/94  Tu  Time:3:30  pm 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  10 
Slider  Speed  (mm/s):  O.S  mm/s 

Burning  Sample  Data 

Material:  3/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  600x152x6.35 

Listrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

j.  (mm) 

File  Column 

Note 

1 

TC 

174 

1.30 

B 

GAS  TCI 

1 

TC 

174 

0 

C 

2 

TC 

351 

1.43 

D 

2 

TC 

351 

0 

E 

SURFACE  TC2 

2 

HG 

.3_5_1  .. 

0 

G 

SERIAL  86862 

3 

TC 

576 

0 

F 

SURFACE  TC3 

3 

HG 

576 

0 

i 

H 

SERIAL  525842 

Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  2.90-10  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  240 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  130 

Speed  Used  (SlM"Steps"):  79 

Number  of  Stq)s  (IlM"Steps"):  80629 

File  Names  (.PRN):  M10K809 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  582.8 
Ignition  Heat  Flux  (Calculated,  kW/m^: 

Ambient  Conditions 
Temperature  (*C):  23. 

Observations  Burning  area  8.5  cm,  taking  some  data  when  marinite  board  cover  is  on.  Small  flux 
gage  is  recessed  at  the  beginning.  Slide  moves  too  fast,  stops  sometime.  There  is  a  lag 
in  front  of  the  small  flux  gage.  Taking  data  when  cover-on  for  about  82.4  seconds 
Personnel  V.M.,  Y.C. 

Processing  of  Data:  SUR  TCI:  l/r=0.18,  Tp=370'C,  Te«=234.3*C, 

1q(>c=0.2798  X 1200  x2200=738672, 

qS =(Tlvc)i'2(T^.T^(y/5)i/2=(73g672)>« x (370.-234.3)(0. 18)‘'2=49.73  kW/m^ 

SURTC2:  1/t=0.11,  Tp=370*C,  T.«=234.3”C, 

qS=(xkjt)c)‘'2(Tp-T^(V/5)>'2=(738672)‘'2x(370-234.3)(0.)‘'2=39.32kW/m2 

'f 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  37  Date:  8/10/94  Tu  Time:  10:00  pm 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  10 
Slider  Speed  (mm/s):  0.2  mm/s 

Burning  Sample  Data 

Material:  3/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  600x152x6.35 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

(mm) 

File  Colunm 

Note 

1 

TC 

177 
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GAS  TCI 
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TC 
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SERIAL  525842 

Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  2.90-10  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  240 

Hme  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  120 

Speed  Used  (SlM"Stq)s"):  32 

Number  of  Stq)s  (IlM"Steps"):  80629 

File  Names  (.PRN):  M10K810,  M10K81a 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  Ae  Sample  (sec):  486.2 
Ignition  Heat  Flux  (Calculated,  kW/m^:  20.1  kW/m^ 

Ambient  Conditions 
Temperature  (*C):  23. 

Observations  Burning  area  10  cm,  good  ignition,  growing  to  - 12.5  c,  size  little  bit  too  Irage.  Data 
acquization  terminates  about  100  seconds.  —20  cm  when  flame  front  @45  cm  location 
Personnel  Yijun  Hu,  Y.C. 

Processing  of  Data:  k^c=0.2798x  1200  x2200=738672 
SDR  TC2:  l/r=0.14,  Tp=370*C,  T«a=230.4*C, 

qS=(Tkpc)‘^(Tp-T^(V/5)‘«=(738672)‘«x(370-230.4)(0.14)‘'2=44.61  kW/m^ 


Ml  0K81  0 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  38  Date:  8/12/94  Fri  Time:  11.00  am 

Fre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^:  10 
Slider  Speed  (mm/s):  0.35  mm/s 

Burning  Sample  Data 

Material:  3/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  600x152x6.35 

Bistrumentation 

Numba*  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

(mm) 

File  Column 

Note 

1 

TC 

178 

1.38 

B 

GAS  TCI 

1 

TC 
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0 

C 

SURFACE  TCI 

2 

TC 

353 

1.43 
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SERIAL  86862 
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HG 
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SERIAL  525842 

Experiment 

Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  3.0-10.3448kW/m^ 

Prdieating  Time  without  Pilot  Flame  (seconds):  240 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Flame  (seconds):  105 

Speed  Used  (SlM"Steps"):  55 

Number  of  Steps  (IlM"Steps"):  80629 

File  Names  (.PRN):  M10K812 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  452.8 
Ignition  Heat  Flux  (Calculated,  kW/m^):  ^  ^  ^ 

Ambient  Conditions  '■  3 

Tempaature  (*C):  22. 

Observations  Burning  area  12.5  cm,  constant 
Personnel  Rich  Kashian,  Y.C. 

Processing  of  Data:  SUR  TC2,l/r=0.14,  Tp=370*C,  T,,«=227.8‘C, 

kpc=0.2798  X 1200  x2200=738672 

qS  =(irkpc)‘'2(Tp-TeJ(V/a)‘'2=(738672)‘" X (370-227.8)(0. 14)‘'2=46.53  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 
Experimental  No.:  32  Date:  7/8/94  Tu  Time:  11. 30  am 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^:  20 
Slider  Speed  (mm/s):  1.27  mm/s 

Burning  Sample  Data 

Material:  3/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  600x152x6.35 

Instrumentation 

Number  of  Thermocouples:  3  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

x  (mm) 

Z  (mm) 

File  Column 

Note 

1 

TC 

N/A 

GAS  TCI 

1 

TC 

mmm 

0 

B 

SURFACE  TCI 

2 

TC 

N/A 

GASTC2 

2 

TC 
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2 
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SERIAL  86862 

3 

TC 

572 

0 

D 

SURFACE  TC3 

3 

HG 
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E 

SERIAL  525842 

Experiment 


Flux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  5.90-20.5156  kW/m^ 

Preheating  Time  without  Pilot  Flame  (seconds):  90 

Time  to  Sample  Ignition  After  Applying  the  Klot  Flame  (seconds):  10 

Speed  Used  (SlM"Steps"):  200 

Number  of  Steps  (IlM"Stfips"):  100000 

File  Names  (.PRN):  M20K708 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  1500 
Time  to  Start  Moving  the  Sample  (sec):  143.6 
Ignition  Heat  Flux  (Calculated,  kW/m^): 

Ambient  Conditions 
Temperature  (*C):  25 
Observations  Burning  width  ~7.5  cm 
Personnel  V.M.,  Y.C 

Processing  of  Data:  Surface  TCI:  l/T=0.65,Tp=370*C,  Te«=269.6*C, 

kpc=0.2798  X  1200x2200=738672,  ’q  "o  =(kpc)‘'2(Tp-T„t)(l/r*^'^  =(738672)*^^  ^  (370- 
269.6)(0.65)‘'2=69.44  kW/m^ 

Surface  TC2:  l/T=0.73,Tp=370*C,  T^=269.6*C,  q"o=(kpc)‘'2(Tp- 

T,J(l/T‘^'2^(73g672)'/2 X (370-269.6)(0.73)‘'2=73.72  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 


Experimental  No.:  39  Date:  8/12/94  Fri  Time:6.30pm 

Fire-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  30 
Slider  Speed  (mm/s):  1.90  mm/s 

Burning  Sample  Data 

Material:  3/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  600x152x6.35 

Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Hux  Gauges:  2 


Station 

Sensor 

X  (mm) 

^  (mm) 

File  Column 

Note 

1 

TC 

176 

1.38 

B 

GAS  TCI 

1 

TC 

176 

0 

C 

2 

TC 

HISSHHI 

1.43 

D 

GAS  TC2 

2 

TC 

330 

0 

E 

2 

HG 

350 

0 

G 

_ 3 

TC 

573 

0 

1 

F 

3 

HG 

573 

0 

H 

SERIAL  525842 

Experiment 

Hux  Gauge  Reading  ®  425  mm  (Serial  No.  27844,  mV):  8.7-29.9976  kW/m^ 

Preheating  Time  without  Pilot  Hame  (seconds):  60 

Time  to  Sample  Ignition  After  Applying  the  Pilot  Hame  (seconds):  0 

Speed  Used  (SlM"Steps"):  300 

Number  of  Stq)s  (IlM"Stq)s"):  80629 

File  Names  (.PRN):  M30K812 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  121.4 
Ignition  Heat  Hux  (Calculated,  kW/m^:  40.1 
Ambient  Conditions 
TempCTature  (*Q:  23.5 
Observations  Burning  area  11  cm,  constant 
Personnel  Rich  Kashian,  Y.C. 

Processing  of  Data:  SUR  TCI,  1/t=0.26,  T  =370*C,  T«a=269.2*C, 

kpc=0.2798  X  1200x2200=738672,  q"o  =(kpc)*'^(T--Te«)(V/6)‘^2=(738672)*'2  x(370- 

269.2) (0.26)‘'2=43.41  kW/m^ 

SUR  TC2,  l/r=0.83,  Tp=370*C,  Te«=306.3*C,  qS=(kpc)‘^(T.-T^(V/5)‘'2=(738672)‘'2x(370- 

306.3) (0.83)>'2=52.00  kW/m^ 
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Data  Sheet 

Constant  Horizontal  Flame  Spread  Experiment 


Experimental  No.:  40  Date:  8/16/94  Tu  Time:  10:30  am 

Pre-Experiment 

Heat  Flux  at  425  mm  Location  (Planned)  (kW/m^):  30 
Slider  Speed  (mm/s):  1.90  mm/s 

Burning  Sample  Data 

Material:  3/4"  PMMA  Sheet  Dimension(LxWxH,  mm):  600x152x6.35 


Instrumentation 

Number  of  Thermocouples:  5  Number  of  Heat  Flux  Gauges:  2 


Station 

Sensor 

X  (mm) 

^  (mm) 

File  Column 

Note 

1 

TC 

176 

1.38 

B 

GAS  TCI 

1 

TC 

C 

SURFACE  TCI 

2 

TC 

350 

D 

^SMSI^^SEIHIIIII^HIII 
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HG 
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SERIAL  525842 

Pbqperiment 


Hux  Gauge  Reading  @  425  mm  (Serial  No.  27844,  mV):  8.7-29.9976  kW/m^ 

Prdieating  Time  without  Pilot  Flame  (seconds):  60 

Time  to  Sample  Ignition  After  Applying  the  Wot  Flame  (seconds):  0 

Speed  Used  (SlM"Steps"):  300 

NumbCT  of  Stq)s  (IlM"Steps"):  80629 

File  Names  (.PRN):  M30K816 

Set-up  File  Name: 

Sampling  Rate  (Hz):  5  Duration  of  Sampling  (sec):  2500 
Time  to  Start  Moving  the  Sample  (sec):  131:4  “L 
Ignition  Heat  Flux  (Calculated,  kW/m^):  40rl-  ^^  .3 
Ambient  Conditions 
Temperature  (*C):  20. 

Observations  Burning  area  15  cm,  original  leading  edge  is  35.5  mm  (?)  ahead  of  ignition  burner. 
Ignition  burner  is  403  mm  from  the  panel  lower  edge 
Personnel  V.M.,  Y.C. 

Processing  of  Data:  SUR  TCI,  l/r=0.35,  T.=370*C,  Toa=270.9*C, 

kpc  =  0.2798 X  1200x2200  =  738672,  q"o  =(kpc)‘'2(T -Te«)(V/5)‘'2=(738672)*'2x(370- 

270.9) (0.35)‘'2=50.45  kW/m^ 

SUR  TC2,  l/r=0.34,  Tp=370*C,  T^=270.9*C,  qS=(kpc)*'2(T.-T^(V/6)‘'2=(738672)*'2x(370- 

270.9) (0.34)‘'2=49.94  kW/m^ 
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Appendix  M  Tables  of  Char  Depth  Measurements 


This  Appendix  contains  the  measurement  of  char  depth  along  the  particle  board  samples 
for  all  3  external  fluxes  and  3  tests  at  each  level.  Measurements  were  made  at  every  5  cm 
along  the  length.  At  each  of  these  locations,  the  center  left  and  right  measurements  refer 
to  the  center  of  the  ^edm^  along  the  width,  right  and  left  were  at  mid  points  between 
the  edge  and  the  center.  The  spedmen  total  thickness  was  also  measured.  Both 
measurements  were  averaged. 

The  test  numbering  system  is  as  follows; 

Letter 

T:  Particle  board 

P:  PMMA  sheet  (6.35  mm  thick  -  Va”) 

M:  PMMA  Sheet  (19  mm  thick  -  V”) 

10k,  20k,  30k  refer  to  the  nominal  external  flux  level  in  kW/m^. 

Last  3  to  four  digits  are  the  month  and  day  of  the  test. 
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Char  depth  raw  data  for  Particle  Board 

Left,  Midde  and  Right  refer  to  approximate  positions  along  the  width  of  the  test  specimen 
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Board  C  :  T20K622.PRN 
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Board  E :  T30K624.PRN 
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